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INTERFERENCE SUPPRESSION IN CDMA SYSTEMS 

DESCRIPTION 

TECHNICAL FIELD: 

The invention relates to Code-Division Multiple Access (CDMA) communications 
5 systems, which may be terrestrial or satellite systems, and in particular to interference 
suppression in CDMA communications systems. 

BACKGROUND ART: 

Code-Division Multiple Access communications systems are well known. For a 

10 general discussion of such systems, the reader is directed to a paper entitled "Multiuser 
Detection for CDMA Systems" by Duel-Hallen, Holtzman and Zvonar, IEEE Personal 
Communications, pp. 46-58, April 1995. 

In CDMA systems, the signals from different users all use the same bandwidth, 
so each user's signal constitutes noise or interference for the other users. On the uplink 

15 (transmissions from the mobiles) the interference is mainly that from other transmitting 
x mobiles. Power control attempts to maintain the received powers at values that balance 
the interference observed by the various mobiles, but, in many cases, cannot deal 
satisfactorily with excessive interference. Where mobiles with different transmission 
rates are supported within the same cells, the high-rate mobiles manifest strong 

20 interference to the low-rate mobiles. On the downlink (transmission towards the 
mobiles) transmissions from base-stations of other cells as well as strong interference 
from the same base-station to other mobiles may result in strong interference to the 
intended signal. Downlink power control may be imprecise or absent altogether. In all 
these so called near-far problem cases, the transmission quality can be improved, or the 

25 transmitted power reduced, by reducing the interference. In turn, for the same 
transmission quality, the number of calls supported within the cell may be increased, 
resulting in improved spectrum utilization. 

Power control is presently used to minimize the near-far problem, but with 
limited success. It requires a large number of power control updates, typically 800 times 

30 per second, to reduce the power mismatch between the lower-rate and higher-rate users. 
It is desirable to reduce the number of communications involved in such power control 
systems, since they constitute overhead and reduce overall transmission efficiencies. 
Nevertheless, it is expected that future CDMA applications will require even tighter 



power control with twice the number of updates, yet the near-far problem will not be 
completely eliminated. It is preferable to improve the interference suppression without 
increasing the number of transmissions by the power control system. 

Multiuser detectors achieve interference suppression to provide potential benefits 
5 to CDMA systems such as improvement in capacity and reduced precision requirements 
for power control. However, none of these detectors is cost-effective to build with 
significant enough performance advantage over present day systems. For example, the 
complexity of the optimal maximum likelihood sequence detector (MLSD) is exponential 
in the number of interfering signals to be cancelled, which makes its implementation 

10 excessively complex. Alternative suboptimal detectors fall into two groups: linear and 
subtractive. The linear detectors include decorrelators, as disclosed by K.S. Schneider, 
"Optimum detection of code division multiplexed signals", IEEE Trans, on Aerospace 
and Electronic Systems, vol. 15, pp. 181-185, January 1979 and R. Kohno, M. Hatori, 
and H. Imai, "Cancellation techniques of co-channel interference in asynchronous spread 

15 spectrum multiple access systems" , Electronics and Communications in Japan, vol. 66-A, 
no. 5, pp. 20-29, 1983. A disadvantage of such decorrelators is that they cause noise 
enhancement. 

Z. Xie, R.T. Short, and C.K. Rushforth, "A family of suboptimum detectors for 
coherent multiuser communications", IEEE Journal on Selected Areas in 

20 Communications, vol. 8, no. 4, pp. 683-690, May 1990, disclosed the minimum mean 
square error linear (MMSE) detector, but such detectors are sensitive to channel and 
power estimation errors. In both cases, the processing burden still appears to present 
implementation difficulties. 

Subtractive interference cancellation detectors take the form of successive 

25 interference cancellers (SIC), as disclosed by R. Kohno et al. f "Combination of an 
adaptive array antenna and a canceller of interference for direct-sequence 
spread-spectrum multiple-access system", IEEE Journal on Selected Areas in 
Communications, vol. 8, no. 4, pp. 675-682, May 1990, and parallel interference 
cancellers (PIC) as disclosed by M.K. Varanasi and B. Aazhang, "Multistage detection 

30 in asynchronous code-division multiple-access communications", IEEE Trans, on 
Communications, vol. 38, no. 4, pp. 509-519, April 1990, and R. Kohno et aL, 
"Combination of an adaptive array antenna and a canceller of interference for 
direct-sequence spread-spectrum multiple-access system", IEEE Journal on Selected 
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Areas in Communications, vol. 8, no. 4, pp. 675-682, May 1990. Both SIC detectors 
and PIC detectors require multi-stage processing and the interference cancellation 
achieved is limited by the amount of delay or complexity tolerated. These detectors are 
also very sensitive to channel, power and data estimation errors. 
5 One particular subtractive technique was disclosed by Shimon Moshavi in a paper 

entitled "Multi-User Detection for DS-CDMA Communications", IEEE Communications 
Magazine, pp. 124-136, October 1996. Figure 5 of Moshavi's paper shows a subtractive 
interference cancellation (SIC) scheme in which the signal for a particular user is 
extracted in the usual way using a matched filter and then spread again using the same 

10 spreading code for that particular user, i.e., the spreading code used to encode the signal 
at the remote transmitter. The spread-again signal then is subtracted from the signal 
received from the antenna and the resulting signal is applied to the next user's 
despreader. This process is repeated for each successive despreader. Moshavi discloses 
a parallel version that uses similar principles. 

15 A disadvantage of this approach is its sensitivity to the data and power estimates, 

i.e., their accuracy and the sign of the data. A wrong decision will result in the 
interference component being added rather than subtracted, which will have totally the 
wrong effect. 

For more information about these techniques, the reader is directed to a paper by 
20 P. Patel and J. Holtzman entitled "Analysis of a Simple Successive Interference 
Cancellation Scheme in a DS/CDMA System", IEEE Journal on Selected Areas in 
Communications, Vol. 12, No. 5, pp. 796-807, June 1994. 

In a paper entitled "A New Receiver Structure for Asynchronous CDMA: STAR - 
The Spatio-Temporal Array-Receiver", IEEE Transaction on Selected Areas in 
25 Communications, Vol. 16, No. 8, October 1998, S. Affes and P. Mermelstein (two of 
the present inventors), disclosed a technique for improving reception despite near/far 
effects and multi-user interference. In contrast to known systems in which the spread- 
again signal is supplied to the input of the despreader of the channel to be corrected, 
Affes' and Mermelstein's proposed system treated all of the users' signals together and 
30 processed them as a combined noise signal. If the components of the received signal 
from the different users were uncorrelated and all had equal power, or substantially equal 
power, this process would be optimal. In practice, however, there will be significant 
differences between the power levels at which the different users' signals are received 



at the base station antenna. The same applies to the downlink. For example, a data user 
may generate much more power than a voice user simply because of the more dense 
information content of the data signal. Also, imperfect power control will result in 
power differences, i.e., channel variations may result in received powers different from 
5 their intended values, despite the best effort of the power-control process to equalize 
them. 

DISCLOSURE OF INVENTION: 

The present invention addresses the need for improved interference suppression 

10 without the number of transmissions by the power control system being increased, and, 
to this end, provides a receiver for a CDMA communications system which employs 
interference subspace rejection to obtain a substantially null response to interference 
components from selected user stations. Preferably, the receiver also provides a 
substantially unity response for a propagation channel via which a corresponding user's 

15 "desired" signal was received. 

According to one aspect of the invention, there is provided a receiver suitable for 
a base station of a CDMA communications system comprising at least one base station 
(11) having a transmitter and a said receiver and a multiplicity (U) of user stations 
(10\...,10 u ) including a plurality (LP) of user stations served by said at least one base 

20 station, each user station having a transmitter and a receiver for communicating with said 
at least one base station via a corresponding one of a plurality of channels (14 1 ,. ..,14 u ), 
the base station receiver for receiving a signal (X(t)) comprising components 
corresponding to spread signals transmitted by the transmitters of the plurality of user 
stations, each of said spread signals comprising a series of symbols spread using a 

25 spreading code unique to the corresponding user station, said base station receiver 
comprising: 

a plurality (U') of receiver modules (20\...,20 NI ,20 d ) each for deriving from 
successive frames of the received signal (X(t)) estimates of said series of symbols 
of a corresponding one of the user stations, 
30 preprocessing means (18) for deriving from the received signal (X(t)) a series of 

observation matrices (F„) each for use by each of the receiver modules (20) in a 
said frame to derive an estimate of a symbol of a respective one of said series of 
symbols, and 



means (19,44;44/l,44/2) for deriving from each observation matrix a plurality of 
observation vectors (V ; Y ; Z l ..Z NI \ Z*) and applying each of the observation 

n — n-\ — n — n — n 

vectors to a respective one of the plurality of receiver modules (20\. . . ,20 NI ,20 d ); 

each receiver module comprising; 

channel identification means (28) for deriving from one of the observation 
vectors a channel vector estimate (H l // N/ y^ -y 1 \ based upon 
parameter estimates of the channel between the base station receiver and 
the corresponding user station transmitter; 

beamformer means (27 1 , . . . ,27 NI ,27 d ; 47 d ) having coefficient tuning means 
(50) for producing a set of weighting coefficients in dependence upon the 
channel vector estimate, and combining means (51,52) for using the 
weighting coefficients to weight respective ones of the elements of a 
respective one of the observation vectors and combining the weighted 
elements to provide a signal component estimate and 
symbol estimating means (29 1 ,...,29 u , 30 1 ,...,30 u ) for deriving from the 
signal component estimate an estimate (h l b Uy ) of a 
symbol (b b?) transmitted by a corresponding one of the user 
stations (10 1 ,...,10 u ), 
wherein said receiver further comprises means (42,43) responsive to symbol 
estimates ($* 9 ... 5 ^ / ;^ 1 ,g 2 5 g 3 ;g /+1 -") and to channel estimates 
(H\...J-C™\J{ i n _ l ) comprising at least said channel vector 
estimates (H 1 ,... 9 H NI ) for channels (14\...,14 NI ) of a first group (I) of said 
plurality of user stations (10\...,10 NI ) to provide at least one constraint 
matrix (c ) representing interference subspace of components of the received 
signal corresponding to said predetermined group, and in each of one or more 
receiver modules (20A d ) of a second group (D) of said plurality of receiver 
modules, the coefficient tuning means (50 A d ) produces said set of weighting 
coefficients in dependence upon both the constraint matrix (c ) and the channel 
vector estimates (//^ so as to tune said one or more receiver modules (20A d ) 
each towards a substantially null response to that portion of the received signal 
(X(t)) corresponding to said interference subspace. 
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Embodiments of the invention may employ one of several alternative modes of 
implementing interference subspace rejection (ISR), i.e. characterizing the interference 
and building the constraint matrix. In a first embodiment, using a first mode 
conveniently designated ISR-TR, each receiver module in the first group generates its 
5 re-spread signal taking into account the amplitude and sign of the symbol and the channel 
characteristics. The re-spread signals from all of the receiver modules of the first group 
are summed to produce a total realization which is supplied to all of the receiver modules 
in the second group. 

Where each receiver module of the second set uses decision feedback, it further 

10 comprises delay means for delaying each frame/block of the observation vector before 
its application to the beamformer. 

Whereas, in ISR-TR embodiments, just one null constraint is dedicated to the 
sum, in a second embodiment, which uses a second mode conveniently designated 
ISR-R, estimated realisations of all the interferers are used, and a null constraint is 

15 dedicated to each interference vector. In this second embodiment, in each receiver 
module of the first set, the symbols spread by the spreader comprise estimated 
realisations of the symbols of the output signal. Also, the constraint waveforms are not 
summed before forming the constraint matrix. Thus, the receiver module estimates 
separately the contribution to the interference from each unwanted (interfering) user and 

20 cancels it by a dedicated null-constraint in the multi-source spatio-temporal beamformer. 
In most cases, estimation of the interference requires estimates of the past, present and 
future data symbols transmitted from the interferers, in which case the receiver requires 
a maximum delay of one symbol and one processing cycle for the lower-rate or low- 
power users and, at most a single null constraint per interferes 

25 In a third embodiment of the invention which uses a third mode conveniently 

designated ISR-D, i.e. the observation vector/matrix is decomposed over sub- 
channels/fingers of propagation path and the beamformer nulls interference in each of 
the sub-channels, one at a time. In most cases, the maximum number of constraints per 
interferer is equal to the number of sub-channels, i.e. the number of antenna elements 

30 M multiplied by the number of paths P. 

In a fourth embodiment using a fourth mode conveniently designated ISR-H 
because it implements null-responses in beamforming using hypothetical realisations of 
the interference, without any delay, each receiver module of the first group further 
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comprises means for supplying to the spreader possible values of the instant symbols of 
the output signal and the spreader supplies a corresponding plurality of re-spread signals 
to each of the receiver modules of the second group. In each receiver module of the 
second group, the despreader despreads the plurality of re-spread signals and supplies 
5 corresponding despread vectors to the beamformer. This embodiment suppresses any 
sensitivity to data estimation errors and, in most cases, requires a maximum of 3 null 
constraints per interferes 

In a fifth embodiment using a fifth mode conveniently designated ISR-RH because 
it uses the past and present interference symbol estimates, in each receiver module of the 

10 first group, the spreader spreads the symbols of the output signal itself and, in each 
receiver module of the second group, the beamformer then implements null-responses 
over reduced possibilities/hypotheses of the interference realization. Conveniently, 
application of the output of the first despreader to the beamformer will take into account 
the time required for estimation of the interferer's symbol. In most cases, the 

15 beamformer will provide a maximum of 2 null constraints per interferes 

In any of the foregoing embodiments of the invention, the channel identification 
unit may generate the set of channel vector estimates in dependence upon the extracted 
despread data vectors and the user signal component estimate. 

For each of the above- identified modes, the receiver modules may employ either 

20 of two procedures. On the one hand, the receiver module may apply the post-correlation 
observation vector to the channel identification unit but supply the observation matrix 
itself directly to the beamformer, i.e. without despreading it. The constraint matrix then 
would be supplied to the beamformer without despreading. 

Alternatively, each receiver module could supply the post-correlation observation 

25 vector to both the channel identification unit and the beamformer. In this case, the 
receiver module would also despread the constraint matrix before applying it to the 
beamformer. 

Where the reception antenna comprises a plurality of antenna elements, the 
beamformer unit may comprise a spatio-temporal processor, such as a filter which has 
30 coefficients tuned by the estimated interference signals. 

The receiver modules may comprise a first set that are capable of contributing a 
constraint waveform to the constraint matrix and a second set that have a beamformer 
capable of using the constraint matrix to tune the specified null response and unity 
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response. In preferred embodiments, at least some of the plurality of receiver modules 
are members of both the first set and the second set, i.e. they each have means for 
contributing a constraint waveform and a beamformer capable of using the constraint 
matrix. 

5 In practice, the receiver modules assigned to the stronger user signals will usually 

contribute a constraint waveform and the beamformer units of the receiver modules 
assigned to other user signals will be capable of using it. 

The receiver module may comprise an MRC beamformer and an ISR beamformer 
and be adapted to operate in multi-stage, i.e., for each symbol period of frame, it will 

10 carry out a plurality of iterations. In the first iteration, the constraints set generator will 
receive the "past" and "future" estimates from the MRC beamformer and the "past" 
symbol estimate, i.e., from the previous frame, and process them to produce a new 
symbol estimate for the first iteration. In subsequent iterations of the current symbol 
period or frame, the constraints-set generator will use the "future" estimate from the 

15 MRC beamformer, the previous estimate from the ISR beamformer and the symbol 
estimate generated in the previous iteration. The cycle will repeat until the total number 
of iterations have been performed, whereupon the output from the receiver module is the 
desired estimated symbol for the current frame which then is used in the similar 
iterations of the next frame. 

20 The ISR receiver module comprising both an MRC beamformer and an ISR 

beamformer may comprise means {101Q d } for extracting from the ISR beamformer 
(47Q d ) an interference-reduced observation vector and reshaping the latter to produce an 
interference-reduced observation matrix for despreading by the despreader. The channel 
identification unit then uses the despread interference-reduced observation vector to form 

25 interference-reduced channel vector estimates and supplies them to the residual MRC 
beamformer for use in adapting the coefficients thereof. 

The ISR beamformer may process blocks or frames of the observation vector that 
are extended by concatenating a current set of data with one or more previous frames or 
blocks of data. 

30 The different receiver modules may use different sizes of frame. 

In order to receive signals from a user transmitting multicode signals, the ISR 
receiver module may comprise a plurality of ISR beamformers and despreaders, each for 
operating upon a corresponding one of the multiple codes. The channel identification 
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unit then will produce a channel vector estimate common to all of the multicodes, spread 
that channel vector estimate with each of the different multicodes and supply the resulting 
plurality of spread channel vector estimates to respective ones of the plurality of ISR 
beamformers. 

5 The channel identification unit of the multicode ISR receiver module may receive 

its post-correlation observation vector from a despreader (19 d ' 6 ) which uses a compound 
code comprising each of the multicodes weighted by the corresponding symbol estimate 
from a respective one of a corresponding plurality of decision-rule units. The despreader 
will use the compound code to despread the observation matrix and supply the 

10 corresponding compound post-correlation observation vector to the channel identification 
unit. The channel identification unit will use that vector to produce the channel vector 
estimate and spread it using the different ones of the multicodes to produce the spread 
channel vector estimates. 

The ISR receiver module may comprise a despreader 19S d ' F " using 

15 a plurality of codes which comprise segments of a main code specified for that user. 
Each segment corresponds to a symbol, and to a symbol duration in a large block of 
data, the number of segments being determined by the data rate, i.e. , number of symbols 
within a block, of that user. Each receiver module may have a different number of 
segments assigned thereto according to the data rate of the corresponding user. 

20 Embodiments of the invention may be adapted for use in a user/mobile station 

capable of receiving user-bound signals transmitted by a plurality of base stations each 
to a corresponding plurality of users, the receiver then comprising a selection of receiver 
modules each corresponding to a different base station and configured to extract a 
preselected number of said user-bound signals. Where the particular user/mobile station 

25 is included in the preselected number, the receiver module may comprise a similar 
structure to the above-mentioned multicode receiver, the plurality of despreaders being 
adapted to despread the observation matrix using respective ones of a set of codes 
determined as follows: (1) a pre-selected number NB of base stations from which the 
mobile receives signals and which have been selected for cancellation - represented by 

30 index v' which ranges from 1 to NB; (2) a preselected number (1 to NT) of interferers 
per base station preselected for cancellation; (3) the data rates of the selected interferers. 

Thus, according to a second aspect of the invention, there is provided a user 
station receiver for a CDMA communications system comprising a plurality (NB) of base 
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stations (11) and a multiplicity (U) of user stations (10\...,10 u ), at least a plurality (U') 
of the user stations being in a cell associated with one of said base stations and served 
thereby, said one base station having a plurality of transmitter modules for spreading 
user signals for transmission to the plurality (U') of user stations, respectively, and a 
5 receiver for receiving spread user signals transmitted by the plurality (U') of user 
stations, the user stations each having a receiver for receiving the corresponding spread 
user signal transmitted by the base station, said plurality (IT) of user stations each 
having a unique spreading code assigned thereto for use by the user station and the 
corresponding one of the base station transmitter modules to spread the user signals of 
10 that user for transmission, 

the spread user signals transmitted from the base station transmitter modules to 
a particular one of the plurality (U') of user stations propagating via a plurality of 
channels (14 1 ,. ..,14 u ), respectively, 

the receiver of a particular one of said plurality (U') of user stations receiving a 
15 signal (X(t)) comprising components corresponding to spread user signals for said 
particular user station and spread user signals transmitted by other transmitter modules 
of said plurality (NB) of base stations for other users, each of said spread user signals 
comprising a series of symbols spread using the spreading code associated with the 
corresponding one of the user stations, 
20 said user station receiver comprising: 

a plurality (NB) of receiver modules (20 v ) each for deriving from successive 

frames of the received signal (X(t)) estimates of sets of said series of symbols 

from a corresponding one of the base stations, 

preprocessing means (18) for deriving from the received signal (X(t)) a series of 
25 observation matrices (Y n ) each for use by each of the receiver modules (20 v ) in 

a said frame to derive estimates of sets of said symbols, and 
means (19,44) for deriving from each observation matrix a plurality of sets of 
observation vectors (F' U1 ... Y*' M * Fm ' Z 1 '' 11 2?' ,ni,Fni ) and applying each of the 
sets of observation vectors to a respective one of the plurality of receiver modules 
30 (20 v '); 

each receiver modules comprising; 

channel identification means (28T V ) for deriving from the respective one 
of the sets of observation vectors a set of spread channel vector 
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estimates (j^'.u f ( ^ f ni\ based upon parameter estimates of the 
channel between the corresponding one of the base stations and said user 
station; 

beamformer means (47T v '- l > 1 477 v / ** / ' F "') having coefficient tuning 
5 means for producing sets of weighting coefficients in dependence upon the 

sets of channel vector estimates, respectively, and combining means for 
using each of the sets of weighting coefficients to weight respective ones 
of the elements of a respective one of the observation vectors and 
combining the weighted elements to provide a corresponding set of signal 

10 component estimates ' u ,..., /yNI,Fm ) and 

symbol estimating means ( 29T v ' xl 29T vf ^ IyFm ) for deriving from the 
set of signal component estimates a set of estimates ( b v n ,1>1 , $ v ^ yNIyFm ) of 
symbols spread by the corresponding one of the transmitter modules and 
transmitted by the base station; 

15 said user station receiver further comprising means (42,43) responsive to said 

symbol estimates (6 v '' l '\... 9 6 v eKg 2 ,? 3 ) and channel 

f n n n n n 

estimates from each of said plurality (NB) of receiver modules, said 

channel estimates comprising at least channel vector estimates (27 v ') for channels 

n 

(14 v ) between the user station receiver and said base stations, for providing at 
20 least one constraint matrix (c ) representing interference subspace of 

components of the received signal corresponding to said spread signals, and in 
each of said receiver modules (20 v ), the coefficient tuning means produces said 
sets of weighting coefficients in dependence upon both the constraint matrix 
(C ) and the channel vector estimates so as to tune said receiver module (20 v ) 
25 towards a substantially null response to that portion of the received signal (X(t)) 

corresponding to said interference subspace. 

Where the signal destined for the particular user/mobile station is not one of the 
preselected number of signals from the corresponding base station, the receiver may 
further comprise an ISR receiver module which has means for updating the ISR 
30 beamformer coefficients using the channel vector estimates from at least some of the 
receiver modules that have generated such channel vector estimates for the preselected 
signals for the same base station. 
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Where the rates of the different users are not known to the instant mobile station, 
the codes may comprise a fixed number of segments N m which is predetermined as a 
maximum data rate to be received. Any slower rates will effectively be oversampled for 
processing at the higher rate. 
5 The complexity of the multicode embodiments may be reduced by reducing the 

number of codes that are used by the despreaders. In particular, the bank of despreaders 
may use a set of codes that represent summation of the codes of the different NI 
interferers, to form a compound code which reduces the total number of codes being 
used in the despreaders. 
10 According to another aspect of the invention, there is provided a STAR receiver 

comprising an MRC beamformer which operates upon an observation vector which has 
not been despread. 

Of course, that does not preclude having all channels feed their interference 
components to all other channels. 
15 Receivers embodying the present invention can operate in a multiple-input, 

multiple-output (MIMO) system, i.e. with multiple transmit antennas and multiple 
receive antennas. 

The foregoing and other objects, features, aspects and advantages of the present 
invention will become more apparent from the following detailed description, in 
20 conjunction with the accompanying drawings, of preferred embodiments of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS: 

Figure 1 is a schematic diagram illustrating a portion of a CDMA 
communications system comprising a plurality of user stations, typically mobile, and a 
25 base station having a reception antenna comprising an array of antenna elements, and 
illustrating multipath communication between one of the user stations and the array of 
antennas; 

Figure 2 is a simplified schematic diagram representing a model of the part of the 
system illustrated in Figure 1 ; 
30 Figure 3 is a detail block diagram of a spreader portion of one of the user 

stations; 

Figures 4(a) and 4(b) illustrate the relationship between channel characteristics, 
power control and signal power; 
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Figure 5 is a simplified block schematic diagram of a base station receiver 
according to the prior art; 

Figure 6 is a detail block diagram of a preprocessing unit of the receiver; 

Figure 7 is a detail block diagram of a despreader of the receiver; 
5 Figure 8 illustrates several sets of users in a CDMA system ranked according to 

data rate; 

Figure 9 is a detail block diagram showing several modules of a receiver 
embodying the present invention, including one having a beamformer operating on data 
that has not been despread; 
10 Figure 10 is a detail schematic diagram showing a common matrix generator and 

one of a plurality of beamformers coupled in common thereto; 

Figure 1 1 is a block diagram corresponding to Figure 9 but including a module 
having a beamformer operating upon data which has first been despread; 

Figure 12 is a schematic diagram of a user-specific matrix generator and an 
15 associated beamformer of one of the receiver modules of Figure 11; 

Figure 13 is a detail block schematic diagram of a receiver using total realisation 
of the interference to be cancelled (ISR-TR) and without despreading of the data 
processed by the beamformer; 

Figure 14 illustrates a respreader of one of the receiver modules of Figure 13; 
20 Figure 15 is a detail block schematic diagram of a receiver using individual 

realisations of the interference (ISR-R) and without despreading of the data processed by 
the beamformer; 

Figure 16 is a simplified block diagram of a receiver which decomposes each 
realisation of the interference over diversity paths (ISR-D) and without despreading of 
25 the data processed by the beamformer; 

Figure 17 is a simplified schematic block diagram of a receiver employing 
interference subspace rejection based upon hypothetical values of the symbols (ISR-H) 
and without despreading of the data processed by the beamformer; 

Figure 18 illustrates all possible triplets for the hypothetical values; 
30 Figure 19 illustrates bit sequences for generating the hypothetical values; 

Figure 20 is a simplified schematic block diagram of a receiver employing 
interference subspace rejection based upon both hypothetical values of the symbols and 
realisations (ISR-RH) and without despreading of the data processed by the beamformer; 
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Figure 21 is a simplified schematic block diagram of a receiver similar to the 
ISR-TR receiver shown in Figure 13 but in which the beamformer operates upon the data 
that has first been despread; 

Figure 22 is a simplified schematic block diagram of a receiver similar to the 
5 ISR-R receiver shown in Figure 15 but in which the beamformer operates upon data that 
has first been despread; 

Figure 23 is a simplified schematic block diagram of a receiver similar to the 
ISR-D receiver shown in Figure 16 but in which the beamformer operates upon data that 
has first been despread; 
10 Figure 24 is a simplified schematic block diagram of a receiver similar to the 

ISR-H receiver shown in Figure 18 but in which the beamformer operates upon data that 
has first been despread; 

Figure 25 illustrates bit sequences generated in the receiver of Figure 24; 

Figure 26 is a simplified schematic block diagram of a receiver similar to the 
15 ISR-RH receiver shown in Figure 20 but in which the beamformer operates upon data 
that has first been despread; 

Figure 27 illustrates an alternative STAR module which may be used in the 
receiver of Figure 5 or in place of some of the receiver modules in the receivers of 
Figures 13-17, 20-24 and 26; 
20 Figure 28 illustrates a receiver module which both contributes to the constraint 

matrix and uses the constraint matrix to cancel interference (JOINT-ISR); 

Figure 29 illustrates a multi-stage ISR receiver module; 

Figure 30 illustrates successive implementation of ISR; 

Figure 31 illustrates a receiver module which uses ISR to enhance channel 
25 identification; 

Figure 32 illustrates extension of the frame size to reduce noise enhancement and 
facilitate asynchronous operation and processing of high data rates; 

Figure 33 illustrates implementation of ISR with mixed spreading factors; 
Figure 34 illustrates an uplink ISR receiver module for a user employing 
30 multicode signals; 

Figure 35 illustrates a modification of the receiver module of Figure 34; 
Figure 36 illustrates how multirate can be modelled as multicode; 
Figure 37 illustrates frame size determination for multirate signals; 
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Figure 38 illustrates grouping of multirate signals to correspond to a specific 
user's symbol rate; 

Figure 39 illustrates an "uplink" multirate ISR receiver module for a base station; 

Figure 40 illustrates one of a plurality of "downlink" multirate receiver modules 
5 for a user station operating as a "virtual base station"; 

Figure 41 illustrates a "downlink" multirate receiver module of the user station 
of Figure 41 for extracting signals for that user station; 

Figure 42 illustrates a multicode alternative to the receiver module of Figure 40; 

Figure 43 illustrates a second alternative to the receiver module of Figure 40; 
10 Figure 44 illustrates an ISR receiver module using pilot symbols; 

Figure 45 illustrates in more detail an ambiguity estimator of the receiver module 
of Figure 44; 

Figure 46 illustrates an alternative ISR receiver module using pilot channels; 
Figure 47 illustrates an alternative ISR receiver module employing symbol 
15 decoding at an intermediate stage; 

Figure 48 illustrates modelling of the downlink as an uplink; and 
Figure 49 illustrates a transmitter having multiple antennas with which receivers 
embodying the invention can operate. 

20 BEST MODE(S) FOR CARRYING OUT THE INVENTION: 

In the following description, identical or similar items in the different Figures 
have the same reference numerals, in some cases with a suffix. 

The description refers to several published articles. For convenience, the articles 
are cited in full in a numbered list at the end of the description and cited by that number 
25 in the description itself. The contents of these articles are incorporated herein by 
reference and the reader is directed to them for reference. 

Figures 1 and 2 illustrate the uplink of a typical asynchronous cellular CDMA 
system wherein a plurality of mobile stations 10 1 . . . 10 u communicate with a base-station 
1 1 equipped with a receiving antenna comprising an array of several antenna elements 
30 12 1 ...12 M . For clarity of depiction, and to facilitate the following detailed description, 
Figures 1 and 2 illustrate only five of a large number (U) of mobile stations and 
corresponding propagation channels of the typical CDMA system, one for each of a 
corresponding plurality of users. It will be appreciated that the mobile stations 10 1 . . . 10 u 
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will each comprise other circuitry for processing the user input signals, but, for clarity 
of depiction, only the spreaders are shown in Figure 2. The other circuitry will be 
known to those skilled in the art and need not be described here. Referring to Figure 
2, the mobile stations 10 1 . . . 10 u comprise spreaders 13 1 . . . 13 u , respectively, which spread 
5 a plurality of digital signals b* ... b„ of a corresponding plurality of users, respectively, 
all to the same bandwidth, using spreading codes c J (t). ..c u (t), respectively. The mobile 
stations lO^.-lO" transmit the resulting user signals to the base station 11 via channels 
14\.A4 V 9 respectively, using a suitable modulation scheme, such as differential binary 
phase shift keying (DBPSK). Each of the mobile stations lO 1 ...^ receives commands 

10 from the base station 11 which monitors the total received power, i.e. the product of 
transmitted power and that user's code and attenuation for the associated channel and 
uses the information to apply power control to the corresponding signals to compensate 
for the attenuation of the channel. This is represented in Figure 2 by multipliers 
lS^-.lS 11 which multiply the spread signals by adjustment factors (t)...^ (t), 

15 respectively. The array of M omni-directional antenna elements 12 1 ...12 M at the base 
station 1 1 each receive all of the spread signals in common. The channels 14 1 . . . 14 u have 
different response characteristics H 2 (t)...H u (t), respectively, as illustrated in more detail 
in Figure 1, for only one of the channels, designated channel 14 u . Hence, channel 14 u 
represents communication via as many as P paths between the single antenna of the 

20 associated mobile station 10 u and each of the base station antenna elements 12 1 ...12 M . 
The other channels are similarly multipath. 

As before, it is presumed that the base station knows the spreading codes of all 
of the mobile stations with which it communicates. The mobile stations will have similar 
configurations so only one will be described. Thus, the mobile station 10 u first 

25 differentially encodes its user's binary phase shift keyed (BPSK) bit sequence at the rate 
1/T, where T is the bit duration, using circuitry (not shown) that is well-known to 
persons skilled in this art. As illustrated in Figure 3, its spreader 13 u then spreads the 
resulting differential binary phase shift keyed (DBPSK) sequence b" (or b u (t) in the 
continuous time domain as represented in Figure 3) by a periodic personal code sequence 

30 C[ u (or cf(t) in the continuous time domain) at a rate 1/T C , where T c is the chip pulse 
duration. The processing gain is given by L = T/T C . For convenience, it is assumed that 
short codes are used, with the period of <?{t) equal to the bit duration T, though the 
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system could employ, long codes, as will be discussed later, with other applications and 
assumptions. Over one period T, the spreading code can be written as: 

c\i) = £ c?<t>(t-lT c ), (1) 

/=0 

5 where c" = ±1 for 1 = 0, L - 7, is a random sequence of length L and <£(f) is the 
chip pulse as illustrated in Figure 3. Also, with a multipath fading environment with P 
resolvable paths, the delay spread at is small compared to the bit duration (i.e. at < 

As illustrated in Figures 4(a) and 4(b), following signal weighting by the power 

10 control factor i/£ c (r) 2 , the spread signal is transmitted to the base station 11 via channel 
14 u . Figure 4(a) shows the "real" situation where the channel characteristics comprise 
a normalized value W(T) and a normalization factor \p u ch (T) which relates to the 
"amplitude" or attenuation of the channel, i.e. its square would be proportional to the 
power divided by the transmitted power. In Figure 4(a), power control is represented 

15 by a multiplier 17 u , and the subscript "pc". Figure 4(b) shows that, for convenience, the 
channel characteristics can be represented (theoretically) by the normalized value W(t) 
and the normalization factor \j/ u ch {t) included in a single power factor \p u (t) which is equal 
to ^pciO^kchiO- ^pdO IS ^ e f actor by which the transmitted signal is amplified or 
attenuated to compensate for channel power gain in \j/ ch (t) and to maintain the received 

20 power (\J/ U (t)) 2 at the required level. 

In such a CDMA system, the signal of each of the mobile stations lO 1 ...!© 0 
constitutes interference for the signals of the other mobile stations. For various reasons, 
some of the mobile stations will generate more interference than others. The components 
of one of these "strongly interfering" user stations and its associated channel are 

25 identified in Figures 1 and 2 by the index "i". The components of one of the other 
"low-power" user stations and its associated channel also are illustrated, and identified 
by the index "d". The significance of this grouping of "interfering" and "low-power" 
user stations will be explained later. 

At the base station 11, the spread data vector signals X 1 (t)...X u (t) from the base 

30 station antenna elements 12 1 ...12 M , respectively, are received simultaneously, as 
indicated by the adder 16 (Figure 2), and the resulting observation vector X(t) is supplied 
to the receiver (see Figure 5). The sum of the spread data vectors (signals) X ] (t).. .^(t) 
will be subject to thermal noise. This is illustrated by the addition of a noise signal 
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component N th (t) by adder 16. The noise signal N th (t) comprises a vector, elements of 
which correspond to the noise received by the different antenna elements. 

Figure 5 illustrates a spatio-temporal array receiver (STAR) for receiving the 
signal X(t) at the base station 11. Such a receiver was described generally by two of the 
5 present inventors in reference [13]. The receiver comprises means, namely a 
preprocessing unit 18 and a plurality of despreaders 19 1 ... 19 u , for deriving observation 
vectors from the observation matrix, and a plurality of spatio-temporal receiver (STAR) 
modules 20*...20 u , each having its input connected to the output of a respective one of 
the despreaders 19 1 ... 19 u . As shown in Figure 6, the preprocessing unit 18 comprises 

10 a matched filter 22, a sampler 23 and a buffer 24. Matched filter 22 convolves the 
antenna array signal vector X(t), which is an M X 1 vector, with a matched pulse <f>(T c - 
t) to produce the matched filtered signal vector Y(t) which then is sampled by sampler 
23 at the chip rate 1/T C , element by element. The sampler 23 supplies the resulting M 
X 1 vectors Y nl , at the chip rate, to buffer 24 which buffers them to produce an 

15 observation matrix Y n of dimension M x (2L-1). It should be noted that, although the 
present inventors' Canadian patent application No. 2,293,097 and United States 
Provisional application No. 60/171,604 had a duplicate of this preprocessing unit 18 in 
each of the despreaders 19 1 ...19 u , it is preferable to avoid such duplication and use a 
single preprocessor 18 to preprocess the received antenna array signal vector X(t). 

20 The despreaders 19 1 ...19 u each have the same structure, so only one will be 

described in detail with reference to Figure 7 which illustrates despreader 19 u . Thus, 
despreader 19 u comprises a filter 25 u and a vector reshaper 26 u . The observation matrix 
Y n is filtered by filter 21" using the pseudo-random number sequence c " corresponding 
to that used in the spreader 13 u of the transmitter, i.e. C/ w , to produce the postcorrelation 

25 observation matrix z" for user u. Vector reshaper 26 u concatenates the columns of the 
M x L matrix z" to form a post-correlation observation vector z u of dimension ML 
X 1. It should be noted that the vector reshaper 26 u need not be a distinct physical 
element but is depicted as such to represent a mathematical function. In practice, the 
function will likely be determined merely by allocation of resources, such as memory. 

30 Referring again to Figure 5 , the post-correlation observation vectors z 1 ...Z u from 

n n 

despreaders 19 1 ...19 u are processed by the STAR modules 20 1 ...20 u , respectively, to 
produce symbol estimates b l n ...b^ corresponding to the transmitted 
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symbols b l n ...b„ (see Figure 2) and power estimates (i// J J 2 _(^) 2 which are supplied 
to subsequent stages (not shown) of the receiver for processing in known manner. 

The STAR modules 20 1 ...20 u each comprise the same elements, so the 
construction and operation of only one of them, STAR module 20 u , will now be 
5 described. 

The STAR module 20 u comprises a beamformer 27 u , a channel identification unit 
28 u , a decision rule unit 29 u and a power estimation unit 30 u . The channel identification 
unit 28 u is connected to the input and output, respectively, of the beamformer 27 u to 
receive the post-correlation observation vector z" and the signal component 

n 

10 estimate ^ respectively. The channel identification unit 28 u replicates, for each frame 
M X L the characteristics H"(t) 9 in space and time, of the associated user's transmission 
channel 14 u . More specifically, it uses the post-correlation observation vector z u and 

n 

signal component estimate to derive a set of parameter estimates h u 9 which it uses 

n 

to update the weighting coefficients w 1 of the beamformer 27 u in succeeding symbol 

n 

15 periods. The symbol period corresponds to the data frame ofMxL elements. 

The beamformer 27 u comprises a spatio-temporal maximum ratio combining 
(MRC) filter which filters the space-time vector z u to produce the despread signal 

n 

component estimate which it supplies to both the decision rule unit 29 u and the 

power estimation unit 30 u . The decision-rule unit 29 u outputs a binary 
20 symbol $^ according to the sign of the signal component estimate The binary 

output signal constitutes the output of the decision rule unit 30 u and is an estimate of the 

corresponding user signal b% spread by spreader 13 u of the corresponding user station 

10 u (Figures 1 and 2). 

The signal component estimate $" is processed in subsequent parts of the 
25 receiver. For example, it may be differentially decoded and, possibly, deinterleaved and 

then data decoded— if the corresponding inverse operations were done before 

transmission. 

The power estimation unit 30 u uses the raw signal component estimate to 
derive an estimate ( xjfj 2 of the power in that user' s signal component of the antenna 
30 array signal vector X(t) and supplies the power estimate ($jj 2 to the subsequent stages 
(not shown) of the receiver for derivation of power level adjustment signals in known 
manner. 
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The receiver shown in Figure 5 will perform satisfactorily if there are no strong 
interferers, i.e., if it can be assumed that all users transmit with the same modulation and 
at the same rate, and that the base-station knows all the spreading codes of the terminals 
with which it is communicating. On that basis, operation of the receiver will be 
5 described with reference to the user channel identified by index u. 

At time t, the antenna array signal vector X(t) received by the elements 12 1 . . . 12 M 
of the antenna array of the one particular cell shown in Figures 1 and 2 can be written 
as follows: 

u 

X(t) = £ X"(r) (2) 
10 " = 1 

where U is the total number of mobile stations whose signals are received at the 
base-station 11 from inside or outside the cell, X u (t) is the received signal vector from 
the mobile station 10 u , i.e., of index w, and N h (t) is the thermal noise received at the M 
antenna elements. The contribution X*(t) of the w-th mobile station 10 u to the antenna- 
15 array signal vector X(t) is given by: 

X u (t) = xlr u (t)H u (i) <8> c u (t)b u (t) 

p (2?l\ 
= # u (t) £ G p u (t) e" p (t) c u (t - T u p (t)) b %t - x u p (t)) 

where IT(t) is the channel response vector of the channel 14 u between the w-th mobile 
20 station 10 u and the array of antenna elements and <8> denotes time-convolution. In the 
right-hand term of the above equation, the propagation time-delays T u p (t) e [0,T] along 
the P paths, p — 1, P, (see Figure 1), are chip-asynchronous, G p (t) are the 
propagation vectors and e u p (t) 2 are the power fractions along each path 
(i.e., ! e p(0 2 = 1) of the total power \p u (t) 2 received from the w-th mobile station 
25 10 u . The received power is affected by path-loss, Rayleigh fading and shadowing. It 
is assumed that G p u (t) , e u p (t) 2 and 0 w (r) 2 vary slowly and are constant over the bit 
duration T. 

In the preprocessing unit 18 (see Figure 6), the antenna array signal vector X(t) 
is filtered with the matched pulse to provide the matched-filtering signal vector Y n (t) for 
30 frame n as follows: 

Y n (t) = — f X(aT/2^nT^t^t / )4>(t / )dt / (3) 
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where /) denotes the temporal support of ^(t) and a G {0,1} stands for a possible 
time-shift by 772 to avoid, if necessary, the frame edges lying in the middle of the delay 
spread (see reference [13]). For the sake of simplicity, it is assumed in the following 
that a = 0. Note that for a rectangular pulse /> is [0,TJ. In practice, it is the 
5 temporal support of a truncated square-root raised- cosine. 

It should be noted that the above description is baseband, without loss of 
generality. Both the carrier frequency modulation and demodulation steps can be 
embedded in the chip pulse-shaping and matched-filtering operations of Equations (1) and 
(3), respectively. 

10 Thus, after sampling at the chip rate 1/T C and framing over 2L - 1 chip samples 

at the bit rate to form a frame, the preprocessing unit 18 derives the M x (2L -1) 
matched-filtering observation matrix: 

where Y nl = Y n (lT c ). 

15 In the despreader 19" (see Figure 7), the post-correlation vector for frame number 

n for user number u is obtained as: 

Ki = ji: Y nMk c;. (5) 

^ k=0 

20 Framing this vector over L chip samples at the bit rate forms the post-correlation 
observation matrix: 

Z n = Z n o , Z n j,..., Z n . (6) 
The post-correlation data model (PCM) (see reference [13]) details the structure of this 
matrix as follows: 

25 K = HX U - K CM ^ (7) 

where z" is the spatio-temporal observation matrix, n u n is the spatio-temporal 
propagation matrix, s " = b"y}/ u n is the signal component and Np CMn is the s P atio " 
temporal noise matrix. Equation 7 provides an instantaneous mixture model at the bit 
rate where the signal subspace is one-dimensional in the M x L matrix space. For 
30 convenience, the vector reshaper 26 u of despreader 19 u transforms the 
matrices z tt n9 H u n a n d N u PCMtn into ( M x Lj - dimensional 
vectors Z M , H u and z u respectively, by concatenating their columns into one 

—n — n —PCM,n 
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spatio-temporal column vector to yield the following narrowband form of the PCM 
model (see reference [13]): 

Z u = H u s" + N u (8) 

~n n n —pCM,n 

To avoid the ambiguity due to a multiplicative factor between h u and s u , the norm 
5 of H u is fixed to Jm 

— n v 

The PCM model significantly reduces inter-symbol interference. It represents an 
instantaneous mixture model of a narrowband source in a one-dimensional signal 
subspace and enables exploitation of low complexity narrowband processing methods 
after despreading. Processing after despreading exploits the processing gain to reduce 

10 the interference and to ease its cancellation in subsequent steps by facilitating estimation 
of channel parameters. 

As discussed in reference [13], the spatio-temporal array-receiver (STAR) can be 
used to detect each user separately at the base-station 11. In addition to exploiting the 
processing gain to reduce interference, the STAR allows accurate synchronization and 

15 tracking of the multipath delays and components and shows inherent robustness to 
interference. The STAR also allows coherent combining of the data. This receiver is 
found to provide fast and accurate time- varying multipath acquisition and tracking. 
Moreover, it significantly improves call capacity by spatio-temporal maximum ratio 
combining (MRC) in a coherent detection scheme implemented without a pilot signal. 

20 For the sake of clarity, the steps of STAR that are relevant to the implementation of the 
present invention will be reviewed briefly below, with reference to receiver module 20 u 
of Figure 5 . 

As shown in Figure 5, the despreader 19 u supplies the post-correlation observation 
vector Z", to both the channel identification unit 28 u and the MRC beamformer 27 u of 
25 STAR module 20 u . Using spatio-temporal matched filtering (W u = H U IM) (i e- 

n — n 

spatio-temporal maximum ratio combining, W U "H U = 1) the STAR module 20 u 
provides estimates of signal component s ? 9 its DBPSK bit sequence b„ and its total 
received power (^/J) 2 as follows: 



30 §* = Real [W^Z 1 ^ } = Real < 



W Z" 



M 



(10) 
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(0 2 = a-«)(Ci) 2 + « KT, (ii) 

where a is a smoothing factor. It should be noted that with ad hoc modifications, 
differential modulation and quasi-coherent differential decoding still apply with DMPSK. 
5 Orthogonal modulation can even be detected coherently by STAR without a pilot 
(references [17] and [18]). Using the post-correlation observation vector Z" and the new 

n 

signal component estimate § u n from the beamformer 27 u , the channel identification unit 
28" provides an estimate H u of the channel 14 u for user station 10 u . The channel 

— n 

identification unit 28 u updates the channel vector estimate h u by means of a decision 

— n 

10 feedback identification (DFI) scheme whereby the signal component estimate § tt is fed 
back as a reference signal in the following eigen-subspace tracking procedure: 

H u =H U + JZ"-HX)S:, (12) 

15 where \i is an adaptation step-size. Alternatively, the product b u could be fed back 
instead of the signal component estimate s u n . It should be noted that, if the modulation 
is complex, the second occurrence of signal component estimate s" should be replaced 
by its conjugate (S ")* This DFI scheme allows a 3 dB coherent detection gain in noise 
reduction by recovering the channel phase offsets within a sign ambiguity without a pilot. 

20 Note that a reduced-power pilot can be used to avoid differential coding and decoding 
(reference [21]). The procedure that further enhances the channel vector 
estimate h u j to obtain f£ u i from the knowledge of its spatio-temporal structure (i.e. 
manifold) allows a fast and accurate estimation of the multipath time-delays 
f" fl in both the acquisition and the tracking modes (both versions of this 

25 procedure can be found in reference [13]). This improved estimation accuracy achieves 
robustness to channel estimation errors, and reduces sensitivity to timing errors, when 
STAR is used in multiuser operation. 

For further information about STAR, the reader is directed to the articles by 
Affes and Mermelstein identified as references [13] and [17] to [21]. 

30 If, as was assumed in reference [13], the spatio-temporal noise vector N u is 

— PCM,n 

spatially uncorrected, power control on the uplink is generally able to equalize the 
received signal powers. However, the assumption that noise is uncorrelated becomes 
untenable on the downlink due to path-loss and shadowing and when the power of 
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particular users {e.g., "priority links", acquisition, higher-order modulations or higher 
data-rates in mixed-rate traffic) is increased intentionally. Within a particular cell, there 
may be users having many different "strengths", perhaps because of different data rates. 
Figure 8 illustrates, as an example, a cell in which there are four different sets of users 
5 arranged hierarchically according to data rate. The first set I comprises users which 
have relatively high data rates, the second set Ml and third set M2 both comprise users 
which have intermediate data rates, and the fourth set D comprises users which have 
relatively low data rates. In practice, the receivers of the high data rate users of set I 
will not need to cancel any "outset" interference from the users in sets Ml, M2 and D, 

10 but their transmissions will contribute to interference for the receiver modules in those 
sets. Intermediate data rate users in sets Ml and M2 will need to cancel "outset" 
interference from the high data rate users of the set I but not from the users in set D. 
They will themselves be contributors of "outset" interference to the users in set D. The 
receivers of users in set D must cancel "outset" interference from sets I, Ml and M2. 

15 It is also possible for a receiver of a user within a particular set to cancel "inset" 

interference from one or more users within the same set; and itself be a contributor to 
such "inset" interference. Embodiments of the invention applicable to these "outset" and 
"inset" situations will be described hereinafter. In the description, where a particular 
user's signal is treated as interference and cancelled, it will be deemed to be a 

20 "contributor" and, where a particular user's receiver module receives information to 
enable it to cancel another user's interference, it will be deemed to be a "recipient". To 
simplify the description of the preferred embodiments described herein, it will be 
assumed that all users employ the same modulation at the same rate. For the purpose 
of developing the theory of operation, initially it will be assumed that, among the mobile 

25 stations in the cell, there will be a first set I of "strong" contributor users, one of which 
is identified in Figures 1 and 2 by index "i", whose received signal powers are relatively 
high and hence likely to cause more interference, and a second set D of "low-power" 
recipient users, one of which is identified in Figures 1 and 2 by index "d", whose 
received signal powers are relatively low and whose reception may be degraded by 

30 interference from the signals from the strong users. In order to receive the low-power 
users adequately, it usually is desirable to substantially eliminate the interference 
produced by the high-power users. For simplicity, most of the preferred embodiments 
of the invention will be described on the basis that the high-power users can be received 
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adequately without interference suppression. It should be appreciated, however, that the 
"strong" user stations could interfere with each other, in which case one could also apply 
to any interfering mobile the coloured noise model below and the near-far resistant 
solution proposed for the low-power user, as will be described later. 
5 Assuming the presence of NI interfering users assigned the indices / = 1 to NI, 

then the spatio-temporal observation vector of any interfering user (u = i €= {1,...,NI}) 
is given from Equation 8 by: 

Z' =#'V +AT 1 " , (13) 

n n n ~~PCM,n 

10 where w can still be assumed to be an uncorrelated white noise vector if the 

— PCM,n 

processing gain of this user is not very low. On the other hand, from the point of view 
of any low-power user (u=d £ {1,...,M}), the spatio-temporal observation vector is: 
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NI 

Z d = H d s* + I d + N d = H d s d n + Y I d * + N d , (14) 

~n — n n ~PCM,n — PCM,n — n n ^ ~PCM,n ~PCM,n 

i = l 



where, in addition to the uncorrelated white noise vector N d there is included a total 

— PCM,n 9 

interference vector j d which sums a random coloured spatio-temporal interference 

-PCM,n 

vector from each interfering mobile denoted by for / = 7, NL At frame 

~PCM,n 

number n, the realization of the vector results from matched-pulse filtering, chip- 

PCM,n 

20 rate sampling, despreading with c f, bit-rate framing, and matrix/ vector reshaping of 
the received signal vector X(t) from the i-th interfering mobile using Equations 
(3) to (6). 

The receiver shown in Figure 5 would receive the signals from all of the user 
stations independently of each other. It should be noted that there is no cross-connection 

25 between the receiver modules 20 1 ,...,20 u , specifically between their STAR modules 
20 1 . . .20 u . . .20 u , for suppression of interference from the signals of mobile stations which 
constitute strong interferers. While the matched beamformer of Equation (9) is optimal 
in uncorrelated white noise, it is suboptimal when receiving the low-power users due to 
spatial-temporal correlation of the interference terms. To allow the accommodation of 

30 additional users in the presence of much stronger interfering mobiles in the target cell, 
in embodiments of the present invention the receiver of Figure 5 is upgraded to obtain 
much stronger near-far resistance, specifically by adapting the beamformer of Equation 
(9) to reject the interference contributions from the interfering strong users. 



26 

In the general case, the total interference i d experienced by a user d in set 

-PCM,n 

D is an unknown random vector which lies at any moment in an interference subspace 
spanned by a matrix, say C PCM , n (i.e., e Vec {C PCM J) with dimension 

depending on the number of interference parameters {i.e., power, data, multipath 
5 components and delays) assumed unknown or estimated a priori. As will become 
apparent from the following descriptions of preferred embodiments, in practice, the 
matrix C PCM ^ which wi H be referred to as the "constraint matrix", can be derived and 
estimated in different ways. To achieve near-far resistance, the beamformer must 
conform to the following theoretical constraints: 

10 



W d "H d = 1 , 

n — n 



W d "H d = 1, 

~~ " ~ n (15} 

w d "r d - o 

IL n I PCM,n ~ U • 



The first constraint provides a substantially distortionless response to the low- 
power user while the second instantaneously rejects the interference subspace and thereby 
15 substantially cancels the total interference. This modification of the beamforming step 
of STAR will be referred to as interference subspace rejection (ISR). 

With an estimate of the constraint matrix C d pCMn available (as described later), 
the ISR combiner {i.e., the constrained spatio-temporal beamformer) \y d after 

n 

despreading is obtained by: 



20 



QpCM,n ~ [CpCM,nCpCM,^) ' 



n = i - r d n d r d " (17) 

u d ti d 

25 W i = ^ (18) 

IL 11 PCM, n H n 
n n 



PCM.n 



where I M * L denotes aM*LxM*L identity matrix. First, the projector II 
orthogonal to the constraint matrix w is formed. It should be noted from Equations 

PCM,n 

(16) and (17) that the inverse matrix 0jL„ is not the direct inverse of constraint 
30 matrix C d PCMn but part of the pseudo-inverse of C PCMn . For convenience, however, 
it will be referred to as the inverse matrix hereafter. Second, the estimate of the low- 
power response vector y 4 is projected and normalized. 
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Whereas, using the above constraints, the ISR beamformer may process the low- 
power user's data vector after it has been despread, it is possible, and preferable, to 
process the data vector without first despreading it. In either case, however, the data 
vector will still be despread for use by the channel identification unit. Although it is 
5 computationally more advantageous to do so without despreading, embodiments of both 
alternatives will be described. First, however, the spread data model of Equation (2) 
will be reformulated and developed and then used to derive various modes that 
implement ISR combining of the data, without despreading, suitable for different 
complementary situations. 

10 

Data Model Without Despreading 

The observation matrix Y n of Equation (4) which provides the post-correlation 
matrix Z n of Equation (7) by despreading and framing at the bit rate, can be expressed 
as: 

15 y = £ Y u n + Nf, (19) 

where each user u contributes its user-otservation matrix y u n , obtained by Equations 
(3) and (4) with X(t) replaced by X"(t) in Equation (3), and where the preprocessed 
thermal noise contributes: 

Nf = [N&(nT),N* h (nT + TJ^N^nT + (2L - 2)T C )\ (20) 
20 Using the fact that any bit-triplet b„_ u b^b" +l contributing to channel 

convolution (see Equation (2a) in y u n can be composed as: 

~K-M,bi\ = [i,o,o] + K [o,i,0] + z>; + i [o,o,i], (21) 

the sequence b u (t) can be locally approximated over the n-th block by means of the 
canonic generating sequences g\t), £ 2 (r) and g 3 (t) in Figure 25 as: 

25 b\t) = Kg_Ht) + K_^{t) + ^1i«H0, (22) 

where the indices L / t / G |1 , 2 3 } are permuted at each block so that the 
corresponding canonic generating sequences locally coincide with [0, 1, 0], [1, 0, 0] and 
fO, 0, 1], respectively. Assuming slow time- variations of \p(t) and H(t) compared to the 
symbol duration: 

30 k = s:n„ + c.n, n - s: + x +Un , (23) 

where the canonic user-observation matrices y^ n are obtained by Equations (3) and (4) 
with X(t) in Equation (3) replaced, respectively for k = -1, 0, -hi, by: 

X k \t) = H\t) ® gH0C(t). (24) 



28 

Good approximations of y , and y can be actually obtained at each iteration by 
L simple backward/forward shifts of the columns of y£„ with zero column inputs. 

It should be noted that the canonic generating sequences allow more accurate 
reconstruction (e.g., overlap-add) of time-varying channels. Also, the resulting 
5 decomposition in Equation (23) holds for long PN codes. 

It should be noted that this decomposition also holds for any complex-valued 
symbol-triplet b n u . u b* 9 b" +x . With ad hoc modifications, therefore, the ISR approach 
according to this invention applies to any complex modulation (e.g., MPSK, MQAM, 
even analog). This new signal decomposition is used to derive the different 
10 implementations of ISR which will be described later. 

With respect to the low-power user assigned the index d and the NI strong 
interfering mobiles assigned the indices / = 1, .„, NI, the observation vector obtained by 
reshaping the observation matrix, before despreading, can now be rewritten as: 



15 



20 



Y = Y d + I d + I + N , (25) 
where the first canonic observation vector y d appears as the "channel" vector of the 

— 0,n 

low-power user d. The total interference vector before despreading: 

NI NI NI 

i = l i = l f = l 7 

is the sum of the interfering signal vectors y i and: 



r -k n = s °-^ hn - (27) 

25 is the intersymbol interference (ISI) vector of user u. In large processing gain situations, 
the self ISI vector i d can be combined with the uncorrelated spatio-temporal noise 

~ISI,n 

vector n , leading to the following data vector model before despreading: 

n 

Y = Yi s d n + / + N . (28) 
Despreading the observation vector in the above equation with the spreading 
30 sequence of the low-power user d provides the data vector model after despreading in 
Equation (14). It is possible to derive a finer decomposition of the data model to allow 
implementation of one or more of the ISR modes over diversities. 



29 

Finer Decomposition of the Data Model Over Diversities 

Thus, Equation (2a) can be further decomposed over the N f = MP diversity 
branches or fingers in such a way that the observation signal contribution l? J (t) received 
by the m-th antenna along the /?-th path for / = (p - 1)M + m = 7,..., N f can be 
5 separated as follows: 



x\t) = j2 xu/ (0- (29) 

The observation signal contribution from the/-th finger is defined as: 

10 * tt/ (0 - r(W»J{t) ® c%t)b\t) 

= r(t)G;\t)s u p {t)b\t - r u p {t))c\t - t*(o), 1 ; 

where the propagation vector from the/-th finger is: 

G u /(t) = y u f (t)R m . (31) 

15 In the above equation, the scalar y u f (t) is the channel coefficient over the/-th finger and 
En = [0,... f 0,l,0,...,0] T is a M x 1 vector with null components except for the m-th 
one. With the above definitions, one can easily check the following decompositions of 
the channel and the propagation vectors: 



20 H\t) = (32) 



M 

G p \t) = £ G^'^it). (33) 



m = \ 



Accordingly, after preprocessing, the matched-filtering observation matrix can be 
25 decomposed as follows: 

u=l u=\ f=\ 

where each user u contributes its user-observation matrices from fingers / = 
30 l,...,N f , obtained by Equations (3) and (4) with X(t) replaced by X* J (t) in Equation (3). 
Note that the complex channel coefficient $(nT) = yf(nT)e u p (nT) is separated from the 
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matrix 1 Y"/ which contains a purely-delayed replica of the spread-data without 
attenuation or phase offset from finger /. This matrix, which is obtained by Equations 
(3) and (4) with X(t) in Equation (3) replaced by: 



X»S(t) = J? m 5(r- r p (t)) ® b u (t)c"(t), (35) 



can be further decomposed over the canonic generating sequences as follows: 



rf = KYI + blX-in + blXin, 06) 

10 

where the canonic user-observation matrices y"/ n from finger / are obtained by 
Equations (3) and (4) with X(t) in Equation (3) replaced, respectively for k = -1, 0, +1, 
by: 

15 X u /{t) = RJ(t- Tp (t)) ® gHOcV), (37) 

where 8(t) denotes the Dirac impulse. Therefore one obtains: 



Y n = EEE C4,yZ + yf . ( 3 «) 

«=1 /=! jfc=-l 



20 

A coarser decomposition over fingers of the total interference vector before despreading 
defined in Equation (26) gives: 

25 /=1 /=1 /=1 

After despreading with the spreading sequence of the low-power user d, it gives: 

NI NI N f 

i* = y i<k = = y y >it n &j*M . (40) 

i=\ i=\ f=\ 



30 1 This matrix is real-valued in the case of a binary modulation. 
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Embodiments of the invention which use the above decompositions of interference, 
denoted as ISR-D implementations before and after despreading, will be described later 
with reference to Figures 16 and 23. 



5 ISR Combining Before Despreading 

As described hereinbefore, the combining step of STAR is implemented without 
despreading by replacing Equation (9) for the low-power user with: 

s d n = Real\w d n "Y^ (41) 
where the spatio-temporal beamformer w d now implements ISR without despreading 

n 

10 to reject / by complying with the following constraints (see Equation (15)): 



W d "Y d = 1, 
W d "C n = 0, 



W d Y d = 1, 

~"7°' n (42) 
W d I =0, 



~n ~n 



and C n is the constraint matrix without despreading that spans the interference subspace 
15 of the total interference vector / (/ e j ^ VeciC \) 

The constraint matrix without despreading, C n , is common to all low-power users. 
Thus, it characterizes the interference subspace regardless of the low-power user. In 
contrast, each constraint matrix after despreading C d PCMn in Equation (15) is obtained 
by despreading C n with the spreading sequence of the corresponding low-power user. 
20 Therefore ISR combining before despreading, although equivalent to beamforming after 
despreading, is computationally much more advantageous. 

In contrast to the "after despreading" case described earlier, when the data vector 
is not despread before processing by the ISR combiner (i.e., the constrained spatio- 
temporal beamformer) w d the estimate of the constraint matrix c is obtained by: 
25 " Q n -[c H n Cy, " (43) 

n„ = W.,) - C n Q n C" n , (44) 

n r 

W = " , (45) 

-n Kn r 

30 where I M * (2L _ 1} denotes a M * (2L - 1) x M * (2L - 1) identity matrix. As before, it can 
be seen from Equations (43) and (44) that the inverse matrix Q n is not the direct inverse 
of constraint matrix c but part of the pseudo-inverse of c - It should also be noted 
that the above operations are actually implemented in a much simpler way that exploits 
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redundant or straightforward computations in the data projection and the normalization. 
As before, the projector orthogonal to the constraint matrix c is formed once for 

n 

all low-power users. This would have not been possible with ISR after despreading. 
Second, the estimate of the low-power response vector y 4 is projected and normalized. 

— 0,n 

5 The estimate is reconstructed by reshaping the following matrix: 

the fast convolution with the channel matrix estimate being implemented row-wise with 
the spread sequence. The symbol <8> denotes overlap - add over the past, current and 
future blocks of the spread sequence to be convolved with a finite-size channel-matrix; 
10 hence is introduced in Equation 46 to isolate the net contribution from the current 

n 

sequence block. The channel vector estimate jj d 9 /.^ H d , is provided by STAR as 

n n 

explained earlier and includes the total contribution of the shaping pulse 4>(t) matched 
with itself [13]. If the channel time- variations are slow, the channel coefficients can be 
assumed constant over several symbol durations [20], thereby reducing the number of 
15 computationally expensive despreading operations required (see Figure 9). 

It should be noted that, although these ISR modes have formulations that are 
analogous whether ISR is implemented with or without first despreading the data vector, 
ISR combining of the data without it first being despread reduces complexity 
significantly. 

20 Receivers which implement these different ISR modes will now be described, 

using the same reference numerals for components which are identical or closely similar 
to those of the receiver of Figure 5, with a suffix indicating a difference. A generic ISR 
receiver which does so without despreading of the data will be described first, followed 
by one which does so after despreading of the data. Thereafter, specific implementations 

25 of different ISR modes will be described. 

Thus, Figure 9 illustrates a receiver according to a first embodiment of the 
invention which comprises a first set I of "strong user" receiver modules 20 1 ...20 NI 
which are similar to those in the receiver of Figure 5, and, separated by a broken line 
34, a second set D of "low-power" user receiver modules which differ from the receiver 

30 modules of set I but are identical to each other so, for convenience, only one, STAR 
module 20A d having a modified beamformer 47 A d , is shown. The outputs of the 
decision rule units 29\...,29 NI and of the channel identification units 28\...,28 NI from 
the set I modules are shown coupled to a constraints-set generator 42A which processes 
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the corresponding symbol estimates and channel vector estimates to produce a set of AT, 
constraints £ n = {Ci,...,C^}. The constraints-set generator 42A may, however, use 
hypothetical symbol values instead, or a combination of symbol estimates and 
hypothetical values, as will be described later. Each individual constraint lies in the 
5 same observation space as the observation matrix Y n from preprocessor 18. The 
constraints-set generator 42A supplies the set of constraints C n to a constraint matrix 
generator 43A which uses them to form a constraint matrix r and an inverse matrix 

n 

Q n which supplies it to the beamformer 47 d and each of the corresponding beamformers 
in the other receiver modules of set D. The actual content of the set of constraints C n 

10 and the constraint matrix c will depend upon the particular ISR mode being 
implemented, as will be described later. 

The observation vector deriving means in the receiver of Figure 9 also comprises 
a vector reshaper 44 which reshapes the observation matrix Y n from the preprocessing 
unit 18 to form an observation vector Y n having dimension M(2L-1) and supplies it to the 

15 beamformer 47 A d and to each of the other beamformers in the other receiver modules 
in set D. 

The STAR module 40 A d comprises a channel identification unit 28 A d , a decision 
rule unit 27 A d and a power estimation unit 30A d which are similar to those of the STAR 
modules 20 1 ...20 u described hereinbefore. The STAR module 20A d is associated with 
20 a despreader 19 d which also is part of the observation vector deriving means. The 
despreader 19 d despreads the observation matrix Y n using the spreading code for user d 
and supplies the resulting post-correlation observation vector z to the channel 

n 

identification unit 28 A d only. The decision rule unit 27 A d and power estimation unit 
30A d produce output symbol estimates b d n and power estimates (iffy 2 , respectively. 
25 The ISR beamformer 47 A d of STAR module 20A d produces corresponding signal 
component estimates § d n but differs from the MRC beamformers 27 i ...27 v/ because it 
operates upon the observation vector which has not been despread. In a manner 
similar to that described with respect to Figure 5, the channel identification unit 28 A d 
receives the post-correlation observation vector z d and the signal component 

n 

30 estimate $ d and uses them to derive the spread channel vector estimates which 
it uses to update the weighting coefficients w* of the beamformer 47 A d in succeeding 
symbol periods. The symbol period corresponds to the spread data frame of M(2L-1) 
elements. The coefficients of the ISR beamformer 47 A d also are updated in response to 
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the constraint matrix c and its inverse Q„, as will be described later. As shown in 

n 

Figure 9, the same matrices C n and Q n are supplied to all of the receiver modules in 
set D, specifically to their beamformers. 

As shown in Figure 10, the constraint matrix generator means 43A comprises a 
5 bank of vector reshapers 4&4 \ 4SA Nc and a matrix inverter 49A. Each of the 
vector reshapers 4SA \ 48A Nc reshapes the corresponding one of the set of 
constraints-set matrices c 1 , to form one column of the constraint 

matrix C n , which is processed by matrix inverter 49 A to form inverse matrix Q n . For 
simplicity of description, it is implicitly assumed that each of the columns of r is 

n 

10 normalized to unity when collecting it from the set of constraints C n . 

As also illustrated in Figure 10, beamformer 47 A d can be considered to comprise 
a coefficient tuning unit 50A d and a set of M(2L-1) multipliers 51?...5li (2£( _ 1) . The 
coefficient tuning unit 50A d uses the constraint matrix c 9 the inverse matrix Q n and the 
channel vector estimates to adjust weighting coefficients w d * W d * 

15 according to Equation (45) supra. The multipliers 51?... 5 1^.^ use the coefficients to 
weight the individual elements y Y respectively, of the observation 

—\,n'"—M{lL-\),n 

vector y . The weighted elements are summed by an adder 52 d to form the raw filtered 

— n 

symbol estimate § d n for output from the beamformer 47 A d . 

An alternative configuration of receiver in which the low-power STAR modules 

20 of set D implement ISR beamforming after despreading of the observation matrix Y n 
from preprocessor 18 will now be described with reference to Figures 11 and 12, which 
correspond to Figures 9 and 10. The receiver shown in Figure 11 is similar to that 
shown in Figure 9 in that it comprises a preprocessing unit 18 which supplies the 
observation matrix Y n to the set I receiver modules 20\...,20 NI , a constraints-set 

25 generator 42B and a constraint matrix generator means 43B. It does not, however, 
include the vector reshaper 44 of Figure 9 and each of the low-power user STAR 
modules in set D has a modified beamformer. Thus, modified beamformer 47B d 
operates upon the post-correlation observation vector z d from the output of the 

n 

despreader 19 d which is supplied to both the channel identification unit 28B d and the 
30 beamformer 47B d . The channel identification unit 28B d generates channel vector 
estimates H d and supplies them to the beamformer 47B d which updates its coefficients 
in dependence upon both them and a user-specific constraint matrix c d and user- 

PCM,n 

specific inverse matrix Q d PCMn . It should be noted that the constraint matrix generator 
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means 43B supplies user-specific constraint and inverse matrices to the other receiver 
modules in set D. 

Referring now to Figure 12, the common constraint matrix generator means 43B 
comprises a bank of user- specific constraint matrix generators, one for each of the 
5 receiver modules of set D, and each using a respective one of the spreading codes of the 
users of set D. Since the only difference between the user-specific constraint matrix 
generators is that they use different spreading codes, only user- specific constraint matrix 
43B d is shown in Figure 12, with the associated beamformer 47 A d . Thus, user-specific 
constraint matrix generator 43B d comprises a bank of 

10 despreaders 55B d - 1 , 55B d,Nc and a matrix inverter 46B d . The 
despreaders 55B d > l 9 1, 55B d ' Nc despread respective ones of the N c matrices in the set 
of constraints C n to form one column of the individual constraint 
matrix C d pcMn implicitly normalized to unity. The matrix inverter 46B d processes 
individual constraint matrix c d to form inverse matrix 0* The user-specific 

15 constraint matrix generator 43B d supplies the constraint matrix c d and inverse 
matrix Q d CMn to the coefficient tuning unit 50B d of beamformer 47B d . As shown in 
Figure 12, the beamformer 47B d has ML multipliers 51?...5l?„ which multiply 

1 ML, 

weighting coefficients w d \ . . W d * by elements z d . . . Z d of the post-correlation 

l,n ML,n 1,/T ' ' ML,n 

observation vector z d . As before, adder 52 d sums the weighted elements to form the 

— n 

20 signal component estimate § d . The beamformer coefficients are timed according to 
Equation (18). 

Either of these alternative approaches, i.e. with and without despreading of the 
data vector supplied to the beamformer, may be used with each of several different ways 
of implementing the ISR beamforming, i.e. ISR modes. It should be noted that all cases 

25 use a constraint matrix which tunes the ISR beamformer to unity response to the desired 
channel and null response to the interference sub-space. In each case, however, the 
actual composition of the constraint matrix will differ. 

Specific embodiments of the invention implementing the different ISR modes 
without despreading of the data will now be described with reference to Figure 13 to 20, 

30 following which embodiments implementing the same ISR modes after despreading will 
be described with reference to Figures 21 to 26; 
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Interference Subspace Rejection over Total Realisation (ISR-TR) 

The receiver unit shown in Figure 13 is similar to that shown in Figure 9 in that 

it comprises a set I of receiver modules 20*\ . . ,20 NI for processing signals of NI strongly 

interfering mobile stations and a set D of receiver modules for signals of other, "low- 

5 power", users. The receiver modules of set D are identical so only receiver module 

21C d , for channel d, is shown in Figure 13. As in the receiver of Figure 9, the 

observation matrix Y n from preprocessor 18 is supplied directly to each of the 

despreaders 19 1 ...19 NI of the set / receiver modules. Before application to each of the 

receiver modules of set D, however, it is delayed by one symbol period by a delay 

10 element 45 and reshaped by vector reshaper 44. The resulting observation vector y 

—n-i 

is supplied to the beamformer 47C d and to each of the other beamformers in the set D 
receiver modules (not shown). STAR receiver module 20C d is associated with 
despreader 19d and, in addition to beamformer 47C d , comprises channel identification 
unit 28C d , decision rule unit 29C d and power estimation unit 30C d which are similar to 
15 those shown in Figure 9. The set of channel estimates < K l n ^.^ e K N n \ which are supplied 
to the constraints-set generator 42C comprise the channel vector estimates H l ^...(f 1 and 

n n 

the power estimates ty 1 9 \J,™ 9 respectively. 

The constraints-set generator 42C comprises a bank of respreaders 57C 1 ...57C NI 
each having its output connected to the input of a respective one of a corresponding bank 

20 of channel replication units 59C 1 . . .59C NI by a corresponding one of a bank of multipliers 
SBC 1 ..^^ 1 . The respreaders 57C 1 ...57C NI are similar so only one, respreader 57C U , 
is illustrated in Figure 14. Respreader 57C U is similar to the corresponding spreader 13 u 
(Figure 3) in that it spreads the symbol £^ from the corresponding decision rule unit 
29C U using a periodic personal code sequence C/ " at a rate 1/T C , where T c is the chip 

25 pulse duration. It differs, however, in that it does not include a shaping-pulse filter. 
The effects of filtering both at transmission with the shaping-pulse (see Figures 2 and 3) 
and at reception with the matched shaping-pulse (see Figures 5 and 6) are included 
baseband in the channel vector estimate ff 1 or h u as disclosed in reference [13]. 

Referring again to Figure 13 and, as an example, receiver module 20C 1 , 

30 replication of the propagation characteristics of channel 14 1 is accomplished by digital 
filtering in the discrete time domain, i.e. by convolution at the chip rate of the channel 
vector estimate ff 1 with the respread data $ 1 C }. This filtering operation immediately 
provides decomposed estimates of the signal contribution of user station 10 1 to the 
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observation matrix Y n . Thus, respreader 57C 1 respreads the symbol $ l from decision 

n 

rule unit 29C 1 , multiplier 58C 1 scales it by the total amplitude estimate and channel 
replication filter 59C 1 filters the resulting respread symbol using the channel vector 
estimate fj l from channel identification unit 28C 1 . The symbol estimates from the other 

n 

5 STAR modules in set I are processed in a similar manner. 

It should be noted that the respreaders 57C 1 . . .57C NI , multipliers 58C 1 . . .58C NI and 
channel replication filters 59C 1 ...59C NI correspond to the elements 13 1 , 15 1 and 14 1 in 
the interfering user channel of Figure 2. The coefficients of the channel replication filter 
units 59C 1 . . .59C NI are updated in successive symbol periods by the channel identification 

10 units 28C 1 ...28C NI using the same coefficients H*..H™, corresponding to the 
transmission channels l4 l . m A4 M 9 respectively, used to update their respective MRC 
beamformers 27C 1 ...27C NI . It will be appreciated that the re-spread 
signals ... y^ 7 from the channel replication filter units 59C 1 ...59C NI , respectively, 

include information derived from both the sign and the amplitude of each symbol, and 

15 channel characteristics information, and so are the equivalents of the set I strong 
interferer's spread signals as received by the base station antenna elements 12 1 ...12 M . 

The constraint-set generator 42C also comprises an adder 60 coupled to the 
outputs of the channel replication units 59C 1 ...59C NI . The adder 60 sums the 
estimates y 1 y N \ of the individual contributions from the different interferers to form 

n-l n-l 

20 the estimate / a of the total interference from the NI interferers in the received 

n-l 

observation matrix Y n . The sum can be called total realization (TR) of the interference. 
In this embodiment, the constraint matrix generator simply comprises a single vector 
reshaper 43C which reshapes the total realization matrix / to form the 

n-l 

vector / which, in this embodiment, constitutes the constraint matrix C„. It should 

n-l 

25 be noted that, because the constraint matrix really is a vector, the inverse matrix Q n 
reduces to a scalar and, assuming implicit normalization, is equal to 1. Hence, no 
matrix inverter is needed. 

The reshaped vector / is supplied to the ISR beamformer 47C d of receiver 
module 20C d and to the beamformers of the other receiver modules in set D. The 

30 beamformer 47C d uses the reshaped vector / and the channel vector 

"n-l 

estimates y^ to update its coefficients, according to Equation (45), for weighting of 

0,n — 1 

the elements of observation vector y 

"~n-l* 
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The beamformer 47C d adjusts its coefficients so that, over a period of time, it will 
nullify the corresponding interference components in the observation vector y from 
the vector reshaper 44 and, at the same time, tune for a unity response to the spread 
channel vector estimate so as to extract the raw signal component 
5 estimate §*_ x substantially without distortion. 

ISR-TR constitutes the simplest way to characterize the interference subspace, yet 
the most difficult to achieve accurately; namely by a complete estimation of the 
instantaneous realization of the total interference vector / in a deterministic-like 
approach. The constraint matrix is therefore defined by a single null-constraint (i.e., 
10 N=l) as: 



I 

~~ n 
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~n 














/=i 



(47) 



15 where each estimate y* is reconstructed by reshaping the following matrix: 

yi = v^h; ® Bid. (48) 

For each interfering user assigned the index i — 1, NI, this mode uses 
estimates of its received power ) 2 and its channel h 1 , both assumed constant over 

n 

the adjacent symbols and made available by STAR. This mode also requires a bit-triplet 
20 estimate bl_ l9 b l n , b l n+l ] of each interfering user (see Equation (23)). To obtain estimates 
of the signs of the interferer bits for both the current and next 
iterations (i.e., b l n and the ISR-TR mode requires that the processing of all the 

low-power users be further delayed by one bit duration and one processing cycle (pc), 
respectively. The one-bit delay is provided by the delay 45 in Figure 13. 
25 In the ISR-TR mode and in the alternative ISR modes to be described hereafter, 

the interference (due to the strongest users) is first estimated, then eliminated. It should 
be noted that, although this scheme bears some similarity to prior interference 
cancellation methods which estimate then subtract the interference, the subtraction makes 
these prior techniques sensitive to estimation errors. ISR on the other hand rejects 
30 interference by beamforming which is robust to estimation errors over the power of the 
interferers. As one example, ISR-TR would still implement a perfect null-constraint if 
the power estimates were all biased by an identical multiplicative factor while 
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interference cancellers would subtract the wrong amount of interference. The next mode 
renders ISR even more robust to power estimation errors. 

The receiver illustrated in Figure 13 may be modified to reduce the information 
used to generate the interfering signal estimates y^ ... y* 7 specifically by omitting the 
5 amplitude of the user signal estimates, and adapting the ISR beamformer 47C d to provide 
more (NI) null constraints. Such a modified receiver will now be described with 
reference to Figure 15. 

Interference Subspace Rejection over Realisations (ISR-R) 

10 In the receiver of Figure 15, the receiver modules in set I are identical to those 

of Figure 13. Receiver module 20D d has the same set of components as that shown in 
Figure 13 but its beamformer 47D d differs because the constraint matrix differs. The 
constraints-set generator 42D differs from that shown in Figure 13 in that it omits the 
multipliers 58C 1 . . .58C NI and the adder 60. The outputs from the power estimation units 

15 SO^.^O™ are not used to scale the re-spread signals from the respreaders 57C 1 ...57C NI , 
respectively. Hence, in the receiver of Figure 15, the signals b 1 b NI from the STAR 

n n 

modules 20 1 ...20 NI , respectively, are re-spread and then filtered by channel replication 
filter units 59C 1 . . .59C NI , respectively, to produce user specific observation matrices 
y* j ... y^ p respectively, as the constraints-set C n . In contrast to the receiver of Figure 

20 13, however, these respread matrices are not summed but rather are processed 
individually by the constraint matrix generator 43D, which comprises a bank of vector 
reshapers 48D 1 ...4SD NI and a matrix inverter 49D (not shown but similar to those in 
Figure 10). The resulting constraint matrix £ comprising the column 
vectors y 1 ... f NI is supplied, together with the corresponding inverse matrix Q„, to 

25 each of the receiver modules in set D. Again, only receiver module 20D d is shown, and 
corresponds to that in the embodiment of Figure 13. Each of the 
vectors y 1 ... y^ 7 represents an estimate of the interference caused by the 

n-l «-l 

corresponding one of the strong interference signals from set I and has the same 

dimension as the reshaped observation vector y 

— /i-f 

30 In this ISR-R mode, the interference subspace is characterized by normalized 

estimates of the interference vectors y\ Consequently, it spans their individual 

— n 

realizations with all possible values of the total received powers (t//,) 2 . The constraint 
matrix is defined by NI null-constraints (i.e., N C =NI) as: 
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— n 





(49) 



where each estimate y* is reconstructed by reshaping the following matrix: 

f m = H n <8> bid. (50) 
5 It should be noted that, in the reconstruction of y* 9 the total amplitude of the 

i-th interferer (see Figure 15) has been omitted intentionally; hence the higher 
robustness expected to near-far situations as well as the enlarged margin for power 
control relaxation. 



10 Interference Subspace Rejection over Diversity (ISR-D) 

The ISR-D receiver shown in Figure 16 is predicated upon the fact that the signal 
from a particular user will be received by each antenna element via a plurality of sub- 
paths. Applying the concepts and terminology of so-called RAKE receivers, each sub- 
path is termed a "finger". In the embodiments of Figures 9, 11, 13 and 15, the channel 

15 identification units estimate the parameters for each finger as an intermediate step to 
estimating the parameters of the whole channel. In the ISR-D receiver shown in Figure 
16, the channel identification units 28E 1 ...28E NI supply the whole channel vector 
estimates fj 1 ... H NI respectively, to the beamformers 27 1 . . .27 NI , respectively, as before. 

n n 

In addition, they supply the sets of channel estimates <yC n ... comprising a sub- 
20 channel vector estimate for each individual sub-channel or finger, to the constraints-set 
generator 42E. The set of channel estimates ^ comprises the sub-channel vector 
estimates ... ,H„ f . The constraints-set generator 42E is similar to that shown in 
Figure 15 in that it comprises a bank of respreaders 57 ! ...57 NI but differs in that the 
channel replication units 59D ! ...59D NI are replaced by sub-channel replication units 
25 59E 1 . . . 59E NI , respectively . The sub-channel replication units 59E 1 . . .59E NI convolve the 
respread symbols with the sub-channel vector 
estimates jy^ 1 ... H^; ... ;H™ ,1 > ... ,H™ Mf respectively, to produce normalized 
estimates y lj \ ... y 1 ^- -y* 7 ' 1 of the sub-channel-specific observation matrices 

n~\ n-\-> ' rc-1 ' 5 n-l ^ 

decomposed over fingers. Hence, the matrices span the space of their realizations with 
30 all possible values of the total received powers (^) 2 and complex channel 
coefficients ff The estimates are supplied to a constraint matrix generator 43E which 
generally is as shown in Figure 10 and produces the constraint matrix accordingly. 
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The constraint matrix c is simply defined by A^W null-constraints (i.e., iV c 
N f x NI = M x P x NI) as: 



C n = 



— n 




— n 




— n 




— n 


n 




n 




— n 




— n 



(51) 



Each estimate f J is reconstructed by reshaping the following matrix: 

— n 

Y'/ = & s n (8) b 1 c{. 

n n n 1 



(52) 



It should be noted that, in the reconstruction of y^, the total amplitude of the 

10 /-th interferer as well as the channel coefficients g (see Figure 1) are 
intentionally omitted; hence the relative robustness of ISR-D to power mismatch, like 
ISR-R. Unlike other modes, it additionally gains robustness to channel identification 
errors and remains sensitive only to the estimated channel parameters remaining, namely 
the multipath time-delays, and to symbol estimation errors. 

15 It should be noted that, in the receivers of Figures 13, 15 and 16, estimation 

errors of the interference bit signs may introduce differences between the estimated 
constraints and the theoretical ones. Hence, although ISR-D, ISR-R and ISR-TR modes 
are satisfactory in most situations, it is possible that the realisation could' be erroneous, 
which would affect the validity of the interference cancellation. Additionally, estimation 

20 of the signs of the interference bits for reconstruction in the ISR-D mode, as in the ISR- 
R and ISR-TR modes, requires that the processing of all of the low-power users be 
further delayed by one bit duration, i.e., by delay 45, and one processing cycle (pc). 
To avoid these drawbacks, alternative ISR approaches to implementation of the 
constraints of Equation (42) are envisaged and will now be described, beginning with 

25 ISR-H which avoids processing delays and is completely robust to data estimation errors. 



Interference Subspace Rejection over Hypotheses (ISR-H) 

It is possible to use a set of signals which represent all possible or hypothetical 
30 values for the data of the interfering signal. Each of the interfering signals constitutes 
a vector in a particular domain. It is possible to predict all possible occurrences for the 
vectors and process all of them in the ISR beamformer and, therefore, virtually guarantee 
that the real or actual vector will have been nullified. As mentioned, the strong 
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interferers are relatively few, so it is possible, in a practical system, to determine all of 
the likely positions of the interference vector and compensate or nullify all of them. 
Such an alternative embodiment, termed Interference Subspace Rejection over 
Hypotheses (ISR-H) because it uses all possibilities for the realisations, is illustrated in 
5 Figure 17. 

The components of the "interferer" receiver modules of set I, namely the 
despreaders ^...lS^ 1 and STAR modules 20 1 ...20 NI , are basically the same as those in 
the receiver of Figure 15 and so have the same reference numbers. In the embodiment 
of Figure 17, however, the constraints-set generator 42F differs because the symbol 
10 estimates b\ ... 5™ from the outputs of the decision rule units 29 1 . . .29 NI are not supplied 
to the respreaders 57F 1 ...57F NI , respectively, but are merely outputted to other circuitry 
in the receiver (not shown). 

Instead, bit sequence generators eSF^.^F* 1 each generate the three 
possibilities g 1 , g 2 , g 3 , which cover all possible estimated values of the previous, current 

n n n 

15 and next bits of the estimated data sequences b\...b™, including the realisation itself (as 
explained later), and supply them to the respreaders 57F 1 . ..57F NI , respectively, which 
each spread each set of three values again by the corresponding one of the spreading 
codes. The resulting re-spread estimates are filtered by the channel replication filters 
59F 1 ...59F NI , respectively, to produce, as the constraint set, the matrix estimates 

20 y\ , F, , F 1 , ; ... ; Y? 1 , Y N * , Y N ! . The bit sequence generators could, of course, be 
replaced by storage units. 

The constraint matrix generator 43F is generally as shown in Figure 10 and 
processes the set of estimate matrices to form the column vectors 
Y x , F , F 1 lY^F* 7 , Y NI of constraint matrix Q which it supplies with 

25 corresponding inverse matrix Q n , in common to the beamformer 47F d and the 
beamformers of the other set D receiver modules. 

Receiver module 20F d comprises similar components to those of the receiver 
module 20E d shown in Figure 16. It should be noted, however, that, because the "next" 
bit is being hypothesized, it need not be known, so the delay 45 is omitted. . 

30 As mentioned above, the two bits adjacent to the processed bit of the i-th 

interferer contribute in each bit frame to the corresponding interference vector (symbol) 
to be rejected. As shown in Figure 18, enumeration of all possible sequences of the 
processed and adjacent bits gives 1? = 8 triplets, each of three bits. Only one of these 
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triplets could occur at any one time at each bit iteration as one possible realization that 
generates the user-specific observation matrix y* These eight triplets can be identified 
within a sign ambiguity with one of the four triplets identified as (a). . . (d) in the left-hand 
part of Figure 18, since the four triplets (e)...(h) are their opposites. 
5 It should be appreciated that the bit sequence generators 63 ^.^S™ (Figure 17) 

each supply only three values, o 1 , g because the dimension of the generated signal 

n n n 

subspace is 3. It should be noted that frames of duration 3T, taken from these sequences 
at any bit rate instant, reproduce the eight possible realisations of the bit triplets of 
Figure 18. Therefore, at any bit iteration, the bit sequence b^ of the interfering mobile 

10 station can be locally identified as the summation of the generating sequences 
k = 1,..., 3 weighted by the bit signs b*_ l9 b* and b* +1 . Replacing the estimate 
B l in Equation (50) by R k fc = i 3 , yields canonic observation matrices that 
span all possible realisations of the received signal vector from the /-th interfering mobile 
within a sign ambiguity. 

15 In the ISR-H embodiment of Figure 17, the interference subspace is characterized 

by normalized estimates of the canonic interference vectors f 1 . Accordingly, it spans 
their individual realizations with all possible values of the total received 
powers (^)2 and bit triplets [b^ u b^b^ +l . The constraint matrix is defined by 3NI 
null-constraints (i.e., N c = 3NI) as: 



20 



c n = 



y 1 y 1 

— 0,ai ~-\,n 



yM yNI 

"^0,n - 1 ,n — +l,n 



—0,n 



--\,n 



— +l,n 



—0,n 



— + 1,/! 



(53) 



25 



where each estimate is reconstructed, respectively, for k = -1, 0, +1 by reshaping 
the following matrix: 

Kn =K ® St?cl. (54) 
It should also be noted that, in the reconstruction above, only the channel vector 
estimates (assumed stationary over the adjacent symbols) are needed for complete 
interference rejection regardless of any 2D modulation employed (see Figure 19); hence 
the extreme robustness expected to power control and bit/ symbol errors of interferers. 
30 The ISR-H combiner coefficients are symbol-independent and can be computed less 
frequently when the channel time- variations are slow. 

Merging of the D mode with the H mode along the decomposition of Equation 
(38) yields ISR-HD (hypothesized diversities) with a very close form to the decorrelator. 
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This ISR-HD mode requires a relatively huge number of constraints (i.e., 3N f NI). 
Consequently, the ISR-HD mode is not considered to be practical at this time. 

In fact, it would be desirable to reduce the number of constraints required by the 
ISR-H receiver described above. This can be done using an intermediate mode which 
5 is illustrated in Figure 20 and in which the receiver modules of both sets I and D are 
similar to those of Figure 15; most of their components are identical and have the same 
reference numbers. In essence, the constraint-set generator 42G of the receiver in Figure 
20 combines the constraint-set generators of Figures 15 and 17 in that it uses estimated 
symbols and hypothetical values. Thus, it comprises a bank of respreaders 57G 1 . . .57G NI , 

10 a corresponding bank of channel replication units 59G 1 ...59G NI and a bank of symbol 
generators 63G 1 ...63G NI . In this case, however, each of the symbol generators 
63G 1 ...63G NI supplies only one symbol to the corresponding one of the respreaders 
57G 1 ...57G NI , which receive actual symbol estimates ... $™ 9 respectively, from the 
decision rule units 29 1 ,. ..,29 NI , respectively. It should be appreciated that, although the 

15 symbol generators 63G 1 ,. . .,63G NI each supply only one symbol for every actual symbol 
or realization from the corresponding one of the decision rule units 29 1 ,...,29 NI , that is 
sufficient to generate two hypothetical values of "future" symbols b^ l9 ... 9 b^ x for 
every one of the symbol estimates b* +l ... b™ x since only two hypothetical values of the 
symbols, namely 1 and -1, are required. The respreaders 57G 1 ,...,57G NI supply the 

20 spread triplets to the channel replication units 59G 1 . . .59G NI which filter them, using the 
channel vector estimates ft 1 ... fj NI respectively, to produce pairs of 

n n 9 

matrices y 1 , f\ : ... Y N \ Y N ! and supply them to the constraint matrix generator 43G 
which is configured generally as shown in Figure 9. The constraint matrix generator 
43G reshapes the matrices y 1 y 1 • y NI v NI to form 

25 vectors y 1 y 1 • ... y N \ Y NI which then are used as the column vectors of the constraint 
matrix c . The constraint matrix generator 43G supplies the constraint matrix c and 
the corresponding inverse matrix Q n in common to the beamformer 47G d and the 
beamformers of other receiver modules in set D. 

Hence, the beamformer 47G d uses the past symbol estimate $ i of the 

30 interference data as well as the present one b l n (delayed by one processing cycle, Le. 
the time taken to derive the interference estimates), and the unknown sign 
of reduces the number of possible bit triplets and the corresponding realisations for 
each interference vector to 2. 
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The receiver of Figure 20, using what is conveniently referred to as ISR-RH 
mode for reduced hypotheses over the next interference bits, rejects reduced possibilities 
of the interference vector realisations. Compared to the receiver of Figure 17 which uses 
the ISR-H mode, it is more sensitive to data estimation errors over h l and h 1 and 

n-\ n 

5 requires only 2 constraints per interferer instead of 3. 

Using the previous and current bit estimates of interferers, uncertainty over the 
interference subspace can be reduced and it can be characterized by the following matrix 
of 2NI null-constraints (i.e., N c = 2NI): 



10 



C n = 



where: 



~r,n 



y 1 



— r,n — +1,/? 



— r,n — 



(55) 



(56) 



15 and where each estimate f 4 is reconstructed by reshaping the matrices in Equation 
(38), respectively for k = -1,0, + 1. It should be noted that this mode requires a delay 
of one processing cycle for the estimation of the current interference bits. 

The ISR-RH mode has the advantage of reducing the number of null-constraints 
as compared to the ISR-H mode. A larger number of null-constraints indeed increases 

20 complexity, particularly when performing the matrix inversion in Equation (43), and may 
also result in severe noise enhancement, especially when the processing gain L is low. 
As the number of strong interferes NI increases in a heavily loaded system, the number 
of null-constraints (2A7and 3NI) approaches the observation dimension M x (2L - 1) and 
the constraint-matrix may become degenerate. To reduce complexity, guarantee stability 

25 in the matrix inversion of Equation (43), and minimize noise enhancement, the constraint 
matrix C n in Equations. (43) and (44) is replaced by the orthonormal interference 
subspace of rank K that spans its column vectors as follows: 

y. = v*c{£) = (57) 

to yield the projector U n = I M<2L .i } ~ V n Vn used in Ec l uation ( 31 )> which replaces 
30 matrix inversion operation by matrix orthonormalization. The Gram-Schmidt 
orthonormalization procedure is used, which implements a cumulative increasing-rank 
cascade of "corank" -1 projections, like in linear SIC [28], except that projections there 
are formed differently for direct cancellation. Ultimately, after orthonormalization, ISR 
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exploits a "corank"-7V c projection to implement N c null-constraints in the combining step. 
This projection could be implemented as a cascade of N c corank-1 projections. It should 
be noted that orthonormalization becomes unnecessary if we check that c is close to 

n 

orthonormal (i.e., c ~ V )• 

n n 

5 In practice, can hardly reflect the real rank of c It corresponds to the 

subspace of reduced rank k < N with the highest interference energy to cancel. To 
further minimize noise enhancement, one can also increase the observation dimension M 
x (2L -1), as will be described later as "X option", and so on. 

Another alternative that reduces noise enhancement constrains the beamformer to 

10 implement a close-to-null response towards each null-constraint instead of an exact null- 
response. This "relaxation" of the null-response leaves more degrees of freedom for 
ambient noise reduction (i.e. , less amplification). It usually amounts to upper-bounding 
the amplitude of the beamformer response towards each null-constraint to be less than 
a maximum threshold. This technique is well-known and classified in the literature as 

15 a "robust beamforming" method. We can combine it with ISR to reduce noise 
enhancement. 

Constraint relaxation in robust beamforming is usually solved as an optimization 
problem under constraints using the Lagrange multipliers technique. Without going into 
the mathematical details of such derivations, we directly provide intuition-based solutions 
20 that extend ISR beamforming in a seamless manner. We extend Equations (43) to (45) 
as follows: 

C„ = C n A n , (ID 

Q n = [c"c n + u N )-\ (12) 

25 

n„ = w_ 0 - C n Q n C" n , (13) 

n y" 

W" = —0Z2f—, (14) 

30 where A n is a ]V c x 7V c diagonal matrix of positive weights and X an additional 
weighting factor. 

The above extended ISR solution reduces to that of Equations (43) to (45) by 
setting A n = I N and \ = 0. It also covers the particular case of MMSE (minimum 
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mean square error) combining also known in the beamforming literature as MVDR 
(minimum variance distortionless response) beamforming. Along the model equation 
(28), the MMSE or MVDR beamforming solution is given by: 

5 W = - , (15) 

~" Y dH Ril N Y d 

where R } n is the correlation matrix of the interference-plus-noise vector I + N . 

- — ~ n n 

Assuming that the additive noise vector is spatially and temporally uncorrelated with 
variance a 2 N , R UN is given by: 

10 = c n V" + (16) 

where C n is the constraint matrix defined along one of the ISR modes described 

previously, but without normalization column-wise. Depending on the signal parameters 

assumed unknown a priori (c/ , discussion above Equation (15)), the conditional 

statistical mean in r gives rise to a particular interference decomposition 

15 in C n along the corresponding ISR mode. The N c X N c diagonal matrix $ holds 

the power of the unknown signal parameters along which interference is decomposed. 

In the TR mode, interference characterization is quasi-deterministic with: 

= til- 

In the R mode, we have: 
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*„ = diag [(V n ) 2 ,...,(t™) 2 



while in the D mode, we have: 



25 *„ = diag [(fcfifi ^ i 2 1 P^. i 2 (?ir) 2 1 is. i 2 1 ^ i 

where i^ 1 ^ and n \ 2 both denote local averages in time of (i//) 2 and |^„| 2 , 
respectively. The H mode gives: 

30 *. - <uag [t?l fcf. (Vl ... («f, (Wf . (K') 2 ]. 



Using the inversion lemma, the inverse of R can be written as: 

I+N 
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J? 1 



2 

<J N 



l M*(2L-\) 



~ C rPn [®n C n C^ n + 0 N I „ } C„ 



(17) 



The MMSE or MVDR beamforming can therefore be identified with the extended ISR 
solution using: 



A = a 



(18) 



and: 



[1] in the TR mode, 

$f in the R mode, 

7 1 7 1 7 1 n NI n NI n NI 



in the D mode, 



in the H mode. 



10 The instantaneous estimates in the equation above can be further averaged or smoothed 
in time. The columns of the constraint matrix estimate c are reconstructed without 

n 

normalization. 

Although the MMSE version of ISR serves as a method to reduce noise 
enhancement, it is still sensitive to hard-decision errors for modes using decision feed- 
15 back. Using weights in a different way is a method to reduce sensitivity to hard-decision 
errors. We notice that the ISR combined signal estimate can be formulated as: 



w dH Y = (Y d "n r \ l r" (Y - CO C H Y V 

— n —n \— 0,n 0,n) —0,n \—n nV£ -» "~n /' 
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\—0,n n —0,nj —0,n \~n ~n 



-I 



(19) 



(HO) 



(111) 



The last reformulation stresses that c v can be understood as a sum of estimated 
25 interference vectors corresponding to the ISR mode applied. For instance each of these 
estimates corresponds to an interfering user in ISR-R. It can be shown that if tentative 
decisions, used to form the /-th constraint, are wrong, the corresponding constraint will 
sometimes be amplified rather than rejected 2 . This amplification is mitigated by 
considering the alternative projector: 



30 2 This is true when the tentative decision, which is in error, is temporally 
overlapping the current bit. 
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H A .„ = / - C n Q n \ n C" n , (112) 
where A n is a diagonal matrix of weights. Using this projection, Equation (111) is 
modified as: 

w""y = (£f n An y" TV" (y - \a\ V - - - {a } i N \ (ii3) 

5 

This modification means that the interferer is never completely rejected although tentative 
decisions are all correct; however, the penalty arising from wrong decisions is reduced 
as well. For ISR-R and ISR-D it can be shown that the optimal weight to be applied to 
columns representing interferer / is (A ) (1 - 2SER(i)) where SER(Z) is the Symbol 

10 Error Rate of the tentative decisions associated with interferer /. 

It should be noted that each of the receivers of Figures 13, 15, 16, 17 and 20 
could be modified to perform ISR "after despreading" of the observation matrix Y n , in 
effect in much the same way that the generic "after despreading" receiver of Figure 11 
differs from the generic "without despreading" receiver of Figure 9. Such modified 

15 receivers will now be described with reference to Figures 21 to 26. 

Thus, in the ISR-TR receiver shown in Figure 21, which corresponds to that 
shown in Figure 13, the delay 45 delays the observation matrix Y n from the 
preprocessing unit 18 by 1 bit period and supplies the resulting delayed observation 
matrix Y n _ }> in common, to each of the low-power user receiver modules in set D. Only 

20 one of these receiver modules, 20H d , is shown in Figure 21, since all are identical. The 
observation matrix Y n _j is despread by despreader 19 d and the resulting post-correlation 
observation vector z d is supplied to both the channel identification unit 28H d and the 
beamformer 47H d . The receiver modules of set I and the constraints-set generator 42C 
are identical to those in the receiver shown in Figure 13, and supply the matrices 

25 f l n _ x ... to an adder 60 which adds them to form the total interference 

matrix / which it supplies to each of the receiver modules in set D. 

Receiver module 20H d is similar to that shown in Figure 13 but has a second 
despreader 43H d which uses the spreading code for user d to despread the total 
interference matrix to form the user-specific constraint matrix as a single column 

30 vector [ d m This despreader 43H d , in effect, constitutes a user-specific constraint 

~PCM,n-\ ^ 

matrix generator because the constraint matrix is a vector and an inverse matrix is not 
needed. Also, in this case, the channel identification unit 28H d supplies the channel 
vector estimate ff d to the beamformer 47H d . 

n-\ 
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It should be noted that the despread data vector z d is equal 

— n-l 

to H d sf + I d , + N* where ff d is the channel response for user station 10 d , c rf is 

n n ~PCM,n ~n — n n 

the signal transmitted by the mobile station 10 d of user d, and j d is the interference 

~PCM y n 

component present in the signal z d as a result of interference from the signals from the 

n 

5 other user stations 10 1 in set I, where is as defined in Equation (14). The 

-PCM,n 

value N d is additional noise which might comprise, for example, the summation of 

PCM t n 

the interference from all of the other users on the system at that time, as well as thermal 
noise. "Other users" means other than those covered by the channels in set I. 

As before, the coefficients of the beamformer 47H d are tuned according to 
10 Equations (16) to (18) and the constraint matrix is defined by a single null-constraint 
(i.e., N c = l) as: 



'PCM,n 



15 



~PCM,n 




NI 

-PCM,n 

i=\ 


i d 






\\NI 1 


-PCM,n 












pr ~ pcM - n 



(58) 



where the estimate f d is obtained by despreading the matrix J n (see Equations (47) 
and (48)) with the spreading sequence of the desired low-power user. 

Figure 22 shows a similar modification to the low-power (set D) receiver modules 
of the "without despreading" ISR-R receiver of Figure 15. In this case, the output of 
20 the constraint-set generator 42D, as before, comprises the matrices y 1 As 

n-l "' n-V 

before, only receiver module 20 J d is shown in Figure 22 and is identical to that shown 
in Figure 21 except that the second despreader 43H d is replaced by a user-specific 
constraint matrix generator 43 J d of the kind shown in Figure 12. The channel 
identification unit 28 J d again supplies the vector u d to the beamformer 47 J d . The 
25 bank of despreaders in the user-specific constraint matrix generator 43J d despread the 
respective ones of the matrices y 1 y NI to form the vectors 

n-l '*' n-l 



-pCM,n-\ 

matrix 



f M which constitute the columns of the user-specific constraint 

~PCM,n-\ 

and the matrix inverter 46B d produces the corresponding inverse 



C d 

QpcMn-v Bot ^ °^ these matrices are supplied to the associated beamformer 47J* 



matrix 

30 which uses them and the channel vector estimate H d to adjust its coefficients that are 

n-l 

used to weight the elements of the post-correlation observation vector z d As before, 

— n-l* 

the coefficients are adjusted according to Equations (16) to (18) and the constraint matrix 
is defined by NI null-constraints (i.e., N C —NI) as: 



'PCM.n 
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~PCM t n 




-PCM,n 








jdM 
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~PCM,n 




-PCM,n 







(59) 



where each estimate is obtained by despreading the matrix y 1 ' of Equation (50) 

PCM,n n 

with the spreading sequence of the desired low-power user. 
5 Figure 23 illustrates the modification applied to the low-power user receiver 

module of the ISR-D receiver of Figure 16. Hence, there is no common matrix inverter. 
Instead, in the receiver of Figure 23, each of the receiver modules of set D has a user- 
specific constraint matrix generator 43K which receives the constraints from the 
constraints-set generator 42E. As illustrated, user-specific constraint matrix generator 
10 43K d processes the sets of matrices y u y 1 ^ - f NI > 1 to form the set of 

vectors ?> l > 1 ... f XN f • ... jW.i jww f which constitute the columns of user- 

-pCM,n-l -pCM,n~l J PCM,n-l 9 ~PCM,n-\ 9 

specific constraint matrix , and the corresponding inverse 

matrix Qp CMn _ x which it supplies to the beamformer 47K d . As before, the beamformer 
47K d tunes its coefficients according to Equations (16) and (18). The constraint matrix 
15 is defined by NjNI null-constraints (i.e., N c = N f x NI = M x P x NI) as: 



'PCM,n 



jfd, 1,1 Jd,\,N f 

~PCM,n ~PCM,n 
, ... , 



fd,MA 
~PCM,n 



Jd,NI,N f 
~PCM,n 



jd,\ t \ 
~PCM,n 



jd.\,N f 
-PCM.n 



£d,NIA 
-PCMm 



-PCM,n 



(60) 



where each estimate is obtained by despreading f J of Equation (52) with the 

PCM,n — n 

20 spreading sequence of the desired low-power user. 

Figure 24 illustrates application of the modification to the ISR-H receiver of 
Figure 17. Again, the common constraint matrix generator (43F) is replaced by a user- 
specific constraint matrix generator 43L d in receiver module 20L d and similarly in the 
other receiver modules of set D. The constraints-set generator 42L differs from 

25 constraints-set generator 42F of Figure 17 because its bit sequence generators 
63L 1 ,...,63L NI use different generating sequences. The sets of matrices 
f\ , f\ , f\ ; yf 7 , yf 7 , f™ from the constraints-set generator 42F are processed 
by the user-specific constraint matrix generator 43 L d to form the 
vectors f' 1 , ?> x j* 1 - ... -f^\ f** 1 f^ 1 which constitute the columns of the user-specific 

-\,n 9 ~2,n ~Xn' 9 ~2,n 9 S s n 

30 constraint matrix cjL w , and the matrix inverter (not shown) produces the 

PCMji 9 

corresponding inverse matrix oi™ The constraint matrix r d and the inverse 

^PCM,n m PCM,n 

matrix Q d PCMn are used by the beamformer 47L d , together with the channel vector 
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estimate fif to adjust its coefficients that are used to weight the elements of the DOSt- 
rt r 

correlation observation vector z d received from despreader 19 d . As before, the 

n 

coefficients are adjusted according to Equations (16) and (18) and the constraint matrix 
is defined by 3NI null-constraints (i.e., N c = 3NI) as: 



'PCM,n 



jdAA 
-PCM,n 



pi f l, 2 ^rf,l,3 
~PCM,n ~PCM,n 



-PCM,n 



-PCM,n 



-PCM,n 



jdAA 
-PCM,n 



pi A, 2 
~PCM,n 



JdA3 
-PCM,n 



pi MA 
-PCM,n 



pi,NI,2 
~PCM,n 



-PCM,n 



(61) 



where each estimate f^ k is obtained by despreading the matrix f* with the 

~PCM,n k,n 

spreading sequence of the desired low-power user. 
10 In this case, each of the bit sequence generators 63L 1 ...,63L NI uses four 

generating bit sequences g l (t), g 2 (t), g\t) and g\t) as shown in Figure 25. 

It should be noted that, in any frame of duration 3T in Figure 25, a bit triplet of 

any of the four generating sequences is a linear combination of the others. Therefore, 

any one of the four possible realisations of each interference vector is a linear 
15 combination of the others and the corresponding null-constraint is implicitly implemented 

by the three remaining null-constraints. The four null-constraints are restricted 

arbitrarily to the first three possible realisations. 

Figure 26 illustrates application of the modification to the ISR-RH receiver of 

Figure 20. Again, the common constraint matrix generator 43G of Figure 20 is replaced 
20 by a set of user-specific constraint matrix generators, 43M d in receiver module 20M d and 

similarly in the other receiver modules of set D. The constraints-set generator 42M 

shown in Figure 26 differs slightly from that (42G) shown in Figure 20 because each of 

the bit sequence generators eSM^^^M* 41 in the receiver of Figure 26 generate the bit 



sequence g 



25 



The user-specific constraint generator 43M d processes the pairs of constraint-set 



matrices f) , f) ; • Y™ , Y? 1 from the channel identification units 59M 1 ,...,59M NI , 



Wl 



respectively, by 



to produce the corresponding set of vectors 



jdAJc x jdA^ jdjujc, piMJ^ W hi c h constitute the columns of the user-specific constraint 

-PCMji -PCM,n 9 ' ~PCM y ri ~PCM,n 



matrix C d pCMn , and to produce the corresponding inverse matrix Q d 



The 



30 constraint matrix and the inverse matrix n 



*PCM,n 



PCMji' 

are used by the beamformer 



47M d , together with the channel vector estimate H d 9 to adjust its coefficients that are 

n 

used to weight the elements of the post-correlation observation vector z d received from 
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despreader 19 d . As before, the coefficients are adjusted according to Equations (16) to 
(18) and the constraint matrix is defined by 2NI null-constraints (i.e., N c = 2NI) as 
follows: 



10 



^PCMm 



-PCMm 



-PCMm 



-PCMm -PCMm 



-PCMm 



fd,Uk 2 
-PCMm 



jd,M,k t 
-PCMm 



jd.NI^ 
-PCMm 



(62) 



where each pair of estimates i dJ - k > and f dJ •*« is obtained by despreading the matrices 

-PCMm -PCMm 

Y* k n and respectively, with the spreading sequence of the desired low-power 

user. 



Inter-Svmbol Interference (ISI) Rejection 

In any of the above-described embodiments of the invention it may be desirable 
to reduce inter-symbol interference in the receiver modules in set D, especially when low 
processing gains are involved. As noted in the PCM model where despreading reduces 
15 ISI to a negligible amount, for a large processing gain, 
Y^Y_ d Xn = 0 and Y^ H Y d « 0. Hence, the before despreading spatio-temporal 
beamformer w d approximately implements the following additional constraints: 



20 



W^Y* ^ 0, 

n -\m 
n — +1m 



(63) 



Accordingly, it rejects interference and significantly reduces ISL Complete ISI 
rejection can be effected by modifying the receiver to make the set of the channel 
parameter estimators available to the constraints-sets generator 42 for processing 

25 in parallel with those of the set I receiver modules. The resulting additional constraint 
matrix and inverse matrix would also be supplied to the beamformer 47 d and taken into 
account when processing the data. 

In such a case, the following matrix can be formed: 
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and the following 2x2 matrix 



inverted to obtain the constrained spatio-temporal beamformer y/ d before despreading 

n 

5 by: 

n /l/,„ = I Mt( 2L-l) ~ ^Sl,nQlSI.n^ sln , (66) 



u d n = n? SI , n n„, (67) 

10 w" = !Z2^_. (68) 

—0,n n —0,n 

The projector n n is produced in the manner described earlier according to 
Equations (43) and (44). The projector n d orthogonal to both c d and c is formed 

n ISI,n n ' 

and then the low-power response vector y 4 is projected and normalized to form the 
15 beamformer which fully rejects ISI from the processed user d and interference from the 
NI users in set I. 

It should be noted that, if the suppression of strong interferers is not needed, ISI 
can still be rejected by the following beamformer: 

20 W « = n ' s '-^, f (69) 

~ n r H ni SIn r 

— 0,/i 0,« 

where projector n n in Equation (67) would be set to identity and hence would have no 
effect. This is the same as setting the matrix c to null matrix. If the 

n 

projector uf SIn in the above equation is replaced by an identity matrix, (equivalent to 
25 setting matrix to null matrix) then a simple MRC beamformer is implemented 

before despreading. A receiver module using such an MRC beamformer is illustrated 
in Figure 27 and could be used to replace any of the "contributor" only receiver 
modules, such as receiver modules 20 1 ,. . . ,20 NI in Figure 9 et seq. The receiver module 
shown in Figure 27 is similar to receiver module 20A d of Figure 9 except that the ISR 
30 beamformer 47 A d is replaced by an MRC beamformer 27N d which implements the 
equation 

T 

W d = °' K (70) 

o,n 
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It is also envisaged that the receiver module of Figure 27 using the MRC 
beamformer denoted y/ d in the following, could be incorporated into a STAR which 

MRC,n 

did not use ISR, for example the STAR described in reference [13]. 

5 Pilot-Assisted ISR 

STAR-ISR performs blind channel identification within a sign ambiguity (or 
quantized-phase ambiguity for MPSK). It hence avoids differential demodulation and 
enables quasi-coherent detection. However, differential decoding is still required at the 
cost of some loss in performance to resolve the sign ambiguity. To implement full 

10 coherent detection and avoid differential decoding, a pilot signal, i.e., a predetermined 
sequence of symbols known to the receiver (usually a constant " 1 " sequence), can be sent 
by the transmitter to enable the receiver to resolve the sign ambiguity. Two pilot types 
are common namely 1) A pilot-symbol will insert pilot symbols in the data sequence at 
predetermined symbol positions or iterations n the other indices n being allocated 

15 to the data symbols with some overhead loss; 2) A pilot-channel will multiplex the pilot 
sequence with the data sequence using a pair of distinct spreading 
codes c ?' u and c /'"; a fraction + f) of the power available being allocated to 

the pilot, the rest i/(i + to the data with some relative power loss. 

Pilot signals are usually used to perform channel identification. STAR-ISR 

20 achieves this task without a pilot and hence reduces the role of the pilot to a simple 
resolution of the phase ambiguity resulting from its blind channel identification approach. 
This new approach to pilot use enables significant reduction of the overhead or power 
fraction allocated to the pilot, as disclosed in references [21] and [31] in the single-user 
context of STAR on both uplink and downlink, respectively. Insertion of a pilot in 

25 STAR-ISR along the new approach is depicted by Figures 44 and 46 for a pilot-symbol 
and a pilot-channel, respectively. 

In Figure 44, a pilot-symbol assisted ISR receiver is shown for user d. The ISR 
beamformer 47V d outputs the signal component estimates § d n and interacts with the 
channel identification and power control units 28V d and 30V d in the regular steps 

30 described earlier. However, an additional pilot/data demultiplexer 35V d at the 
beamformer output isolates the symbol positions n and n allocated to the data and 
the pilot. It hence extracts at the corresponding index positions two distinct streams of 
signal component estimates s d n and 9 respectively. 
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On one hand the data signal component estimates s d ^ a d \jf n b d delivered within 
a quantized phase ambiguity a d (i.e., belongs to the constellation) are fed to the 
decision rule unit 29V d to estimate the corresponding symbol 
estimates t> d nb — a d b d within the same phase ambiguity a d . 
5 On the other hand the pilot signal component estimates § d — a d $ d n , which 

basically hold the constant phase ambiguity within little fluctuation of the received 
power |^j 2 9 are fed to an ambiguity estimator 3 1 V d shown in Figure 45 . This estimator 
buffers the pilot signal component estimates over a predetermined number of symbols. 
Once the buffer 33 V d is full, the estimator averages (or smooths) the buffered values to 

10 significantly reduce the residual noise and hence minimize estimation errors using the 
averaging or smoothing unit 34 V d . It then feeds the resulting average pilot signal 
component estimate to a decision rule unit 29V/2 d identical to the one used for the data 
(i.e., 29V d ) to quantize and estimate the phase ambiguity a d . Once a new estimation 
is made, the buffer 33V d is flushed and the phase ambiguity estimate is kept constant 

15 until a new estimate is derived when the buffer is filled again. 

In the pilot-symbol assisted ISR receiver of Figure 44, a conjugator 32 V d 
conjugates the ambiguity estimate and a multiplier 15V d multiplies the resulting phase 
conjugate with the data symbol estimate $ d ~ a d b„ for phase ambiguity 

compensation. The resulting symbol estimate b d — b is ambiguity-free and no longer 

— n n 

20 needs differential decoding. 

In Figure 46, a pilot-channel assisted ISR receiver is shown for user d. To better 
understand this structure, it is necessary to develop the corresponding data model a step 
further beforehand. Taking into account the fact that user d uses a pair of spreading 
codes for pilot and data multiplexing with relative powers f 2 and 1, respectively, the 
25 data model writes: 

Y = ya, rf w + Y*4 s ^<t + j + N = Y ^ d b d + + / + N , 

— n — 1),/7 " —0,n n ~n ~n —0,n Tn n —0,n n ~n ~n 

where s ^ d and s *> d denote the data and pilot signal components, respectively. 

The received pilot and data signals can be seen as two separate users received 
from the same physical channel (i.e., H d spread by a pair of codes to yield two spread 

— n 

30 channel versions y 8 ^ and y**A 

In Figure 46, a data ISR beamformer 47V/ l d tuned to the corresponding data 
spreading code outputs the data signal component estimates and interacts with the 
channel identification and the power control units 28V d and 30V d in the regular steps 
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described earlier. However, an additional pilot ISR beamformer 47V/2 d tuned to the 
corresponding pilot spreading code simultaneously provides pilot signal component 
estimates §*' d m Since the pilot power fraction is weaker than the data power, the channel 
identification unit 28V d uses a more reliable and stronger feedback from the data signal 
5 component estimates. The channel estimate is spread by the pair of codes and the 
resulting estimates y 5 ^ and f- d are fed to the data and pilot ISR beamformers, 47 V/ l d 
and 47V/2 d , respectively. The data ISR beamformer 47V/ l d may not need to steer a null 
towards its own weak-power pilot channel. 

On one hand, the data signal component estimates ~ a d ip d b d are derived 

10 within a quantized phase ambiguity a d . They are fed the decision rule unit 29V d to 
estimate the data symbol estimates b d — a d b d within the same phase ambiguity. On 
the other hand the pilot signal component estimates ~ a d \p d n $ hold the constant 
phase ambiguity within little fluctuations of the received power fraction f [^jj . Hence, 
in the same way described above for pilot-symbol assisted ISR, the ambiguity estimator 

15 31V d estimates the phase ambiguity a d n , and forwards it to the conjugator 32V d to 
have {d*}* , then to the multiplier 15 V d for phase compensation of the data symbol 
estimate ~ a d b d . The resulting symbol estimate b d b n is again ambiguity-free 
and no longer needs differential decoding. 

20 Joint ISR Detection 

In the foregoing embodiments of the invention, ISR was applied to a selected set 
D of users, typically users with a low data-rate, who would implement ISR in respect 
of a selected set I of high-rate users. Although this approach is appropriate in most 
cases, particularly when the number of high-rate users is very low, there may be cases 

25 where the mutual interference caused by other high-rate users is significant, in which 
case mutual ISR among high- rate users may be desired as well. Such a situation is 
represented by user sets Ml and M2 of Figure 8. Hence, whereas in the foregoing 
embodiments of the invention, the receiver modules of set I do not perform ISR but 
merely supply constraints sets for use by the receiver modules of set D, it is envisaged 

30 that some or all of the receiver modules in set Ml and M2 also could have beamformers 
employing ISR. Such a Joint ISR (J-ISR) embodiment will now be described with 
reference to Figure 28, which shows only one receiver module, 20* , as an example. In 
any symbol period, each such receiver module 20' (i) receives a constraint 
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matrix C n _ x and an inverse matrix Q n _ 2 and uses them in suppressing interference, 
including its own interference component, and (ii) contributes constraints to the 
constraint matrix c and inverse matrix Q n which will be used in the next symbol 
period. In the case of ISR-H mode receivers, which use hypothetical symbols, it is 
5 merely a matter of replacing the receiver modules in set I with receiver modules 20 d 
having ISR beamformers, since the constraints sets are generated by the hypothetical 
symbols from the bit sequence generators 63 1 ,...,63 NI . Contrary to other ISR modes 
which require decision-feedback, in the ISR-H mode receiver module, no processing 
delay is required for one user to cancel another. Hence, ISR-H can be implemented to 
10 cancel strong interferers without successive interference cancellation or multi-stage 
processing, which will be described later. 

Using Q n and Q n already computed, the ISR combiner for each interferer can 
be obtained readily by: 

W l = C n QR , (71) 

15 

where r = [0,...,0, l,0,...,0] r is a (3NI) -dimensional vector with null components 
except for the k-th one. This implementation has the advantage of implicitly rejecting 
ISI among strong interferers with a single 3NI x 5A7-matrix inversion. 

For the ISR-TR, ISR-R and ISR-D modes, each receiver module, in effect, 

20 combines a receiver module of set I with a receiver module of set D, some components 
being omitted as redundant. Referring again to Figure 28, which shows such a combined 
receiver module, the preprocessor 18 supplies the observation matrix Y n to a 1-bit delay 
45 and a first vector reshaper 44/1, which reshapes the observation matrix Y n to form 
the observation vector X,. A second vector reshaper 44/2 reshapes the delayed 

25 observation matrix Y n _ } to form delayed observation vector Y„_ 7 . These matrices and 
vectors are supplied to the receiver module 20P* and to others of the receiver modules, 
together with the constraint matrix c , and the inverse matrix Q n , from a common 
constraint matrix generator 43P, which generates the constraint matrix C„„ 7 and the 
inverse matrix Q n , 2 from the constraints-set C n _ 2 produced by constraint set generator 42P. 

30 The receiver module 20P 1 comprises a despreader 19 s , a channel identification unit 

28P 1 , a power estimation unit 30P', and a decision rule unit 29P 1 , all similar to those of 
the above-described receiver modules. In this case, however, the receiver module 20P' 
comprises two beamformers, one an ISR beamformer 47P 1 and the other an MRC 
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beamformer 27P 1 , and an additional decision rule unit 29P 72 s which is connected to the 
output of MRC beamformer 27P 1 . The ISR beamformer 47P processes the delayed 
observation vector X,_ 7 to form the estimated signal component estimate and 
supplies it to the first decision rule unit 29P\ the power estimation unit 30P\ and the 
5 channel identification unit 28P 1 , in the usual way. The decision rule unit 29P* and the 
power estimation unit 30P 1 operate upon the signal component estimate to derive 
the corresponding symbol estimate and the power estimate $ and supply them 

to other parts of the receiver in the usual way. 

The despreader 19* despreads the delayed observation matrix Y n , t to form the post- 
10 correlation observation vector z* and supplies it to only the channel identification unit 

n-i 

28P 1 , which uses the post-correlation observation vector z 1 and the signal component 

/i-i 

estimate to produce both a spread channel vector estimate f 1 and a set of channel 

— o,«-i 

vector estimates J-^_ v At the beginning of the processing cycle, the channel 

identification unit 28P' supplies the spread channel vector estimate y* to both the ISR 

— o,«-i 

15 beamformer 47P* and the MRC beamformer 27P 1 for use in updating their coefficients, 
and supplies the set of channel vector estimates j£ to the constraints-set generator 
42P. 

The MRC beamformer 27P 1 processes the current observation vector to 
produce a "future" signal component estimate s i MRCn for use by the second decision rule 
20 unit 29P/2* to produce the "future" symbol estimate b l MRCn , which it supplies to the 
constraints-set generator 42P at the beginning of the processing cycle. The constraints- 
set generator 42P also receives the symbol estimate & 1 from the decision rule unit 29* , 

n-i 

but at the end of the processing cycle. The constraints-set generator 42P buffers the 
symbol £ * from the decision rule unit 29P/2 5 and the symbol estimate h 1 from 

25 the decision rule unit 29P' at the end of the processing cycle. Consequently, in a 
particular symbol period n-1, when the constraints-set generator 42P is computing the 
constraints-set c it has available the set of channel vector estimates the "future" 
symbol estimate £ the "present" symbol estimate b 1 „ and the "past" symbol 
estimate $ 1 „ the latter two from its buffer. 

30 Each of the other receiver modules in the "joint ISR" set supplies its equivalents 

of these signals to the constraints-set generator 42P. The constraints-set generator 42P 
processes them all to form the constraints set c _ and supplies the same to the 
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constraint matrix generator 42P, which generates the constraint matrix c and the 

n 

inverse matrix Q n and supplies them to the various receiver modules. 

The constraints-set generator 42P and the constraint matrix generator 43P will be 
constructed and operate generally in the same manner as the constraints-set generator 43 
5 and constraint matrix generator 42 of the embodiments of the invention described 
hereinbefore with reference to Figures 9 to 27. Hence, they will differ according to the 
ISR mode being implemented. 

When the constraints-set generator 42P of the receiver of Figure 28 is configured 
for the ISR-D mode, i.e. like the constraints-set generator shown in Figure 16, the 
10 constraint matrix c supplied to the ISR beamformer 47P* contains enough information 
for the beamformer 47P 1 to estimate the channel parameters itself. Hence, it forwards 
these estimates to the channel identification unit 28P 1 for use in improving the channel 
vector estimation and the set of channel vector estimates produced thereby. 

An ISR-RH receiver module will use a similar structure, except that the one-bit 
15 delay 45 will be omitted and the constraints-set generator 42P will use the previous 
symbol, estimate the current MRC symbol estimate h l and the two 

n-l' MRC,n 

hypothetical values for "future" symbol to produce current symbol estimate S l n . 
Modification of the receiver module shown in Figure 28 to implement such a "ISR-RH 
mode" will be straightforward for a skilled person and so will not be described hereafter. 
20 In order to implement J-ISR, a more general formulation of the constraint matrix 

is required. The general ISR constraint matrix counting N c constraints, is as follows: 

c , c c N ' 

c J \\c c 

— 11 11 — nj 11 11 — n,N c 11 

25 where the y'-th constraint c is given by: 

C ^ y f / (73) 

— nj ^— ' — k,n 

where Sj defines a subset of diversities which form the y-th constraint when summed. 

30 As shown in Table 2, the sets Sj, j = 7, N c are assumed to satisfy the following 
restrictions: 

S = S l \jS 2 \J.-.[JS N ={(uJ,k)\u = U^M 
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and 



0 being the empty set. Table 1 defines the sets Sj, j=l, N c for all presented ISR 
5 modes of operation. 

The objective signal belongs to the total interference subspace as defined by the 
span of the common constraint matrix c Therefore, to avoid signal cancellation of 
the desired user d by the projection: 

W = / - c o c d " 

10 

the desired-signal blocking matrix (f is introduced, as given by: 



(74) 



C d = 

n 



C 



c 

—1,1 



c 



c d 
c 

—n,N r 



(75) 



15 where: 



—nj ^ — ky 



(76) 



with = {(u,f,k)\u = d; f = i, A^* k = 0J. Normally is a small subset of S 
and c d is very close to C . 



20 



7<?/>if Multi-User Data and Channel Gain Estimation in ISR-D 

Neglecting the signal contributions from the weak-power low-rate users, and 
limiting to the signals of the NI interferers, y 9 can be formulated as: 



25 



NI N f 



(77) 



1=1 /=1 



= C V + N" h , 

"—n —n 



(78) 
(79) 



where r is a N f NI x 1 vector which aligns channel coefficients from all fingers over 

n 

30 all users. Estimation of r , may be regarded as a multi-source problem: 
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t = Q n C" n Y . (80) 

— n n n — n 

This constitutes one step of ISR-D operations and allows joint multi-user channel 
identification. 

5 Multi-Stage ISR Detection 

Multi-stage processing may be used in combination with those of the above- 
described embodiments which use the above-described joint ISR, i.e. all except the 
receivers implementing ISR-H mode. It should be appreciated that, in each of the 
receivers which use decision-feedback modes of ISR (TR,R,D,RH), coarse MRC symbol 

10 estimates are used in order to reconstruct signals for the ISR operation. Because they 
are based upon signals which include the interference to be suppressed, the MRC 
estimates are less reliable than ISR estimates, causing worse reconstruction errors. 
Better results can be obtained by using multi-stage processing and, in successive stages 
other than the first, using improved ISR estimates to reconstruct and perform the ISR 

15 operation again. 

Operation of a multi-stage processing receiver module which would perform 
several iterations to generate a particular symbol estimate is illustrated in Figure 29, 
which depicts the same components, namely constraint-sets generator 42P, constraint 
matrix generator 43P, ISR beamformer 47P 1 and decision rule unit 29P/1 1 , MRC 
20 beamformer 27P' and decision rule unit 29P/2', in several successive symbol periods, 
representing iterations 1 , 2, . . . ,N S of frame n which targets the symbol estimate $ 1 for 

71-1 

user station 10 1 . Iteration 1, if alone, would represent the operation of the receiver 
module 20' of Figure 28 in which the constraints-set generator 42P uses the coarse 
symbol estimates £ 1 previously received from the second decision rule unit 29P/2' 

MKL ,n-\ 

25 (and others as applicable) and buffered. In each iteration within the frame, the other 
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variables used by the constraints-set generator 42P remain the same. These variables 
comprise, from at least each "contributor" receiver module in the same joint processing 
set, the previous symbol estimate b l n _ v the set of channel parameters and the 

current MRC symbol estimate $* Likewise, the spread channel vector 

MRC,n r 

5 estimate y 1 ' and the delayed observation vector Y^^ used by the ISR beamformer 47P 1 

0,n-l 

will remain the same. 

In iteration 1, the constraint matrix generator 42P generates constraint 
matrix C^Q) and the inverse matrix Q n . 2 (l) and supplies them to the beamformer 47F 
which uses them, and the spread channel vector estimate f 4 to tune its coefficients 

~0,n - 1 

10 for weighting each element of the delayed observation vector Y^ l9 as previously 
described, to produce a signal component estimate which the decision rule unit 29V IV 
processes to produce the symbol estimate at iteration 1, which would be the 

n-l v / 

same as that generated by the receiver of Figure 28. This symbol estimate $ 1 m is 
more accurate than the initial coarse MRC estimate b 1 so it is used in iteration 2 

MRC,n-l 

15 as the input to the constraints-set generator 42P 1 , i.e., instead of the coarse MRC 
estimate of beamformer 27P 1 . As a result, in iteration 2, the constraint matrix generator 
42P produces a more accurate constraint matrix C n _ x (2) and inverse matrix Q (2). 
Using these improved matrices, the ISR beamformer 47P* is tuned more accurately, and 
so produces a more accurate symbol estimate 6*^(2) in iteration 2. This improved 

20 symbol estimate is used in iteration 3, and this iterative process is repeated for a total 
of N s iterations. Iteration N s will use the symbol estimate fa 1 X (N -\) produced by the 
preceding iteration and will itself produce a symbol estimate b 1 (N) which is the target 
symbol estimate of frame n and hence is outputted as symbol estimate b ' 

This symbol estimate $* t will be buffered and used by the constraints-set 

71 — 1 ^ 

25 generator 42P in every iteration of the next frame (n + 1) instead of symbol 
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estimate b Other variables will be incremented appropriately and, in iteration 1 of 

frame + a new coarse MRC beamformer 27P 1 symbol estimate b 1 will be used 

MRC,n+l 

by the constraints-set generator 42P. The iterative process will then be repeated, 
upgrading the symbol estimate in each iteration, as before. 
5 It should be noted that, in Figure 29, the inputs to the channel identification unit 

28P 1 use subscripts which reflect the fact that they are produced by a previous iteration. 
These subscripts were not used in Figure 28 because it was not appropriate to show the 
transition between two cycles. The transition was clear, however, from the theoretical 
discussion. 

10 One stage ISR operation can be generalized as follows: 

H(X) = StiKCn ~ ^^(1)^(1), vjtl) = 0„(l)C„(l)»r , (81) 
where s d n (l) is the ISR estimate from first ISR stage, §^ RCn is the MRC signal estimate, 
and the constraint matrices C n (l), C*(l), and 0 (1) are formed from MRC estimates 
at the first stage. Generalizing notation, the signal estimate at stage N s may be derived 
15 after the following iterations: 

SUV = W,„ - W^ C^vp), v_p.) - Q n (2)C„(2)»Y n , 

(82) 

&N S ) = s d MR c. n ~ W^ RCn C d n (N s )v n (N s ), v n (N s ) = Q n {N)C n {NfY n . 

20 The multistage approach has a complexity cost; however, complexity can be 

reduced because many computations from one stage to the next are redundant. For 
instance, the costly computation v(j) could instead be tracked because v (j) *~ v(j-l) if 
the number of symbol estimation errors does not change much from stage to stage, which 
can be expected in most situations. 
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It should be noted that, in the embodiment of Figure 29, the channel identification 
unit 28P 1 updates the channel coefficients after the last iteration N s . Hence, for the next 
cycle, the inputs to the channel identification unit 28P' will be h i z i Y and H l 
It is envisaged, however, that the channel coefficients could be updated more frequently, 
5 conveniently at each iteration. Hence, after the first iteration, the interim symbol 
estimate 6^(1) would be used; in an interim signal component feedback to the channel 
identification unit 28P 1 , after the second iteration, interim symbol estimate &' (2) would 

n-l ^ * 

be used, and so on. The corresponding interim channel estimates for a given iteration 
would be supplied to the constraints-set generator 42P for use in the next iteration. 

10 

ISR using Multistage with Intermediate Channel Decoding (MI CD) 

The TURBO channel encoder has recently attracted researchers interest as a new 
efficient coding structure to achieve Information Bit Error Rates (IBER) close to the 
Shannon limit. Basically, the strength of the TURBO coding scheme is the concatenation 

15 of two convolutional decoders each transmitting the same information data, however, 
data is temporally organized differently (different interleaving) before it is encoded. 
From one of the data streams, the decoder provides likelihood estimates to be used as 
a sort of extrinsic information for decoding of the second stream. 

The TURBO idea has recently been generalized to TURBO multiuser receivers. 

20 Like the TURBO decoder, the TURBO multiuser principle concatenates detection stages. 
This idea applies to ISR, and we will name this extension of ISR Multistage with 
Intermediate Channel Decoding, ISR-MICD. Contrary to multistage ISR, ISR-M, ISR- 
MICD performs channel decoding between stages as illustrated in Figure 47. 

First MRC beamforming 27W 1 is as usual performed in a preliminary stage to 

25 provide coarse signal estimates (s^c J from wh i ch tentative estimates of the transmitted 
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data is derived using the decision rule unit 29W/1 1 . The tentative decision is fed to the 
constraint set generator 42W, which reconstructs the signals of the objective users as 
well as all other interfering users. The constraint matrix generator 43 W assembles the 
constraint matrix (C^U)) to be applied to the ISR beamformer 47W\ The ISR 
5 beamformer 47 W 1 outputs an improved signal estimate (^_j(l)). As usual is, the ISR 
symbol estimate (6^(1)) is computed using the decision rule unit 29W/T, and fed back 
to the constraint set generator 3 42 W for construction of later constraints. More 
importantly, though, the ISR signal estimate is passed on for intermediate channel 
decoding (horizontal branch on figure). First the buffer 90 1 collects n f symbols 

10 corresponding to a code frame 4 . The code frame is passed to the de-interleaving unit 
91 1 which de-interleaves data using a predefined rule 5 . The de-interleaved data is 
Viterbi decoded in the channel decoding unit 92 j to provide a block of NjJR information 
bits where R is the rate of the code. Having performed the channel decoding, the 
decoded information bit sequence is next re-encoded in the re-encoding unit 93', then re- 

15 interleaved in the re-interleaving unit 94' to arrive at the improved channel decoded 
symbol sequence t> l n . This estimate is usually much better than the ISR 
estimate (5 (1)) available before the channel decoding step because redundancy of 
coding and time diversity due to interleaving is exploited 6 . ISR-MICD, as opposed to 
ISR-M, therefore gains from improved estimates in the second stage. As usual 

20 constraints are generated by the constraint set generator 42W and fed to the constraint 



3 This feed back is not explicitly shown on the figure 

4 It is assumed that frame synchronization has been achieved. 

5 Usually block interleaving. 

6 Only at very low SNR, channel decoding does not provide better estimates. 
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matrix generator 43 W to form the improved constraint matrix (c ^2)). This matrix 
is then passed to the ISR beamformer 47W 1 to provide the improved ISR-MICD second 
stage signal estimate (5 n _ 1 (2)) which is fed to decision rule unit 29W/V to provide the 
ISR-MICD second stage data symbol estimate (ft ^2)). ISR-MICD can also be 
5 performed in more stages by repeating the whole process, which amounts to repeating 
the part of the block diagram limited by the broken lines using § (2) for channel 
decoding in the third stage, ^(3) in the fourth stage, etc. 

Group ISR Detection 

10 In practice, the receiver of Figure 28 could be combined with one of the earlier 

embodiments to create a receiver for a "hierarchical" situation, i.e., as described 
hereinbefore with reference to Figure 8, in which a first group of receiver modules, for 
the weakest signals, like those in set D of Figure 8, for example, are "recipients" only, 
i.e., they do not contribute to the constraint matrix at all; a second group of receiver 

15 modules, for the strongest signals, like the receiver modules of set I in Figure 8, do not 
need to cancel interference and so are "contributors" only, i.e., they only contribute 
constraints-sets to the constraint matrix used by other receiver modules; and a third set 
of receiver modules, for intermediate strength signals, like the receiver modules of sets 
M2 of Figure 8, are both "recipients" and "contributors", i.e. they both use the 

20 constraint matrix from the set I receiver modules to cancel interference from the 
strongest signals and contribute to the constraint matrix that is used by the set D receiver 
modules. Generally, this approach is referred to as "Group ISR" (G-ISR) and the 
equations for the constraint matrices and inverse matrices comprising the set 
K = iCU ,£>o C r O 1 use< i by the ISR beamformers in the different 

n l*' Outset ,n* ^Outset, n' ^ Inset, *Zlnset,n) J 

25 receivers are as follows: 
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Qoutset,n = {^Ouset,n^Outset,nj ' (83) 

^■Outset.n ~ ^MmQL-l) ~ ^ Outset, nQ Outset \iA 'Outset. n> (84) 

^Inset.n ~ ^Outset.n^ ' Inset, n' (85) 

Qlnset,n 'inset, 'inset, n) ' (86) 

^ ^Inset.n ~ ^Af*(2L-l) ~ ^ Inset, nQlnset,n^ Inset, n^ 

= nLet, n n 0 utse,,n' ( 88 > 

W = = x ~°" . (89) 

It should be noted that normalization of the columns of r and r is 

Inset, n Outset, n 

10 implicit. 

A receiver module for set D will set II, in Equation (88) to identity which 
means that only "outset" interference will be cancelled. Otherwise, the processing will 
be as described for other receivers of set D. 

A receiver in set Ml does not need to cancel "outset" interference, but does need 
15 to cancel "inset" interference. Consequently, it will set n o in Equation (88) to 
identity so that only inset interference will be cancelled. This corresponds to the joint 
ISR embodiment described with reference to Figure 28. 

Finally , a receiver in set I does not need to cancel any interference . 
Consequently, it will set both n /#lMf and H Qutset to identity, which means that nothing 
20 will be cancelled. This corresponds to the group I receiver modules 20 1 ,...,20 NI 
described with reference to Figures 9, 11, 13, 15-17, 20-24 and 26. 



Successive versus Parallel ISR Detection 

Although the embodiments of ISR receivers described hereinbefore use a parallel 
25 implementation, ISR may also be implemented in a successive manner, denoted S-ISR, 
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as illustrated in Figure 30. Assuming implementation of successive ISR among NI 
interferes, U users, and assuming without loss of generality that users are sorted in 
order of decreasing strength such that user 1 is the strongest and user NI is the weakest 
user, when processing user i in S-ISR, the ISR estimate can be computed as: 

5 ^ = ^-iCc/U«- ".(o-fii^y. (90) 

where £ in spans only the subspace of users i, , /-7 7 , Q l n is the corresponding inverse 
and where c\ n is the user specific constraint matrix. Clearly, c\ n is no longer 
common for all users, which entails expensive matrix inversion for each user. However, 
with ISR-TR this inversion is avoided, since c? "c? is a scalar, and S-ISR-TR is a good 
10 alternative to its parallel counterpart, ISR-TR. Other ISR modes may take advantage of 
the common elements of C in from one processing cycle to the next using matrix 
inversion by partitioning. 

Hybrid ISR Detection 

15 It should also be appreciated that the different ISR modes may be mixed, 

conveniently chosen according to the characteristics of their signals or transmission 
channels, or data rates, resulting hybrid ISR implementations (H-ISR). For example, 
referring to Figure 8, the sets I, Ml and M2 might use the different modes ISR-H, ISR- 
D and ISR-TR, respectively, and the receiver modules in set D would use the different 

20 modes to cancel the "outset" interference from those three sets. Of course, alternatively 
or additionally, different modes might be used within any one of the sets. 



7 And also user i if ISR rejection is desired. 
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ISR Projection for Enhanced Channel Identification 

In all of the above-described embodiments of the invention, the channel 
identification units 28 d in the ISR receiver modules use the post-correlation observation 
vector z' to generate the spread channel vector estimate (by spreading H d )• 

5 Unfortunately, the interference present in the observation matrix Y n is still present in the 
post-correlation observation vector z 1 (see Equation (14)) and, even though it is 

n 

reduced in power by despreading, it detracts from the accuracy of the spread channel 
vector estimate . As has been discussed hereinbefore, specifically with reference 
to Equations (83) to (89), the ISR beamformer 47 d effectively constitutes a 
10 projector and a tuning and combining portion 



r 



r H n d n r 



fl4 

which, in effect, comprises a residual MRC beamformer w d = ~°' n 

Figure 31 illustrates a modification, applicable to all embodiments of the invention 
described herein including those described hereafter, which exploits this relationship to 
improve the spread channel vector estimate (or channel vector estimate H d ) by 
using the projector nf to suppress the interference component from the observation 
20 vector X,. In the receiver of Figure 31, the ISR beamformer 47Q d is shown as 
comprising a projector 100 d and a residual MRC beamformer portion 27Q d . The 
projector 100 d multiplies the projection by the observation vector to produce the 
"cleaned" observation vector ^ d and supplies it to the residual MRC beamformer 

n 

27Q d , which effectively comprises a tuner and combiner to process the "cleaned" 
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observation vector and produce the signal component estimate § d from which 

n n 

decision rule unit 29Q d derives the symbol estimate in the usual way. 

The "cleaned" observation vector ^ d is reshaped by matrix reshaper 102Q d to 

n 

form "cleaned" observation matrix which despreader 19 d despreads to form the 

5 "cleaned" post-correlation observation vector z n,rf for application to the channel 

n 

identification unit 28Q d for use in deriving the spread channel vector estimate 

~~0 t n 

The new "cleaned" vector resulting from the projection of the observation 
vector y by j\ d is defined as follows: 

— n n 

F/ - nfr = n d „ ^ £ ^ y« + jv j = (n'r'X Hn d n N n ) = - . (91) 

10 The new observation vector is free from the interferers and ISI and contains a 

projected version of the channel vector y* 1 Without being a condition, it is reasonable 

— 0,n 

to assume that the projector is almost orthogonal to the channel vector, especially 
in high processing gain situations and/or in the presence of few interferers, and therefore 
consider that Y n - d = Y d . Otherwise an oblique projection can be formed which 

— 0,/7 — 0,n 

15 guarantees = Y d . When despreader 19 d despreads y^ m with the spreading 

sequence of the desired user d, it produces an interference-free projected post-correlation 
observation vector z 11 -" which the channel identification unit 28Q d uses to create the 

— n 

channel vector estimate to use in updating the coefficients of the residual MRC 

—Q,n 

beamformer portion 27Q d . 
20 With respect to the new observation vectors Y* 14 and zF- n before and after 

— n — n ' 

despreading, respectively, the ISR and DFI steps in STAR are modified as follows: 

yd 

W d = -0,n ^ -Q,„ (92) 

-« || 2 r H n d n r 9 

" —0,n 11 —0,n n —0,n 

s" n = Real[w d n H F n "}, (93) 
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H" = ft + JZ** - tfsiSC (94) 
The equivalence between the two expressions of the beamformer coefficients in 
Equation (92) due to the nilpotent property of projections should be noted. In more 
adverse near-far situations, the modification illustrated in Figure 31 allows more reliable 
5 channel identification than simple DFI and hence increases near-far resistance. If 
necessary, this new DFI version will be termed II-DFI. It is expected to be suitable for 
situations where the interferers are moderately strong and when the null constraints cover 
them all. For simplicity of discussion, projection of the observation will become implicit 
without reference to y d , Z d or to the corresponding modifications in STAR-ISR 
10 operations. 

Expanding Dimensionality (X-option) 

When the number of users becomes high compared to the processing gain, the 
dimension of the interference subspace becomes comparable to the total dimension (M(2L 

15 -1)). The penalty paid is an often devastating enhancement of the white noise. Unlike 
ISR-TR, which always requires a single constraint, other DF modes, namely ISR-R and 
ISR-D, may suffer a large degradation because the number of constraints these modes 
require easily becomes comparable to the total dimension available. However, the 
dimension may be increased by using additional data in the observation. This option also 

20 allows for asynchronous transmission and for the application of ISR to Mixed Spreading 
Factor (MSF) systems. 

The matched-filtering observation vector X, is generated to include additional past 
spread data which has already been processed. If the model is expanded to include past 
processed N x symbols and arrive at a total temporal dimension N T = (N X +1)L-1 9 the 

25 observation becomes: 
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Y = 

— ft 



' Y 


U N f 








-n-N x *\ 






— n~N x +l 




= ££ 




+ 
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u = \ f=l 








— n 




— n 




— n 



U N f 

= 5^ J2 ¥ J + n? 



(95) 



« = 1 /=1 



where double underlining stresses the extended model. It should be noted that Y*/ is 

~j 

5 overlapping temporally y*/ and only the first ML samples of the past frames n-1, n-2, 

— y'±i 

.. etc. are used; however the same syntax is used for simplicity of notation. 

As an example, application of the X option to ISR-D, referred to as ISR-DX, 
requires the following constraint matrix: 



10 



C n = 



£.1 
















— n 


" l?"'f" 

— n 


" nr'ii " 

— n 


— n 



(96) 



The extended vectors in Equation (96) have been treated in the same way as those in 
Equation (95), i.e., by concatenating reconstructed vectors from consecutive symbols in 

15 the extended frame and by implicitly discarding overlapping dimensions in the 
concatenated vectors. Clearly, extension of the observation space leaves additional 
degrees of freedom and results in less white noise enhancement. However, it may exact 
a penalty in the presence of reconstruction errors. 

Although the X-option was illustrated in the case of ISR-D, its application to the 

20 remaining DF modes is straightforward. It should also be noted that the X-option allows 
for processing of more than one symbol at each frame while still requiring one matrix 
inversion only. The duration of the frame, however, should be small compared to the 
variations of the channel. 

In the above-described embodiments, ISR was applied to a quasi- synchronous 

25 system where all temporal delays were limited to 0 < r < L. Although this model 
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reflects well the large processing gain situation, where the limit (L -» oo) f allows for 
placing a frame of duration 2L - 1 chips which fully cover one bit of all users, including 
delay spreads. With realistic processing gains, and in particular in the low processing 
gain situation, this model tends to approach a synchronous scenario. Using the X-option 
5 serves as a method supporting complete asynchronous transmission. 

Referring to Figure 32, assuming that the users of the system have processing 
gainL as usual, the transmitted signal of any user is cyclo-stationary and a possible time- 
delay of the primary path r, is therefore 0 < r 7 < L where possible time delays of 
remaining paths are r } < r 2 < . . . < L + Ar where At is the largest possible delay 

10 spread considered. To ensure that the frame covers at least one bit of all users, the 
frame must at least span L + Ar in the despread domain and therefore 2L 4- Ar in the 
spread domain. The observation should be extended slightly beyond that to ease 
interpolation near the edges of the frame. 

Multi-Modulation (MM), Multi-Code (MC), and Mixed Spreading Factor (MSF) 

15 are technologies that potentially can offer mixed-rate traffic in wideband CDMA. MSF, 
which has become very timely, was shown to outperform MC in terms of performance 
and complexity and is also proposed by UMTS third generation mobile system as the 
mixed-rate scenario. Application of ISR to MSF as the mixed rate scenario considered 
herein will now be discussed. 

20 In MSF, mixed rate traffic is obtained by assigning different processing gains 

while using the same carrier and chip-rate. In a system counting two groups of users, 
a low-rate (LR) and a high rate (HR) group, this means that every time a LR rate user 
transmits 1 symbol, a HR user transmits 2r + 1 HR symbols, r = L/L h being the ratio 
of the LR processing gain to HR processing gain. This is illustrated in Figure 33 with 

25 r = 2. 
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Therefore, fitting the ISR frame subject to LR users or in general the lowest-rate 
users ensures that also at least r HR symbols are covered when HR and LR have the 
same delay spread. The ISR generalizes readily to this scenario regarding every HR user 
as r LR users. In Figure 33, the grey shaded HR/LR bits symbolize the current bits to 
5 be estimated; whereas, former bits have already been estimated (ISR-bits) and future bits 
are unexplored. It should be noted that current HR bits should be chosen to lie at the 
end of the frame. 

Multi-Code ISR 

10 It is envisaged that a user station could use multiple codes, N m in number, each 

to transmit a different stream of symbols. Figure 34 illustrates this modification as 
applied to a "without despreading" receiver module 20R d for receiving such a multicode 
signal and using ISR cancellation to cancel interference from other users. The receiver 
module shown in Figure 34 is similar to that shown in Figure 9 except that, instead of 

15 a single ISR beamformer 47 d , the receiver module of Figure 34 has a bank of ISR 
beamformers 47R d * l ,...,47R d,N '" for extracting signal component estimates s*' 1 ,..^^" 1 , 
respectively, and supplying them to a bank of decision rule units 29R d >\. .,29R d ' N % 
respectively which produce a corresponding plurality of symbol estimates S djl \.^b^ m . 
Likewise, the receiver module 20R d has a bank of despreaders 19 rfl 9 ...,19 rf ^ each of 

20 which uses a respective one of the multiple spreading codes of the corresponding user 
d to despread the observation matrix Y n from the preprocessing unit 18 to produce a 
corresponding one of a multiplicity of post-correlation observation 
vectors Z rf ' 1 5 ... 5 2 rf, * m which are supplied to a common channel identification unit 28R d . 

n n 

It should be appreciated that the post-correlation observation vectors share the same 
25 channel characteristics, i.e., of the channel 14 d between user station 10 d and the base 
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station antenna array. Consequently, only one channel identification unit 28R d is 
required, which essentially processes the plural signal component 
estimates s*' 1 ,...,^" 1 and the post-correlation observation vectors and, in essence, 
averages the results to produce a single channel vector estimate H d representing the 

n 

5 physical channel 14 d . The channel identification unit 28R d has a bank of spreaders (not 
shown) which spread the channel vector estimate H d using the multiple spreading codes 

n 

to create a set of spread channel vector estimates y^ 1 yf^m which it supplies to the 
ISR beamformers 47/? rf * 1 ,...,47/? rf,A % respectively. Likewise, the power estimation unit 
30R d is adapted to receive plural signal component estimates s d ' 1 ,...,s„ tNm and essentially 
10 average their powers to produce the power estimate ^ 

While using all of the multiple codes advantageously gives a more accurate 
channel vector estimate, it requires many expensive despreading operations. In order to 
reduce the cost and complexity, the receiver module 20M d may use only a subset of the 
spreading codes. 

15 It can be demonstrated that the multiple spreading codes can be replaced by a 

single spreading code formed by multiplying each of the multiple spreading codes by the 
corresponding one of the symbol estimates b* l ,... 9 6^ Nm an d combining the results. 
Figure 35 illustrates a receiver module which implements this variation. Thus, the 
receiver module 20R d shown in Figure 35 differs from that shown in Figure 34 in that 

20 the bank of despreaders \9 d >\...,\9 dtNm are replaced by a single despreader 19 d - 6 which 
receives the symbol estimates \..,fi d f m and multiplies them by the multiple spreading 
codes to form a compound spreading code, which it then uses to despread the 
observation matrix Y n and form a single post-correlation observation vector z^ 6 . The 

n 

channel identification unit 28R d does not receive the signal component estimates but 
25 instead receives the total amplitude from the power estimation unit 30R d . This 
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serves as a compound signal component estimate because the use of the compound code 
is equivalent to modulating a constant "1" or a constant "-1" with that code, as will be 
formulated by equations later. The channel identification unit 28R d processes the single 
post-correlation observation vector z^ 6 to produce a single channel vector 

n 

5 estimate /y rf and spreads it, as before, using the multiple spreading codes to form the 
multiple spread channel vector estimates f** 1 y 4 ^ for use by the 

— 0,rt v ""' — 0,n 

beamformers 47R d > 1 ,...,47/J d,iV » as before. 

The theory of such multicode operation will now be developed. Assuming for 
simplicity that each user assigned the index u transmits N m streams of DBPSK 
10 data b ul (t),...,b u,Nm (t), using N m spreading codes c uA (t), ...,c u ' Nm (t), each spread stream 
can be seen as a separate user among a total of U X N m access channels, assigned the 
couple-index (u,l). The data model can then be written as follows: 

Y « = E E E + K = EEEE tWifiJtf + K\ ( 97 > 

M =l i=i k=-l u=\ 1=1 f=\ k=-l 

where the canonic w-th user l-th code observation matrices from finger / are 

15 obtained by Equations (3) and (4) of with X(t) in Equation (3) replaced, respectively for 
k = -i,0,+7, 
by: 

X u / f (t) = RJ(t-T p (t)) ® &Ht)c"Xty (98) 
In the equation above, = [0, 0, 1, 0, Of is an M-dimensional vector with 
20 null components except for the m-th one and 5(t) denotes the Dirac impulse. Reshaping 
matrices into vectors yields: 

u=l 1=1 k=-l ' w =l 1=1 f=l k=-l 

The particularity of the above multi-code model, where N m codes of each user 
share the same physical channel h u and the same total received power (\l/ u ) 2 should be 
25 noted. Exploitation of these common features will be discussed hereinafter in relation 
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to adapting of the power-control and the channel-identification procedures to the multi- 
code configuration. The ISR combining step will now be explained. 

Considering first joint ISR combining among the group of N interferers, the 
regular ISR modes, namely TR, R, D, H and RH easily generalize to the new multi-code 
5 configuration of NJJI users instead of A7, as shown in Table 3. ISR combining 
operations are carried out as usual using the constraint and blocking 
matrices £ n and C n l \ respectively. It should be noted, however, that a further 
dimension of interference decomposition and rejection arises over the codes of each user, 
yielding two additional ISR modes. The new modes depicted in Table 3 and referred to 
10 as MCR and MCD (multi-code R and D) characterize interference from the entire set of 
codes of each user by its total realization or by the decomposition of this total realization 
over diversities, respectively. They combine the R and D modes, respectively, with the 
TR mode by summing the corresponding constraints over all the multi-codes of each 
user. 

15 Although these modes partly implement TR over codes, they are still robust to 

power estimation errors. Indeed, the fact that the received power of a given user is a 
common parameter shared between all codes enables its elimination from the columns 
of the constrain matrices (see Table 3). The MCR and MCD modes inherit the 
advantages of the R and D modes, respectively. They relatively increase their sensitivity 

20 to data estimation errors compared to the original modes, since they accumulate symbol 
errors over codes. However, they reduce the number of constraints by N m . 

For a desired user assigned the index d, the constraint matrix c is used to form 

n 

the projector n n . The receiver of the data stream from a user-code assigned the couple- 
index (dj) can simply reject the NI interfering multi-code users by steering a unit 
25 response to j^ ,z and a null response to the constraint matrix r with the 

— Q,n n 
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projector n . It can further reject ISI by steering nulls to y*' and f*- 1 . However, 
.the signals received from other multi-codes contribute to self-ISI. This interference, 
referred to here as MC-ISI, is implicity suppressed when receiving an interfering user. 
It can be suppressed too when receiving the desired low-power user by joint ISR among 
5 the codes of each mobile with any of the ISR modes. The multi-code constraint and 
blocking matrices C d MCn and C^ l Cn , respectively, as shown in Table 4 are formed, and 
derive the ISR beamformer coefficients for user-code (d f l) derived, as follows: 

flfic = {c d M "cAj\ (100) 

W*> 1 = n ^ n . (103) 

The projector n^' that is orthogonal to both MC-ISI and to the NI interferers is formed 
15 and then its response normalized to have a unity response to y*' 

The above processing organization of ISR among the high-power or low-power 
user-codes themselves or between both subsets is a particular example that illustrates G- 
ISR well. The fact that joint ISR among the high-power users and joint ISR among the 
codes of a particular low-power user may each implement a different mode is another 
20 example that illustrates H-ISR well. In the more general case, ISR can implement a 
composite mode that reduces to a different mode with respect to each user. For instance, 
within the group of NT interferers, each user-code assigned the index (i,l) can form its 
own multi-code constraint and blocking matrices C i MCn and Cf^ Cn along a user-specific 
mode (Il n should be set to identity in Table 4). The constraint and blocking matrices 
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then can be reconstructed for joint ISR processing by aligning the individual constraint 
and blocking matrices row-wise into larger matrices as follows: 

& l = \c x Cf' 1 (? J Cf +l Cf 1 1 (105) 
5 This example illustrates the potential flexibility of ISR in designing an optimal 
interference suppression strategy that would allocate the null constraints among users in 
the most efficient way to achieve the best performance/complexity tradeoff. It should 
be noted that, in the particular case where the TR mode is implemented, the matrices in 
Equations (104) and (105) are in fact vectors which sum the individual multi-code 
10 constraint vectors r l and ff* 1 respectively. 

MC,n ^MC,n^ ^ J 

After deriving the beamformer coefficients, each MC user assigned the index u 
estimates its N m streams of data for / = 1, N m as follows (see Figure 34): 

?' 1 = Real{wfY n y (106) 

15 K l = Signfa 1 }, (107) 

and exploits the fact that its N m access channels share the same power, and hence 
smooths the instantaneous signal power of each data stream over all its codes as follows: 

ffln) 2 = (l-a)(tf-i) 2 + *J±^ ( 108 > 



N. 



20 It should be noted that the multi-code data-streams can be estimated using MRC, 

simply by setting the constraint matrices to null matrices. This option will be referred 
to as MC-MRC. 

After despreading of the post-correlation observation vector y by the N m 

— n 

spreading codes of a user assigned the index u, the following post-correlation observation 
25 vectors for / = 7, N are obtained as follows: 
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Z"-' = H u tfjb u n l + N"- 1 = H u s" 1 + A--' . 

—n — n rn —PCM.n — n " —PCM.n 



(109) 



The fact that all user-codes propagate through the same channel is exploited in the 
following cooperative channel identification scheme (see Figure 34): 

H\ = ft + JLTlZ"* 1 - H u s?M, (HO) 
5 which implements a modified DFI scheme, referred to as multi-code cooperative DFI 
(MC-CDFI). MC-CDFI amounts to having the user-codes cooperate in channel 
identification by estimating their propagation vectors separately, then averaging them 
over all codes to provide a better channel vector estimate. It should be noted that 
implicit incorporation of the II-DFI version in the above MC-CDFI scheme further 
10 enhances channel identification. 

Since the STAR exploits a data channel as a pilot, it can take advantage of a 
maximum of N m expensive despreading operations. To limit their number in practice, 
MC-CDFI can be restricted to a smaller subset of 1 to N m user-codes. A compromise 
can be found between channel estimation enhancement and complexity increase. 
15 Another solution that reduces the number of despreading operations reconstructs 

the following data-modulated cumulative-code after ISR combining and symbol 
estimation in Equations (106) and (107): 

1 

«,5 1 £ u -l u,l /1 1 1\ 

C * = M-22 b n C lc - OH) 



20 A single despreading operation with this code yields: 



Z"- a = \p u „H u 
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* •> 

4-? —PCM,n 
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r "^—n — PCM,n 



(112) 
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It has the advantage of further reducing the noise level by N m after despreading, while 
keeping the signal power practically at the same level 8 . The data-modulated cumulative- 
code can be used to implement channel identification as follows: 



This CDFI version is referred to as S-CDFI (see Figure 35). 

Whereas multicode operation involves user stations transmitting using multiple 
spreading codes, but usually at the same data rate, it is also envisaged that different users 
within the same system may transmit at different data rates. It can be demonstrated that 
10 the receiver modules shown in Figures 34 and 35 need only minor modifications in order 
to handle multirate transmissions since, as will now be explained, multicode and 
multirate are essentially interchangeable. 

Multi-Code Approach to Multi-Rate ISR 
15 Reconsidering now the conventional MR-CDMA, in this context, STAR-ISR 

operations previously were implemented at the rate IIT where T is the symbol duration. 

As described earlier, with reference to Figure 32, the "X option" extensions, enables 

reduction of noise enhancement by increasing the dimension of the observation space and 

provides larger margin for time-delay tracking in asynchronous transmissions. A 
20 complementary approach that decomposes the observation frame into blocks rather than 

extends it using past reconstructed data will now be described. 

This block-processing version of STAR-ISR will still operate at the rate IIT on 

data frames barely larger than the processing period T However, it will decompose 




— n + 1 — n r \ — n 



(113) 



There is a small power loss due symbol estimation errors (very low in practice). 



83 

each data stream within that frame into data blocks of duration T r where T r is a power-of- 
2 fraction of T. The resolution rate 1/T r can be selected in the interval [l/7\l/r c ]. 
Hence, a receiver module that processes data frames at a processing rate 1/T with a 
resolution rate 1/T r can only extract or suppress data transmissions at rates slower than 
5 or equal to \IT r . Also, the channel parameters of the processed transmissions must be 
almost constant in the interval T, the processing period. This period should be chosen 
to be much larger than the delay spread At for asynchronous transmissions, but short 
enough not to exceed the coherence time of the channel. 

In one processing period, STAR-ISR can simultaneously extract or suppress a 
10 maximum N m = T/T r blocks (N m is a power of 2). In the n-th processing period of 
duration T, a stream of data b u (t) yields N m samples b^\... 9 b" ,Nm sampled at the 
resolution rate. Over this processing period, therefore, the spread data can be developed 
as follows: 

N m 

c\t)b%i) = XXtfcC* - nT)C(t), (114) 

1=1 

15 where |J^(r) is the indicator function of the interval [(/ - 1)7;, lT r ). This equation can 
be rewritten as follows: 

N m 

c u (t)b\i) = Y^b ul {t)c ul {t), (115) 
where b u l (t) 9 .„,b u,Nm (t) represent N m data-streams at rate 1/T spread by N m virtual 
orthogonal codes c w/ (f ),..., c"^ m (0 (see Figure 36). 

20 With the above virtual decomposition, one arrives at a MC-CDMA model where 

each of the processed users can be seen as a mobile that code-multiplexes N m data- 
streams over N m access channels. This model establishes an equivalence between MC- 
CDMA and MR-CDMA and provides a unifying framework for processing both 
interfaces simultaneously. In this unifying context, codes can be continuous or bursty. 

25 Use of bursty codes establishes another link with hybrid time-multiplexing CDMA (T- 
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CDMA); only the codes there are of an elementary duration T r that inserts either symbols 
or fractions of symbols. A larger framework that incorporates MR-CDMA, MC- 
CDMA, and hybrid T-CDMA can be envisaged to support HDR transmissions for third 
generation wireless systems. 
5 Exploiting this MC approach to MR-CDMA, the data model of MR-CDMA will 

be developed to reflect a MC-CDMA structure, then a block-processing version of 
STAR-ISR derived that implements estimation of a symbol fraction or sequence. 

The multi-code model of Equation (97) applies immediately to MR-CDMA. 
However, due to the fact that codes are bursty with duration T r < T, the self-ISI 

10 vectors f 14 ' 1 and p* 1 and the spread propagation vector y"'' of a given user-code do 
not overlap with each other. If 27 denotes an arbitrarily enlarged delay-spread 
(reference [20]) to leave an increased uncertainty margin for the tracking of time- varying 
multipath-delays (i.e., Ar < Sr < 7), and if jsf r = fSr/r"] denotes the maximum 
delay-spread in T r units, then only the last N r symbols b"^ m ~ N ' +1 9 ... 9 b"^ m among the past- 

15 symbols in the previous frame may contribute to self-ISI in the current processed frame 
(see Figure 37): 



u 



Y n = E 



« = 1 



E VKtxin + E + e 

i=n-n+\ i=\ i=\ 



\ + <\ (116) 



20 In this frame of duration 2T - T c , the desired signals' contribution from the N m current 
symbols is contained in the first interval of duration T + St, whereas the remaining 
interval of the frame contains non-overlapping interference from the last N m - N r future 
symbols in the next frame, namely b"^' +1 ,...,b"^ m (see Figure 37). The remaining part 
of the frame can be skipped without any signal contribution loss from the current bits. 

25 Hence, the duration of the processed frame can be reduced to y + ~Sr - T as follows: 
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y n = Ko.^p--w 2 ]' (117) 

where L A = [ ~KtIT^\ is the maximum length of the enlarged delay-spread in chip 
samples. With the data block-size reduced to M x (L + L A - 1), the matched-filtering 
5 observation matrix reduces to: 



u 



W=l 



l=N -N+\ 1=1 l=l 



\ + *f , (118) 



where tf% h is the noise matrix reduced to the same dimension. This data model 
10 equation can be rewritten in the following compact vector form: 

z„ = EEE Mtetfi (H9) 

«=1 1=1 k=-l 



where \[ = 0 if k = -1 and / G {l,...,N r } or if* = +7 and / G {N m - N r + 1,..., 
N„}, and 1 otherwise. 

15 The constraint matrices can be formed in an MC approach to implement joint or 

user-specific ISR processing in any of the modes described in Tables 3 or 4, 
respectively. In contrast to the conventional MC-CDM A, the factor \ l k discards all non- 
overlapping interference vectors in the processed frame and somehow unbalances ISI 
contribution among the virtual multi-code streams. In the DF modes, only the central 

20 streams of each user (i.e., I = N r + 1, N m - N r ) sum symbol contributions from the 
previous, current and future symbols; whereas the remaining streams sum signal 
contributions from either the current and the previous or the current and the future 
symbols. Indeed, the 2(N m - N r ) ISI terms discarded from summation contribute with 
null vectors to the processed frame. In the ISR-H mode, the columns previously 

25 allocated to individually suppress these vectors are eliminated from the constraint 
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matrices, thereby reducing the number of its columns to N H = N m + 2N r constraints per 
user 9 (see Tables 3 and 4). ISR-H hence approaches ISR-R in computational complexity 
when N r is small compared to N m . 

After derivation of the beamformer coefficients of each virtual user-code assigned 
5 the couple-index (u,l) f its signal component s £ l is estimated using Equation (106). In 
this process, each ISR combiner rejects the processed interferers regardless of their exact 
data rates, which only need to be higher than the resolution rate. This feature finds its 
best use when implementing ISR at the mobile station on the downlink where data rates 
of suppressed interferers are not necessarily known to the desired mobile-station. For 

10 instance, orthogonal variable spreading factor (OVSF) allocation of Walsh spreading 
codes is no longer necessary. On the uplink, each transmission rate is known to the base 
station. However, one can still gain from this feature by allowing joint and well 
integrated processing of mixed data traffic at a common resolution rate. 

Indeed, the estimation of the signal components provides sequences oversampled 

15 to the resolution rate 1/T r Hence, after a given data stream is decomposed at this 
common rate, its signal component estimate must be restored to its original rate in an 
"analysis/synthesis" scheme. To do so, the data rate 1/T U < 1/T r of user u is defined 
and it is assumed temporarily that it is faster than the processing rate (i.e., 1/T U > 1/T). 
Hence, one can extract from each frame F u = T/T u < N m signal component estimates out 

20 of N m by averaging the oversampled sequence s uJ over consecutive blocks of size B„ = 
N m /F u = TJT r for n f = 0,..., F u - 1 as follows: 



9 ISR may be equally reformulated with N m + 2N T generating sequences that process 
all the contributing symbols as if they were independent streams without MC-ISI. Only 
25 the N m current symbols are estimated then; the 2N T remaining symbols being corrupted 
by the edge effect. 
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(fl' + l)B, 

E # 

B 



^ + „, = (120) 



u 



= M^.-'}' (121) 



= (1 - a)0 + • ( 122 > 

F u 



In the particular case where the data rate is equal to the processing rate (i.e., 1/T U = 
10 7/7), the equations above have simpler expressions with F u = 7 and i? w = Af m : 

E 

j; = ±» d23) 
N m 

S u n = S/^}, (124) 

15 

(^) 2 = (1 " «)(C-i) 2 - « |^| 2 . (125) 

If the data rate is slower than the processing rate, the signal component 
estimate § u n of Equation (123) is further averaged over consecutive blocks of 
20 size p = J IT to yield the following subsampled sequence: 

^In/F'jFi+n' 

e" - (126) 

Symbol and power estimations in Equations (124) and (125) are on the other hand 
25 modified as follows: 
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^Uj = Sign{sl /F .j}, 



(127) 



K*if - <* - «>(*w.) 2 + « i^j i 2 - < 128 > 

5 It should be noted that a higher value is needed for the smoothing factor a to adapt to 
a slower update rate of power estimation. If the channel power variations are faster than 
the data rate, then it is preferable to keep the power estimation update at the processing 
rate in Equation (125). In this case, Equation (126) is modified as follows 10 : 



E7 u $ u 
V ln/Fj F' u+ n' S [n/Fj F u + n f 

10 = , (129) 



F a -l 



5^ foln/FjF^n') 



to take into account channel power variations within each symbol duration. 

It should be noted that the multi-rate data-streams can be estimated using MRC, 
15 simply by setting the constraint matrices to null matrices. This option may be referred 
to as MR-MRC. 

It should be also noted that combination of Equations (106) and (120), along with 
Equation (128) for data rates slower than the processing rate, successively implements 
the processing gain of each user in fractioned ISR combining steps. 
20 In general, regrouping the symbol-fractions back to their original rate can be 

exploited in the design of the constraint matrices; first by reducing reconstruction errors 
from enhanced decision feedback; and secondly by reducing the number of constraints 



10 This signal component estimate is not used for power estimation. Only its sign 
is taken in Equation (127) as the estimate of the corresponding bit. Hence, power 
25 normalization given here for completeness is skipped in practice. 
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of a given user u from N m to F u in the modes implementing decomposition over user- 
codes {i.e., R, D, and H). For these modes, the common factor A^V/, appearing in the 

NI 

total number of constraints N c reduces to ^ p., by regrouping the constraint vectors 

1=1 

over the user-code indices that restore a complete symbol within the limit of the 
5 processing period 11 . 

Regrouping the constraints of user u to match its original transmission rate 
amounts to regrouping the codes of this user into a smaller subset that corresponds to a 
subdivision of its complete code over durations covering its symbol periods instead of 
the resolution periods. In fact, user u can be characterized by F u concatenated multi- 

10 codes instead of N m . Overall, MR-CDMA can be modeled as a mixed MC-CDMA 
system where each user assigned the index u has its own number F u of multi-codes (see 
Figure 38). Therefore, the ISR-combining and channel- identification steps can be carried 
out in one step along the MC formulation of the previous section, using user-codes 
simply renumbered from 1 to F u for simplicity. Hence, as shown in Figure 39, the only 

15 change needed to the receiver module of Figure 34 is to the bank of despreaders. In the 
receiver module shown in Figure 34, the spreading codes used by the 
despreaders \9 d ^...,\9 d,Fd comprise segments of the spreading code of user d, i.e., the 
segments together form the part of the code used in a particular frame. The number of 
code segments F u corresponds to the number of symbols b u / ,...,b u „ Fu transmitted in the 

20 frame. The estimates of these symbols, and the signal component 
estimates ^\...J u n F " map with those of Equations (120), (121), (123) and (124) within 
a parallel/serial transform. 



11 Feedback of symbols with rates slower than the processing rate to the constraints- 
set generator is feasible. 
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This illustrates again the flexibility afforded by using ISR in designing optimal 
interference suppression strategies that suit well with MR-CDMA. It enables 
simultaneous processing of blocks of symbols or fractions of symbols in an integrated 
manner at two common resolution and processing rates. 

To carry out channel identification operations, the M X L reduced-size post- 

A 

correlation observation matrix of user-code (u,l) is defined as follows: 



(130) 



where the columns of this matrix are given for 7 = 0,...,L A - 1 by: 

z "j = Y njv4' 1 = y £ (131 > 

10 

This correlation with the virtual user-code (u,l) amounts to partial despreading by a 
reduced processing gain l = TIT = LIN , using the l-th block of length L r of the 

r re m ' r 

user's code Cj 1. It should be noted that, in contrast to conventional MC-CDMA, the 
above partial despreading operations are less expensive in terms of complexity per user- 
15 code. 

The reduced-size post-correlation observation vector z*- 1 resulting from vector- 

n 

reshaping of j?* has the same model expression of Equations (109), except that vectors 
there all have reduced dimension (ML A ) x 1. It should be noted that the post- 
correlation window length L A was fixed long enough to contain the delay-spread with an 
20 enlarged margin for asynchronous time-delay estimation from the reduced-size 
propagation vector u u (reference [20]). Identification with post-correlation windows 

— n 

shorter thanL, investigated in [6], reduces complexity and proves to work nearly as well 
as the original full-window version of STAR (Le, L A = L). 

Channel identification with the MC-CDFI scheme of Equation (110) can be 
25 readily implemented using the user-code post-correlation observation vectors z? J . 
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However, this procedure would feed back symbol fractions without taking full advantage 
of the complete processing gain. Instead, the vectors Z"-' are regrouped and averaged 

n 

in the same way the signal component estimates are restored to their original rate in 
Equations (120), (123) or (126), and z w , Z u or z? ■ respectively 12 are obtained. 
5 Hence, the CDFI channel identification procedure, renamed MR-CDFI, is implemented 
as follows: 



F-l 



when the data-rate is faster than the processing rate 13 , or by: 

10 ft = ft + JZ U - frsftSl (133) 

in the particular case where the data rate is equal to the processing rate, or by: 

S-lnlF'J = KnlFA + ^ ^Uj " ) • (134) 

when the data-rate is slower than the processing rate. It should be noted that channel 
identification at data rates faster than the processing gain in Equation (132) has a 

15 structure similar to MC-CDFI. Averaging over F u despread observations there can be 
reduced to a smaller subset to gain in complexity like in MC-CDMA. Use of the 5- 
CDFI version described in Equations (111) to (113) instead of, or combination with, the 
above scheme are other alternatives that reduce the amount of complexity due to 
despreading operations. 

20 By regrouping codes to match the original data transmission rates as discussed 

earlier (see Figure 38), channel identification can be easily reformulated along a mixed 



12 In practice, these vectors are computed directly from F n in regular despreading 
steps which exploit the entire spreading sequences in one step along a mixed MC-CDMA 
scheme. 



25 ^Implementation of F u channel updates (with time-delay tracking) instead of 
averaging is computationally more expensive. 
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MC-CDMA model where each user is characterized by F u multi-codes and F u despread 
vectors Z w l ,... J Z" ,F % as shown in Figure 39. 

n n 

To reduce further the number of expensive despreading operations, slower 
channel identification (reference [20]) can update channel coefficients less frequently if 
5 the channel can still show very weak variations over larger update periods. However, 
high mobility can prevent the implementation of this scheme and faster channel 
identification update may even be required. For data rates faster than the processing 
rate, updating at a rate higher than the processing rate is not necessary. The processing 
period T is chosen to guarantee that the channel parameters are constant over that time 

10 interval. For data rates slower than the processing rate, the channel update rate could 
be increased above the data rate up to the processing rate using Equation (133) and 
partial despreading to provide Z". In Equation (133), $1 , from Equation (126) 
should be fed back instead of § u n to benefit from the entire processing gain in the 
decision feedback process. 

15 Although the foregoing embodiments of the invention have been described as 

receiver modules for a base station, i.e., implementing ISR for the uplink, the invention 
is equally applicable to the downlink, i.e., to receiver modules of user stations. 

Downlink ISR 

20 Figure 48 illustrates how the downlink can be modelled like an uplink, so that the 

ISR techniques developed for the uplink can be employed. Figure 48 shows a single 
(desired user) mobile 10 d which will be one of many, receiving signals from all base 
stations within range. Only the serving base station ll v and three main interfering base 
stations ll 1 ,..., ll u , ... 11 M are shown. The mobile station may communicate with two 

25 or three base stations at a given time, so two of the interferers could be such other 
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stations. The bold arrow illustrates that only one signal, the desired user signal d from 
server base station ll u is the signal which is to be received. The serving base station 11" 
will also transmit to other mobiles and those signals will also be received by mobile 10 d 
as interfering signals. Also, each of the other three base stations will be transmitting. 
5 Even though one of the signals from each of the other base stations are intended for the 
mobile d, they constitute interference so far as the reception of signal 10 d from serving 
base station ll u is concerned. Basically, signals #1 to #NI of BS IV (excepting d) are 
in-cell interferers and the interfering signals from the other base stations are out-cell 
interferers. 

10 The base stations range from #1 to # NB. BS \\ v is a generic one of them. BS 

11° is a specific one of the NB stations is the serving base station. 

There will be other signals from other base stations, some of which are shown 
in broken lines, and other signals also shown in broken lines, from the base stations 
shown in full, since figure 48 shows only the strongest signals transmitted by each base 

15 station. I.e. , there will be more than NI mobile stations in cell v but their transmissions 
weaker. Those mobiles close to the base station will require weak transmissions whereas 
those far from the base station will require more power, and the base station power 
control will increase transmit power to achieve it. Also, the data rates could vary and 
hence affect power levels. Consequently, many signals are not represented in Figure 48. 

20 They are ignored because they are relatively weak. Of course, they will be part of the 
noise signal represented in Figure 2. 

Each base station may be transmitting multiple codes, so its signal will have a 
multi-code structure. Because these signals are being transmitted via the same antenna, 
they are similar to a multicode signal. Hence, mobile station 10 d receives from base 

25 station ll v the signals #1 ...d...NI via channel v, so those signals appear to be a 
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multicode signal. The same applies to signals received from other base stations for 
processing. 

The "zoom" inset shows the signal transmitted to a mobile station i by its serving 
base station v' but received as interference by mobile 10 d . It could be a multirate and/or 
5 multicode signal, since it is a summation of different signals using different codes and 
possibly different rates as described with reference to Figures 40-42, for example, though 
their despreaders handle all interferers whereas here we select only one. The composite 
signal comprises components transmitted from base station v' to user / using codes 
ranging from 1 to F h i.e., from first component #(v', /, 1) to a final component (v', i, 

10 F,). Hence the multicode structure model developed for the uplink applies to the 
downlink as well at the mobile receiver. 

To implement ISR rejection, the user/mobile station needs to identify the group 
of users {i.e., interferers) to suppress. Assuming temporarily that suppression is 
restricted to in-cell users, served by base-station v, and that the number of suppressed 

15 interferers is limited to NI to reduce the number of receivers needed at the desired base- 
station to detect each of the suppressed users, in order to identify the best users to 
suppress, the user station can probe the access channels of base-station v, seeking the M 
strongest transmissions. Another scheme would require that the strongest in-cell 
interfering mobiles cooperate by accessing the first NI channels (i.e., u = i G 

20 {1, . . .,NI}) of base-station v. 

Once the NI suppression channels have been identified, the desired user-station 
can operate as a "virtual base-station" receiving from A7 "mobiles", i.e., base station 
transmitter modules, on a "virtual uplink". If the desired user is not among the NI 
interferers, an additional user station is considered. Similar NI channels may be 

25 identified for transmissions from the neighbouring base-stations. Accordingly, 
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consideration will be given to the NB base-stations,. assigned the index v' G {1,...,NB}, 
which include the desired base-station with index v' = v without loss of generality. This 
formulation allows the user-station to apply block-processing STAR-ISR with specific 
adaptations of ISR combining and channel identification to the downlink. 
5 In essence, each "virtual base station" user station would be equipped with a set 

of receiver modules similar to the receiver modules 20 1 ,...,20 u , one for extracting a 
symbol estimate using the spreading code of that user station and the others using 
spreading codes of other users to process actual or hypothesized symbol estimates for the 
signals 'of those other users. The receiver would have the usual constraints-set generator 

10 and constraint matrix generator and cancel ISR in the manner previously described 
according to the mode concerned. 

It should be appreciated, however, that the signals for other users emanating from 
a base station are similar to multicode or multirate signals. Consequently, it would be 
preferable for at least some of the user station receiver modules to implement the 

15 multicode or multirate embodiments of the invention with reference to Figures 34 and 
39. Unlike the base station receiver, the user station's receiver modules usually would 
not know the data rates of the other users in the system. In some cases, it would be 
feasible to estimate the data rate from the received signal. Where that was not feasible 
or desired, however, the multirate or multicode receiver modules described with 

20 reference to Figures 34 and 39 could need to be modified to dispense with the need to 
know the data rate. 

Referring to Figure 40, the user station receiver comprises a plurality of receiver 
modules similar to those of Figure 39, one for each of the NB base stations whose NI 
strongest users' signals are to be cancelled, though only one of them, receiver module 
25 20 v ' for base station v', is shown in Figure 40. Recognizing that one or more of those 
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NI signals could be multirate or multicode, and hence involve not only different 

i=NI 

spreading codes but also different code segmentations, the number of despreaders equals p., 

1=1 

i.e. 19 V ' ' \„, 19 V ' 1F ',..., 19 v '* u ,..., 19 V '-'' F ',..., 19 v ' i9*w» In any 
given base station, the NI users are power-controlled independently and so are received 
5 by the mobile/user station with different powers. Consequently, it is necessary to take 
into account their power separately, so the power estimates from power estimation means 
30T V ' are supplied to the channel identification unit 28T V '. The channel identification unit 
28T V ' processes the data in the same way as previously described, spreading the resulting 
channel vector estimate /y' to form the spread channel vector 

— 0,n 

10 estimates y*' ,u ,.. , f w « and supplying them to the ISR 

—0,n ' — 0,/i 

beamformers 477 ,v ' ,u ,..., 477 TV '* MF ^ ) respectively, foruse in processing the observation 
vector y 

n * 

O 

The resulting signal component estimates S^ ni,Fn ' are similarly fed 

back to the channel identification unit 28T V ' to update the channel vector estimates and ' 
15 to the decision rule units 30T v ' l \ 30T V ' ,A7,Fw for production of the corresponding 
symbol estimates ft*' 1,1 ,..., ^' A7F -. In all modes except ISR-H, these symbol 
estimates are supplied to the constraints-set generator (not shown in Figure 40), together 
with the set of channel vector estimates < y^ 9 from channel identification unit 28T V ' for 
use by the constraint matrix generator (not shown) in forming the set of constraints c . 

n" 

20 The set of channel parameter estimates includes the power estimates from the power 
estimation units. It should be noted that the constraints-set generator and constraint 
matrix generator may be described hereinbefore, the actual configuration and operation 
being determined by the particular ISR mode selected. 

If the desired user, i.e. of the user station receiver 20 v , is among the NI users, 

25 its symbols will be extracted and output to the subsequent parts of the receiver in the 
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usual way. If, however, the desired user is not among the NI strong users of the serving 
base station v, the user station receiver will include not only one of the receiver modules 
of Figure 40 for each base station but also a separate receiver module specifically for 
extracting the signal for the desired user and which could be similar to that shown in 
5 Figure 39. 

Bearing in mind, however, that the channel vector estimate derived by the 
receiver module for the serving base station's strong users in Figure 40 will be for the 
same channel, but more accurate than the estimate produced by the channel identification 
unit of Figure 39, it would be preferable to omit the channel identification 

10 unit (29 F *) and despreaders 19* 1 ,..., 19 d,F " (Figure 39), and supply the spread 
channel vector estimates from the channel identification unit of the receiver module for 
serving base station v. As shown in Figure 41, in such a modified receiver module 
20T d , the beamformer means comprises a bank of ISR beamformers 
(47T v > d > 1 y ... 9 47T v,d ' Fd ), , a bank of decision rule means (29T v ' d '\.„ 9 29T v,d ' Fd ), and 

15 power estimation means 30T v d which process elements of the observation vector Y n from 
vector reshaper 44 in a similar manner to the receiver module of Figure 39. The bank 
of beamformers (47r v ^ 1 ,... 5 47T v ' d ' Fd ), are tuned, however, by a set of channel vector 
estimates (y^' 1 ,..., Y^^), produced by the channel estimation means (28T V ) (Figure 
40) corresponding to serving base station v, which produces the channel vector estimate 

20 h v and spreads it to produce the channel vector estimates ( y^ 1 x?> d * F d\ . 

The receiver module shown in Figure 40 is predicated upon the data rates of each 
set of NI users being known to the instant user station receiver. When that is not the 
case, the receiver module shown in Figure 40 may be modified as shown in Figure 42, 
i.e., by changing the despreaders to segment the code and oversample at a fixed rate that 

25 is higher than or equal to the highest data rate that is to be suppressed. 
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It is also possible to reduce the number of despreading operations performed by 
the receiver module of Figure 42 by using a set of compound segment codes as 
previously described with reference to Figure 35 to compound over segments. However, 
as shown in Figure 43, a set of different compound codes could be used to compound 
5 over the set of NI interferers. It would also be possible to combine the embodiment of 
Figure 43 with that of Figure 35 and compound over both the set of interferers and each 
set of code segments. 

A desired user station receiver receiving transmissions on the downlink from its 
base-station and from the base-stations in the neighbouring cells will now be discussed. 

10 Each base-station communicates with the group of user stations located in its cell. 
Indices v and u will be used to denote a transmission from base-station v destined for 
user w. For simplicity of notation, the index of the desired user station receiving those 
transmissions will be omitted, all of the signals being implicity observed and processed 
by that desired user station. 

15 Considering a base-station assigned the index v, its contribution to the matched- 

filtering observation vector y of the desired user station is given by the signal vector 

n 

of the v-th base-station Y" defined as: 

— n 

yv = y T- u , (135) 

— u,n — u,n 

u = \ 

20 

where the vector denotes the signal contribution from one of the U v users 

— u,n 

communicating with base-station v and assigned the index u. Using the block-processing 
approach described in the previous section, the vector y.« can be decomposed as 

— u,n 

follows: 

25 T- = £ btf&JgK. (136) 

l=k k=-\ 
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It should be noted that the channel coefficients f fn just hold the index of the base-station 
v. Indeed, transmissions from base-station u to all its mobiles propagate to the desired 
user station through a common channel. Base-station signals therefore show a multi-code 
structure at two levels. One comes from the virtual or real decomposition of each user- 
5 stream into multiple codes, and one, inherent to the downlink, comes from summation 
of code-multiplexed user- streams with different powers. As will be described 
hereinafter, this multi-code structure will be exploited to enhance cooperative channel 
identification at both levels. 

In a first step, the desired user-station estimates the multi-code constraint and 
10 blocking matrices of each of the processed in-cell users (i.e., u G {1,...,NI} U {d}). 
Table 4 shows how to build these matrices, renamed here as £ v,u and r v,uJ to show 

MC,n MC,n 

the index v of the serving base-station. Indexing the symbol and channel vector 
estimates with v in Table 4 follows from Equation (136). In a second step, the user- 
station estimates the base-specific constraint and blocking matrices c v and C v ' aJ using 

BS,n BS,n 

15 Table 5. These matrices enable suppression of the in-cell interferers using one of the 
modes described in Table 5. For the downlink, a new mode BR, for base-realization, 
replaces the TR mode of Table 3. Suppression of interfering signals from multiple base- 
stations adds another dimension of interference decomposition and results in TR over the 
downlink as shown in Table 6. Therefore, in a third step the mobile-station estimates 

20 the base-specific constraint and blocking matrices c v an d C v w,/ from the interfering 

BS t n BS,n 

base-stations and concatenates them row-wise to form the multi-base constraint and 
blocking matrices denoted as £ n and C v ^\ respectively. In the TR mode, the base- 
specific constraint and blocking vectors in the BR mode now are summed over all 
interfering base-stations, leaving a single constraint. For the other modes, the number 
25 of constraints AT. in Table 3 is multiplied by the number of interfering base-stations NB. 
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The receiver module dedicated to extracting the data destined to the desired mobile- 
station # d from the serving base station # v is depicted in Figure 41. 

It should be noted that the multi-rate data-streams can be estimated using MRC 
on the downlink, simply by setting the constraint matrices to null matrices 14 . This 
5 option will be termed D-MRC. 

If the user-station knows 15 the data rates of the suppressed users, it can estimate 
their symbols 16 as long as their symbol rate does not exceed the processing rate. As 
mentioned hereinbefore, this block-based implementation of the symbol detection 
improves reconstruction of the constraint matrices from reduced decision feedback 

10 errors 17 . Otherwise, the user-station can process all interfering channels at the common 
resolution rate regardless of their transmission rate. It should be noted that estimation 
of the interferers' powers is necessary for reconstruction in both the BR and TR modes, 
for channel identification as detailed below, and possibly for interference-channel probing 
and selection. It is carried out at the processing rate. 

15 Identification of the propagation channels from each of the interfering base- 

stations to the desired user station is required to carry out the ISR operations. 
Considering the in-cell propagation channel, its identification from the post-correlation 
vectors of the desired user is possible as described hereinbefore with reference to Figure 

14 In this case, ISR processing is not needed and the desired signal is expected to be 
20 strong enough to enable reliable channel identification for its own. 

15 Data-rate detection can be implemented using subspace rank estimation over each 
stochastic sequence of N m symbol fractions. 

16 In the ISR-H mode, only the signal component estimates are needed for power and 
channel estimation (see next subsection). 

25 17 Recovery of the interfering symbols at data rates slower than the processing gain 
could be exploited in slow channel identification. However, selection of a user as a 
strong interferer suggests that its transmission rate should be high. 
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39. It exploits the fact that the multi-codes of the desired user propagate through the 
same channel. However, the in-cell interfering users share this common channel as well. 
Therefore, the MC-CDFI and MR-CDFI approaches apply at this level as well. Indeed, 
the user-station has access to data channels which can be viewed as NI x N m virtual 
5 pilot-channels with strong powers. It is preferable to implement cooperative channel 
identification over the interfering users whether the desired user is among the in-cell 
interferers or not. The same scheme applies to the neighbouring base-stations and 
therefore enables the identification of the propagation channel from each out-cell 
interfering base-station using its NI interfering users. 
10 If the data rates are known to the base-station, identification of the propagation 

channel from a given base-station v' G {!,..., NB) can be carried out individually from 
each of its NI interfering users, as described in the previous section. To further enhance 
channel identification, the resulting individual channel vector estimates are averaged over 
the interfering users. Both steps combine into one as follows: 



This downlink version of MR-CDFI, referred to as DMR-CDFI, is illustrated in Figure 
40. It should be noted that averaging over the interferers takes into account 
normalization by their total power. To reduce the number of despreading operations, 
20 averaging over interferers can be limited to a smaller set ranging between 1 and NI. 

If the data rates of the interfering users are unknown to the user-station, 
identification can be then carried out along the steps described with reference to Figure 
34 to process interfering signals at the common resolution rate as follows: 




F.-l 



15 



(137) 



NI AT 




(138) 
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This downlink version of MC-CDFI, referred to as DMC-CDFI, is illustrated in Figure 
42. To reduce the number of despreading operations, averaging over interferers and 
user-codes can be limited to smaller subsets ranging between 1 and NI and 1 and N m9 
respectively. 

5 An alternative solution that reduces the number of despreading operations utilizes 

the following cumulative multi-codes for / = l,...,N m : 



NI 



E c*' u - 



(139) 



10 Despreading with these cumulative codes yields: 



n — n 



NI 


+ 


NI 

yw'.u 

— PCM,n 

i = l 


sjNI 







= ITS 



—PCM,n 



(140) 



Averaging the user-codes over interferers does not reduce noise further after 
15 despreading. However, the composite signal § v ^ xu collects an average power from the 
NI interferers and therefore benefits from higher diversity. The cumulative multi-codes 



can be used to implement channel identification as follows: 



H v ' =H V ' + 

n + 1 — n 



(141) 



20 where: 



NI 

E 



(142) 



25 



N m NI 

EE 



= (1 -a)(^'-f) + al±±l 



NJJI 



(143) 
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This downlink version of MC-CDFI, referred to as DSMC-CDFI, is illustrated in Figure 
43. Again, averaging over a smaller set of user-codes reduces the number of 
despreading operations. Use of the 5-CDFI version described in Equations (111) to 
(113) instead of, or combination with, the above scheme 18 , are other alternatives that 
5 reduce the amount of complexity due to despreading operations. 

In situations where a pilot code is transmitted, it can be incorporated into the 
cumulative code. This version which we denote St-MC-CDFI, uses a data modulated 
cumulative code over multi-codes, interferers, and pilot(s) for despreading. We also 
introduce normalization by powers in order to meet situations with significant differences 
10 in powers (such as on the downlink). 57T-MC-CDFI uses the following code for 
despreading: 

! N m NI 

(144) 



1 w - 

c * = f =r£ 



AT 



NI 



/=1 



NI 

S n C k + X*T k 

i-l 



i = l 

where 7r, is the spread pilot signal and \* is a weight factor which may serve to signify 
15 the pilot since no decision errors are associated with it. In the special case 

where \ x = 0, the scheme amounts to 5-MC-CDFI with normalization by powers. 

Moreover, these downlink versions could be combined with the pilot assisted ISR 

embodiments described with reference to Figures 44 and 46. It will be appreciated that 

the pilot channel will not be user-specific but rather specific to the serving base station 
20 or to a group of user stations served by that serving base station. In such cases, the pilot 

power is relatively strong so the corresponding beamformer (see Figure 46) may 

implement simple MRC instead of ISR. 



Summing user-codes over resolution periods does not increase diversity. 
However, use of the 5-CDFI version further reduces noise after despreading. 
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The embodiments of the invention described so far use one transmit (Tx) antenna. 
Adding spatial dimension using multiple Tx antennas provides a means of supporting 
more users. In the following section we present a transmitter structure especially suited 
for high-rate low-processing gain downlink transmission, which can potentially provide 
5 high capacity when ISR is employed at the receiver. 

Multiple-Input. Multiple Output (MIMO) ISR with Space-Time Coding (STC) 

On the uplink, increasing the number of receive (Rx) antennas from one to two 
almost doubles capacity of the system when ISR is employed at the BS. The 

10 improvement results from the additional spatial dimension which allows users to be 
distinguished not only by their code but also by their spatial signature. When, on the 
downlink, a single Tx antenna is employed, all signals originating from one specific base 
station (BS) antenna have the same spatial signature at the antenna array of the receiver. 
It is therefore a demand that the BS transmitter is equipped with multiple antennas and 

15 that a wise space-time coding strategy for transmitting signals is employed. 

Such a MIMO transmitter using STC is shown in Figure 49 for serving base 
station v. This figure shows base station v serving the mobiles indexed from 1 to N u 
located within its cell coverage. It transmits to each of these mobile a corresponding 
stream of data (b„), each power-controlled with the total amplitude using a 

20 multiplier It should be noted that is not the real amplification factor, but 

a product obtained at the receiver and put at the transmitter to ease both data modelling 
and algorithmic description (please refer to description of Figure 4). 

The power-controlled data streams are then fed to a group-selector 110 v to divide 
the N u served mobiled in N G groups (see discussion below). Each group cannot be large 

25 than L mobiles in number. Depending on the cell load N u , a given number n <. L of 
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mobiles will be allocated within each group. For convenience of notation, N IG is fixed 
in Figure 49 to L, i.e., the maximum. 

For each group, identical operation will be performed; hence only operation for 
group # 1 will be described. The n = L data streams in group # 1 will be spread 
5 by codes c u (r), c hL (t) using multipliers 15X/2 1 - 1 , l5Xj2 hL , respectively. The 
spread signals are then summed using adder 16X/1 V>1 . The resulting signal is summed 
to a pilot signal ir x (t) specific to group # 1 using adder 16X/2 V, \ then spread with 
code c v (i) specific to base station v using multiplier 15X/3 V>1 . The resulting spread 
signal is denoted as G x (r). Identical operations on other groups will result in similar 

10 signals. The N G signals resulting from such operations in the N G groups of mobiles are 
denoted as Gl (f), G N Jf). 

These signals are fed to a spatial mapper 120 v . This mapper operates like a 
N G x M T coupler to linearly transform the N G input streams inot M T output streams 
equal in number to the Tx antennas and denoted as a (t) A (t) The first 

15 stream A x (t) is delayed by delay 45X 1 and filtered by the shaping pulse 13-y- 1 (which 
includes carrier frequency modulation), then transmitted over antenna 14X vl as 
signal S (*). Identical operations performed on the mapper output streams result in 
signals S x (t) 9 S M (t) being separately transmitted on antennas 14X V, \ l4X VyMr , 
respectively. The structure of this transmitter is particularly interesting in high data-rate, 

20 low-processing gain situations. First, the data sequences b VtU -(t) of the N u users, are 
scaled with V"(r), where <p v < u (t) 2 , u = 1,- 9 N is the desired transmit power. These 
signals are then grouped in js/ c groups. A user belonging to group g is assigned a user 
specific short channelization code drawn from a fixed set of N JG L-chip codes (L is the 
processing gain) which we conveniently organize column-wise in the 

25 matrix c = [c 8 > l (t) T , - 9 c 8tN ' G (t) T ] T and denote it the code-set of group g. Codes 
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belonging to same code-set are all chosen mutually orthogonal such that c H C is 

& S 

diagonal, which in its turn means that the number of fixed channelization codes per 
group is limited to n ig < L\ at the same time code-sets across groups should have low 
cross-correlation, as will be explained later. It should be noted that code-sets are 
5 assumed to be reused in all sectors. When user signals of the same group are coded by 
their respective channelization code, they are summed to provide a one streamed signal. 
A group-specific pilot code is added, and the resulting signal is scrambled by a BS 
specific PN code to arrive at the total group signal G v g {i). Using a linear mapping 
function M, the N G group signals are mapped onto M T antenna branches, to arrive at 
10 A\f) = MG\f) 9 (145) 

where A%t) = [A?(t) T , -^O 7 ], and = [^(O^-G^Wn. Branch signals are 

c 

finally delayed to allow for transmit delay diversity, then formed by the chip-pulse 
matched filter. 

The physical channel matrix which defines transmission from m t Tx antennas 
15 of the BS indexed v to the M receive antennas of the mobile with couple 
index is: 



(146) 



20 

where the (ij)-th element of the matrix is the channel between the y-th Tx antenna of 
the base-station and the /-th Rx antenna of the mobile. This definition allows the signal 
transmitted from the base-station U) when received at the antenna array of the desired 
mobile, to be written as: 
25 X v - W (f) = £T W (0 ® S\f) (147) 
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= £r w (0 ® MG\f) ® £>(*) (148) 

= IT M (t)M ® [G v (r) ® (149) 

= D(t) ® W(t)M ® G v (0 (150) 



where £>(f) = - </>(f -A M )] r represent the transmit delays and chip-pulse 

5 forming. From Equation (149) it is clear that the effective spatial mapping function due 
to the channel becomes H(t)M, which is generally not orthogonal which in its turn means 
that groups will interfere with each other. Therefore, unless the transmitter has 
knowledge of the channel, orthogonality at reception cannot be promised. Design of the 
mapping function will be described later. Equation (150) further shows that the delay 
10 elements can be regarded as part of the channel, and it is clear that they provide virtual 
multi-path. This is normally referred to as "delay transmit diversity." 

If the channel as seen by the vector of group signals G\t) is denoted 
by r v '"(0 = D(t) ® H v > u (t)M, then the signal received at the antenna array may be 
written as: 

!5 X v *(t) = ^T vu {t) ® G\t) (151) 

v=l 

= £ IT ® g/(0 + ... + r- <8> g; c (0 (152) 

v=l 

where Y v / U is the j th column of r v u . From Equation (20) it is noted that each group 
20 of users' signals propagate through the same channel, and that the receiver 
sees N B N G sources. When the proposed MIMO downlink transmitter is employed, the 
model for the received matched filtered signal (Equation 3) can be written as: 



25 



Y & = E E yu ' u '(') + w 



(153) 




(154) 



v=l u=l 



= EE E <p v *<fW2*(f) <s» b v "(t)c«'%f)) + aw 



(155) 



^=1 5=1 « gg; 



5 



(156) 



from which it is clear that users belonging to the same group of same BS, have the same 
channel response H£ 8 {f). This allows for improved identification of the channel at the 
mobile receiver. The most apparent difference from the uplink model is hence that users 

10 belonging to the same group experience the same channel, a feature well exploited on 
the downlink using CDFI (see previous descriptions). 

Appropriate choice of codes across groups is important since the resultant channel 
is not generally orthogonal. Since channelization codes are chosen from a fixed set, sets 
with good properties can be found by optimization. It should be noted that, since the 

15 same scrambling code is used across groups, cross-correlation properties, once set by 
proper choice of channelization code-sets, are preserved after scrambling and hence after 
transmission. Here, only the situation with two groups (N = 2) is considered since 
it is particularly simple. c x is first chosen as an orthogonal matrix, for instance the 
Hadamard matrix (or part of it). It is noted that, if A is a diagonal matrix having ± 1 

20 entries, then AC l will still be an orthogonal span because A is unitary. Therefore, the 
second set of codes is defined by c = AC , where A satisfies 




(157) 



where 



25 



A = C/VC, 



(158) 
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which, in most cases, is easily solved by search. When sets are 
full (i.e., N IG = L), this amounts to a 45° rotation between sets. In the more general 
case where N G > 2, the optimization becomes difficult especially when the processing 
gain is large. Here a blend of guessing and search must be invoked. 
5 Experience shows that reusing the same or opposite code in different groups is 

not attractive. The number of users (total number of possible channelization codes) is 
therefore limited by 2 L ~ l . This limit clearly suggests that a low processing gain also 
encompasses a very limited potential capacity. For instance if l = 2, only two users 
can be supported, and it can be verified that there is no profit from multiple Tx antennas, 

10 since one Tx antenna can already provide near orthogonal transmission (using orthogonal 
codes). With L = 4, the limit is N u < 8, and for this situation, there is benefit from 
going from one to two Tx antennas, but not three. In general, the number of Tx 
antennas for a cost-effective system is limited by m < ^ * because adding more 
antennas will not provide potential increase of capacity. 

15 The purpose of the spatial mapping function is to assign a unique M T -dimensional 

spatial signature to each of the w groups. In situations where the transmitter has some 
knowledge of the channel 19 , the mapping function may be chosen as a function of time 
such that the resultant channel r v,w (r) becomes orthogonal. In mobile cellular CDMA, 
however, this is not an option because the physical channels are user specific, and 

20 therefore an orthogonality condition can only be provided for one receiving user. 

Considering therefore the design of a fixed mapping function, the rank properties 
of the resultant channel r v "(0 = D(f) ® H v > u (i)M are intuitively optimized if the 
mapping function M is chosen to have the best rank properties. In the case where the 



19 E.g., feed-back from the receivers or time-division duplex. 
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number of groups (pj^ equals the number of Tx antennas (M T ), optimal rank 
properties are obtained by choosing an orthogonal mapping function. Identity 
mapping (M = I) is orthogonal, and simply maps group 1 to antenna 1, group 2 to 
antenna 2, etc. However, it causes unequal load on antennas, and delay transmit 
5 diversity is not exploited. These problems are both avoided if the Hadamard matrix is 
used, since it distributes signals equally on antenna branches and therefore exploits 
"delay transmit diversity" and avoids power imbalance. 

In the more general case, power control can be distributed among the transmit 
antennas. Such power control distribution techniques are well known and will not be 
10 described in detail here. 

It should be noted that embodiments of the invention are not limited to using the 
space-time coding scheme described with reference to Figure 49 but could use other 
known space-time coding schemes. 

It will be appreciated that the base station transmitter described with reference to 
15 Figure 49 does not require a modification to the receiver at the user station, i.e., any of 
the receivers described with reference to Figures 40-43 may be used to receive its 
signals. Actually, with a MIMO system, the receiver will "see" the MIMO base station 
transmitter as N G sub-base stations corresponding to the group partitioning of Figure 49. 

It should be noted that each user station could have a plurality of transmit 
20 antennas and use a MIMO transmitter similar to that of the base station and described 
with reference to Figure 49. Of course, the corresponding receiver at the base station 
will not require modification for the reasons given above. 
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TABLE 5 





BR 


c v BS , n - 


|_«=i i=\ f=\ k=-i 




NI N m Nj +1 




i 



10 TABLE 5 

Table 5 shows base-specific constraint and blocking matrices c v and c v ' which 

BS,n BS,n 

will apply to the modes shown in Table 3 except for TR, replaced by BR. Indices of 
remaining modes in Table 3 should be modified to include the index of the base-station 
v as shown for the TR mode. It should be noted that channel coefficients t u hold the 
15 index of the base-station u instead of the user /. Transmissions to all user-stations from 
base-station u propagate to the desired user station through a common channel. It should 
also be noted that summation over users is weighted by the estimate of the total 
amplitude due to user- independent power control. Definitions of 0 and \[ are 
given in Table 3. 

20 



TABLE 6 





TR 




NB NI N m N r + i 

E E C"E E E 

|_v=l C/=l 1=1 f=l k=-\ 




NB NI N m N f + i 

E E *«"E E E 

v=l u=l 1 = 1 f=\ k=-\ 



118 





l 




BR 




M N m N r + i 
w=l /=1 / =1 k--\ 




NI N m N f + i 1 

E ttT E E &;XXuf^ ••• 

|_ M=l /=1 /=1 *=-l 




A© 



TABLE 6 

10 

Table 6 shows multi-base constraint and blocking matrices C n and C?' 4 ' 4 ' which apply 
to the modes of Table 5 by row-wise aligning the constraint and blocking 
matrices and C"'/ '' from base-stations into larger matrices r and fT"' ''*'' in 

Bo,n BS,n n n 

the way suggested by Equations (104) and (105). The number of constraints in Table 
15 5 is multiplied by NB as shown here for the BR mode. The additional TR mode sums 

the constraint- vectors of the BR mode over all base-stations. The definition of \\ is 

given in Table 3 and 1 V = o if (v,i,l,k) = (v', V, V,k') and 1 otherwise. 

It should be appreciated that, whether ISR is used for the uplink or the downlink, 

it will function either as a single antenna or multiple antenna for reception or 
20 transmission. 

Embodiments of the invention are not limited to DBPSK but could provide for 
practical implementation of ISR in mixed-rate traffic with MPSK or MQAM modulations 
without increased computing complexity. Even orthogonal Walsh signalling can be 
implemented at the cost of a computational increase corresponding to the number of 
25 Walsh sequences. Moreover, different users could use different modulations. Also, one 
or more users could use adaptive coding and modulation (ACM). 

It is also envisaged that embodiments of the invention could employ carrier 
frequency offset recovery (CFOR). It should be appreciated that the decision rule units 
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do not have to provide a binary output; they could output the symbol and some other 
signal state. 

It should also be noted that, although the above-described embodiments are 
asynchronous, a skilled person would be able to apply the invention to synchronous 
5 systems without undue experimentation. 

The invention comprehends various other modifications to the above-described 
embodiments. For example, long PN codes could be used, as could large delay-spreads 
and large inter-user delay-spreads. 

For simplicity, the foregoing description of preferred embodiments assumed the 
10 use of short spreading codes. In most practical systems, however, long spreading codes 
would be used. Because the portion of the long code differs from one symbol to the 
next, certain operations, which were unnecessary for short codes, will have to be 
performed, as would be appreciated by one skilled in this art. For further information, 
the reader is directed to references [22] and [23]. It is envisaged, however, that short 
15 codes could still be used during acquisition and the long codes used once a link has been 
established. 

During the acquisition step, a user station could be required to connect using one 
of a plurality of predetermined (short) codes. The null-constraints used by the receiver 
then would be preselected to cancel signals using such predetermined codes. This would 
20 avoid problems arising when a user station begins to transmit and for which the receiver 
has not derived any constraints. Such a modification would be applicable to the 
downlink situation and use ISR-H. 
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CLAIMS: 



1. A receiver suitable for a base station of a CDMA communications system 
comprising at least one base station (11) having a transmitter and a said receiver and a 
5 multiplicity (U) of user stations (10\...,10 u ) including a plurality (U') of user stations 
served by said at least one base station, each user station having a transmitter and a 
receiver for communicating with said at least one base station via a corresponding one 
of a plurality of channels (14 1 ,.. .,14 u ), the base station receiver for receiving a signal 
(X(t)) comprising components corresponding to spread signals transmitted by the 
10 transmitters of the plurality of user stations, each of said spread signals comprising a 
series of symbols spread using a spreading code unique to the corresponding user station, 
said base station receiver comprising: 

a plurality (U') of receiver modules (20 1 ,...,20 NI ,20 d ) each for deriving from 
successive frames of the received signal (X(t)) estimates of said series of symbols 
15 of a corresponding one of the user stations, 

preprocessing means (18) for deriving from the received signal (X(t)) a series of 
observation matrices (Y n ) each for use by each of the receiver modules (20) in a 
said frame to derive an estimate of a symbol of a respective one of said series of 
symbols, and 

20 means (19,44;44/1 ,44/2) for deriving from each observation matrix a plurality of 

observation vectors (Y\Y ; Z l ...Z m : Z**) and apply ing each of the observation 

— n — n-l — n — n ~n 

vectors to a respective one of the plurality of receiver modules (20 1 ,...,20 NI ,20 d ); 
each receiver module comprising; 

channel identification means (28) for deriving from one of the observation 
25 vectors a channel vector estimate (H l fj 1 ^ 1 ^ Y* ) based upon 

parameter estimates of the channel between the base station receiver and 

the corresponding user station transmitter; 

beamformer means (27 \ . . . ,27 NI ,27 d ; 47 d ) having coefficient tuning means 
(50) for producing a set of weighting coefficients in dependence upon the 
30 channel vector estimate, and combining means (51,52) for using the 

weighting coefficients to weight respective ones of the elements of a 
respective one of the observation vectors and combining the weighted 
elements to provide a signal component estimate an d 
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symbol estimating means (29\...,29 u , 30\...,30 u ) for deriving from the 
signal component estimate an estimate ($ l 9 ... 9 b u ) of a 
symbol (b*, ...,b?) transmitted by a corresponding one of the user 
stations (10 1 ..,10 u ), 

5 wherein said receiver further comprises means (42,43) responsive to symbol 

estimates (b*,...,b™;g\g 2 ,g 3 ;g l + l > n ) and to channel estimates 
(M l R --.Mn ! >M i n _d comprising at least said channel vector 
estimates (H\.^H NI ) for channels (14 1 ,...,14 NI ) of a first group (I) of said 

n n 

plurality of user stations (10 1 ,...,10 NI ) to provide at least one constraint 
10 matrix (c ) representing interference subspace of components of the received 

signal corresponding to said predetermined group, and in each of one or more 
receiver modules (20A d ) of a second group (D) of said plurality of receiver 
modules, the coefficient tuning means (50 A d ) produces said set of weighting 
coefficients in dependence upon both the constraint matrix (c ) and the channel 
15 vector estimates so as to tune said one or more receiver modules (20 A d ) 

n 

each towards a substantially null response to that portion of the received signal 
(X(t)) corresponding to said interference subspace. 

2. A receiver according to claim 1, wherein said coefficient tuning means (50A d ) 
20 also tune said one or more receiver modules (20 A d ), respectively, each towards a 

substantially unity response for the component of the received signal (X(t)) from the 
transmitter of the corresponding one of the user stations. 

3 . A receiver according to claim 1 , wherein the observation vector deriving means 
25 (19;44) comprises first reshaping means (44) for reshaping the observation matrix (YJ 

from the preprocessing means (18) and supplying the resulting observation vector (YJ 
to said beamformer means (47 A d ) of each of said one or more receiver modules (20 A d ), 
and wherein the means (42,43) for providing said at least one constraint matrix 
comprises constraints-set generating means (42A) responsive to said channel 
30 estimates (^...X;^,) symbol estimates (b l n ,..J™;g\g\g 3 ;g 1 ^) 
corresponding to said first group for generating a plurality of constraints-set 
matrices (cj,,..., £? n c ) together characterizing the subspace of interference attributable to 
said first group of user stations and the constraint matrix generating means (43A) 
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comprises a bank of vector reshapers (4&4 48.4 Nc ) for reshaping the constraints-set 
matrices (C^,..., df n c ) to form columns, respectively, of the constraint matrix (C ), the 
constraint matrix generating means (43 A) supplying the constraint matrix to each of said 
coefficient tuning means (50A d ) of said one or more receiver modules (20 d ) of said 
5 second group, and wherein, in each of said one or more receiver modules (20 A d ), the 
channel estimation means (28A d ) supplies spread channel vector estimates (y 4 ^ to the 

v — 0,n J 

coefficient tuning means (50A d ) for use in updating said weighting coefficients. 

4. A receiver according to claim 3, wherein said constraint matrix generating means 
10 (43A) comprises transformation means (49A) for forming an inverse matrix (QJ in 

dependence upon the constraint matrix (C n ) and supplying said inverse matrix to said 
coefficient tuning means (50 A d ) of said one or more receiver modules (20 A d ), and 
wherein said coefficient tuning means (50A d ) computes said weighting coefficients in 
dependence upon said constraint matrix, said inverse matrix and said channel vector 
15 estimate. 

5 . A receiver according to claim 1 , wherein the observation vector deriving means 
comprises a plurality of despreaders (19 1 , . . . , 19 NI , 19 d ) each for using a corresponding one 
of the user spreading codes to despread the observation matrix (TJ using a respective 

20 one of the spreading codes to form a user-specific post-correlation 
vector (Z\...,Z N * Z**) and supplying same to a respective one of the channel 

— n — n — n 

identification means (28\. ..,28 NI , 28 d ). 

6. A receiver according to claim 3, further comprising means (30 1 ,...,30 NI ) for 
25 deriving amplitude estimates (^,...,#* 7 ) of signal components from said first group of 

user stations and supplying the amplitude estimates to said constraints-set generating 
means as parts of said channel estimates, and wherein the constraints-set generating 
means (42C) comprises a plurality of respreaders (57C 1 ,. ..,57C NI ) each for using a 
corresponding one of the user spreading codes to respread a respective one of the symbol 
30 estimates (^ /j 1 ,...,foj v/ ) from the receiver modules corresponding to said first group (I) 
of user stations, scaling means (58C 1 , . . . ,58C NI ) for scaling the respread symbol estimates 
by the amplitudes (^,...,^ / ) of the signal components corresponding to the symbol 
estimates 
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(b l n , B™) 9 respectively, and a plurality of channel replication means (59C 1 , . . . ,59C NI ) 
having coefficients adjustable in dependence upon the channel vector 
estimates (H 1 ,...,//* 7 ), respectively, for filtering the corresponding respread and scaled 

n n 

symbol estimates to provide user-specific observation matrix 
5 estimates .(f^_ p f^), respectively, and means (60) for summing the user-specific 
observation matrices to form an observation matrix estimate (/^ ) s and supplying same 
to the constraint matrix generator means (43C), the constraint matrix generator means 
(43 C) comprising vector reshaping means for reshaping the observation matrix 
estimate (I n l ) to form an observation vector estimate (/ ) as a single column 
10 constraint matrix (C ) for application to the coefficient tuning means (50A d ) of each of 
said receiver modules (20 d ) of said second group (D). 

7. A receiver according to claim 3, wherein the constraints-set generating means 
(42D) generates a number of constraints (N c ) equal to the number (NI) of user stations 

15 in said first group (I) and comprises a plurality of respreaders (57D\...3,57D NI ) each for 
using a corresponding one of the user spreading codes to respread a respective one of 
the symbol estimates b™) from a predetermined group (I) of said receiver 

modules corresponding to said selected ones of said components of the received signals, 
and a plurality of channel replication means (59D 1 ,...,59D NI ) having coefficients 

20 adjustable in dependence upon the channel vector estimates (^ 1 ,...,fl A7 ) ) respectively, 

n n 

for filtering the corresponding respread symbol estimates to provide a plurality of user- 
specific observation matrix estimates (f^ V -,V^ 1 ) respectively, and wherein, in the 
constraint matrix generator (43D), the bank of vector reshapers (4SA l 9 .,48/1 Nc ) reshape 
the user-specific observation matrix estimates to form a plurality of user-specific 
25 observation vector estimates (y 1 respectively, as respective columns of the 

constraint matrix (C n ) for supply to each of the coefficient tuning means (50A d ) of each 
of said receiver modules (20 d ) of said second group (D). 

8. A receiver according to claim 3, wherein the plurality of channel identification 
30 means (28 1 , . . . ,28 NI ) of said first group of receiver modules (20E 1 , . . . ,20E NI ) provide both 

the channel vector estimates (/y 1 ,... } 7y M ), respectively, and a plurality of sets of sub- 

n n 

channel vector estimates (fl]; x ,...,H^\...,H N J' 1 ,...,H N n ^ f \ respectively, each of the 
channel estimates (J^ n , J^ 7 ) comprising a respective one of the sets of sub-channel 
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vector estimates, each of the sets of sub-channel vector estimates representing an estimate 
of the channel parameters of N f sub-channels of said channel between the base station 
and the transmitter of the corresponding one of the NI user stations in said first group, 
the constraints-set generator means (42E) comprises a plurality of respreaders 
5 (57E 1 ,...,57E NI ) and a plurality of channel replicators (59E 1 ,.. .,59E NI ) coupled to the 
plurality of respreaders (57E\...,57E NI ), respectively, for filtering the plurality of 
respread symbols _ 9 &™), respectively, using respective ones of the sub-channel 
vector estimates (^'^...^^..^i/^'^...,^ 7 ^), to form a plurality (NJ of constraints 
equal in number to NjNI and forming a plurality of user-specific sub-channel observation 

10 matrix estimates (F^, f^ 9 F^,...,^fo corresponding to the sub-channels, 
respectively, and, in the constraint matrix generator (43E), the bank of reshapers 
(4SA l 9 ... 9 4SA Ne ) reshape the sets of user-specific observation matrix estimates to form 
a corresponding plurality of sets of user-specific sub-channel observation vector 
estimates (f 1 - 1 y 1 y* 7 - 1 . j^"/) as respective columns of the constraint 

15 matrix (C n ) for supply to each of the coefficient tuning means (50 A d ) of each of said 
receiver modules (20 d ) of said second group (D). 

9. A receiver according to claim 3, further comprising hypothetical symbol estimate 
generating means (63F 1 ,.. .,63F NI ) for generating for each of said selected ones of said 
20 components a series of hypothetical symbol estimates (g 1 ,^ 2 ,^ 3 ), wherein the 

n n n 

constraints-set generating means (42F) comprises a plurality of respreaders 
(57F\...,57F NI ) each for respreading, using a corresponding one of the user spreading 
codes, selected sets of said hypothetical symbol estimates, and a plurality of channel 
replicator means (59F 1 ,...,59F NI ), respectively, for filtering the sets of respread symbol 
25 estimates, each to form a plurality (NJ of constraints equal in number to 3 NI and 
forming a plurality of user-specific observation matrix estimates 
(l^,f^,f^^ the plurality of channel replication means 

(59F 59F NI ) having coefficients adjustable in dependence upon the channel vector 
estimates (H l , ...,fl^ 7 ), respectively, and, in the constraint matrix generating means 

n n 

30 (43F), the bank of reshapers (4&4 4&4 * c ) reshape the sets of user-specific 
observation matrix estimates to form a plurality of user-specific observation vector 
estimates (f f f* 1 f 41 y* 7 \ respectively, as respective columns of the 



129 

constraint matrix (c ) for supply to each of the coefficient tuning means (50A d ) of each 
of said receiver modules (20 d ) of said second group (D). 

10. A receiver according to claim 3, further comprising means (63G\...,63G NI ) for 
5 providing hypothetical symbol estimates (g /+1 - n ) and wherein the constraints-set 
generating means uses a combination of said symbol estimates (h 1 h NI } and said 
hypothetical symbol values in producing said set of user-specific observation matrix 
estimates. 

10 11. A receiver according to claim 10, wherein the constraints-set generating means 
(42G) comprises a plurality of respreaders (57G 1 ,...,57G NI ) each for respreading, using 
a corresponding one of the user spreading codes, a respective one of the symbol 
estimates (b*,... 9 b„ r ) and a said hypothetical symbol estimates to provide a 

plurality of respread symbol estimates, and a plurality of channel replicator means 

15 (59G\...,59G NI ), respectively, for filtering the respread symbol estimates to form a 
plurality (NJ of constraints equal in number to 2NI and forming a plurality of user- 
specific observation matrix estimates ( f l Y l Y* 7 Y* 7 ^ the plurality of channel 
replication means (59F 1 , . . . ,59F NI ) having coefficients adjustable in dependence upon the 
channel vector estimates (H 1 ,...,!^ 1 ), respectively, and, in the constraint matrix 

n n 

20 generating means (43G), the bank of vector reshapers (4&4 1 9 ... 9 4%A Nc ) reshape the 
user-specific observation matrix estimates to form a corresponding plurality of user- 
specific observation vector estimates (V f 1 ... I^ 7 Y* 1 ) respectively, as respective 

— r,n — + '"' — — +l,n 

columns of the constraint matrix (C n ) for supply to each of the coefficient tuning means 
(50A d ) of each of said receiver modules (20 d ) of said second group (D). 

25 

12. A receiver according to claim 1, wherein the observation vector deriving means 
comprises a plurality of despreaders (19 1 ,...,19 NI ,19 d ) each for despreading the 
observation matrix (Y n ) using a corresponding one of the user spreading codes to 
produce a corresponding one of a plurality of post-correlation observation 
30 vectors (Z\...,Z NI , Z*) and supplying the post-correlation observation vectors to both 

n n — n 

the channel identification means and the coefficient tuning means of the beamformer 
means of each of said one or more receiver modules, said elements of the observation 
vector weighted by said combining means of said one or more receiver modules being 
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elements of the corresponding post-correlation observation vector, and wherein the 
constraint matrix providing means (42B, 43B) comprises means (42B) responsive to said 
channel estimates (j^,..., J^ 7 ; and to symbol estimates 

(b*,...,b™\g\g 1 ,g 3 ;g l+hn ) corresponding to said first group (I) of user stations for 
5 providing a plurality of constraints-sets matrices (CJ together characterizing the subspace 
of interference attributable to said spread signals of said first group of user stations, and 
constraint matrix generating means (43B) comprising a plurality of user-specific 
constraint matrix generators (43B d ) each associated with a respective one of said one or 
more receiver modules (20 d ), each user-specific matrix generator (43B d ) having 
10 despreading means (55B d >\ ...,55B d ' N <) for using the corresponding user spreading code 
of the specific user to despread each of the user-specific constraints-set matrices 
respectively, to form a respective column of a corresponding one of a plurality of user- 
specific post-correlation constraint matrices (C d pcMr ), the plurality of user-specific 
constraint matrix generating means (43B d ) supplying said plurality of user-specific post- 
15 correlation constraint matrices to the coefficient tuning means of the respective one of 
said one or more receiver modules (20B d ). 

13. A receiver according to claim 12, further comprising means (SO 1 ,...^^ 1 ) for 
deriving amplitude estimates #* 7 ) of signal component estimates of said first 

20 group of user stations and supplying the amplitude estimates to said constraints-set 
generating means as parts of said channel estimates, and wherein the constraints-set 
generating means (42C) comprises a plurality of respreaders (57C\...,57C NI ) each for 
using a corresponding one of the user spreading codes to respread a respective one of 
the symbol estimates from receiver modules (20 1 ,...,20 NI ) corresponding to said first 

25 group of user stations, scaling means (58C 1 ,...,58C NI ) for scaling each of the respread 
symbol estimates by said amplitudes respectively, a plurality of channel 

replication means (59C 1 , . . . ,59C NI ) having coefficients adjustable in dependence upon the 
channel vector estimates (h\ ..^ff 11 ), respectively, for filtering the corresponding 

n n 

respread and scaled symbol estimates to provide user-specific observation matrix 
30 estimates (F^ p f^), respectively, and means (60) for summing the user-specific 
observation matrix estimates to form an observation matrix estimate (/ ) and supplying 
same to each of the user-specific constraint matrix generators (43H d ) of said one or more 
receiver modules of said second group, each of the user-specific constraint matrix 
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generators (43H d ) comprising despreading means (55B d ) for despreading the observation 
matrix estimate (/ ) using the corresponding user spreading code to form a respective 
one q_ d ) of a plurality of post-correlation user-specific observation vector estimates 

PCM,n - 1 

each as a single column respective constraint matrix (C d ) for use by the associated 
5 one of the coefficient tuning means (50B d ) in said second group (D). 

14. A receiver according to claim 12, wherein the constraints-set generating means 
(42D) comprises a plurality of respreaders (57D 1 ,.. .,57D NI ) each for respreading a 
respective one of the symbol estimates (b*,...,b™) from a predetermined group (I) of 

10 said receiver modules corresponding to said selected ones of said components of the 
received signals, and a plurality of channel replication means (59D\...,59D NI ) having 
coefficients adjustable in dependence upon the channel vector estimates (H l , 

n n 

respectively, for filtering the corresponding respread symbol estimates to provide a 
plurality of user-specific observation matrix estimates (f l Y* 7 ) and wherein, in 

15 each of the user-specific constraint matrix generating means (43I d ), the despreading 
means (55B d *\...,55B dyNc ) despreads the user-specific observation matrix estimates to 
form a plurality of user-specific observation vector estimates (y* ,.„,$ A7 ) 
respectively, as respective columns of a respective user-specific constraint 
matrix (C PCMn ) for supply to the associated one of the coefficient tuning means (50B d ) 

20 of said one or more receiver modules (20 d ) in said second group (D). 

15. A receiver according to claim 12, wherein the plurality of channel identification 
means (28 1 , . . . ,28 NI ) of said first group (I) of receiver modules (20E 1 , . . . ,20E NI ) generate 
both said channel vector estimates (H 1 , ..., J H A7 ), respectively, and a plurality of sets of 

n n 

25 sub-channel vector estimates (^^^...^^..^T/^^.,.^^ 7 ^), respectively, each of the 
channel vector estimates (//*,.. comprising a respective one of the sets of sub- 

n n 

channel vector estimates, each of the sets of sub-channel vector estimates representing 
an estimate of the channel parameters of sub-channels of said channel between the base 
station and the transmitter of the corresponding user station, the constraints-set generator 
30 means (42E) comprises a plurality of respreaders (57E\...,57E NI ) and a plurality of 
channel replicators (59E 1 ,...,59E NI ) coupled to the plurality of respreaders 
(57E\...,57E NI ), respectively, for filtering the plurality of respread symbols 
(£ n \..., b™) 9 respectively, of said first group using respective ones of the sub-channel 
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vector estimates (T/^V^tfJ/V to form a plurality of sets of user- 

specific observation matrix estimates (f^\... ,Y l n ' N {, , f"^ 1 , tf^f'), ^ sets 
corresponding to the sub-channels, respectively, and, in each of the user-specific 
constraint matrix generators (43K d ), the despreading 
5 means (55B d > 1 ,... ,55B d,Nc ) despreads the sets of user-specific observation matrix 
estimates to form a corresponding plurality of sets of user-specific observation vector 
estimates (j u jdM,N f . the sets f orm i n g respective columns of a respective 

user-specific constraint matrix (C d pcM n ) for supply to the associated one of the coefficient 
tuning means (50B d ) of the one or more receiver modules (20 d ) in said second group (D). 

10 

16. A receiver according to claim 12, further comprising hypothetical symbol estimate 
generating means (63L 1 ,...,63L NI ) for generating for each of said selected ones of said 
components a series of hypothetical symbol estimates (g\g 2 ,g 3 ), wherein the 

n n n 

constraints-set generating means (42L) comprises respreading means (57L 1 , . . . ,57L NI ) for 
15 respreading, using a corresponding one of the user spreading codes for said first group, 
selected sets of said hypothetical symbol estimates, and channel replicator means 
(59L 1 ,...,59L NI ), respectively, for filtering the sets of respread hypothetical values, each 
to form one of a plurality of sets of user-specific observation matrix 
estimates ^ f[ { n , f\ x n , w , n ), the plurality of channel replication means 
20 (SQL 1 ,...^?)!/") having coefficients adjustable in dependence upon the channel vector 
estimates (H l \ respectively, and, in each of the user-specific constraint matrix 

n n 

generating means (43L d ), the despreading means (55B dA ,...,55B d,Nc ) despreads the sets 
of user-specific observation matrix estimates to form a corresponding plurality of sets of 
user-specific observation vector estimates f d - U2 t* 1 - 3 f dMA f d > Nia f d > NI > 3 \ 

y -PCM,n ' -PCM,n ' ~PCM,n 9 " ' ' ~PCM,n ' ~PCM,n ' ~PCM,n * 9 

25 respectively, as respective columns of a respective user- specific constraint 
matrix (C d pcMt ) for supply to the coefficient tuning means (50B d ) of the associated one 
of said one or more receiver modules (20 d ) in said second group (D). 

17. A receiver according to claim 12, further comprising means (63M 1 , . . . ,63M NI ) for 
30 providing hypothetical symbol estimates (V^) and wherein the constraints-set 

generating means (42M) uses a combination of said symbol estimates (& 1 ,...,b NI ) and 
said hypothetical symbol estimates in producing said sets of user-specific observation 
matrix estimates. 
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18. A receiver according to claim 17, wherein the constraints-set generating means 
(42M) comprises a respreading means (57M 1 ,...,57M NI ) for respreading, using the 
symbol estimates and said hypothetical symbol estimates (g^) using 

respective ones of the user spreading codes, to provide a plurality of respread symbol 
5 estimates, and channel replicator means (59M 1 ,...,59M NI ) for filtering the respread 
symbol estimates to form a plurality of pairs of user-specific observation matrix 
estimates (^J^^.^^fjj, each pair corresponding to one of said first group of 
user stations, the channel replication means (59M 1 ,. ,59M NI ) having coefficients 
adjustable in dependence upon the channel vector estimates (H\ ... J ff /I ) 9 respectively, 

n n 

10 and, in each of the user-specific constraint matrix generating means (43M d ), the 
despreading means (55B d ' l ,.,. 9 55B dtNc ) despreads the user-specific observation matrix 
estimates to form a corresponding plurality of user-specific observation vector estimates 
(I d,l ' ki Z^ 1 ^, respectively, as respective columns of a corresponding 

K -PCM,n -PCMy " ~PCM,n -PCM,n'" r & 

user-specific constraint matrix (C PCMn ) for supply to the coefficient tuning means (50B d ) 
15 of the associated one of said one or more receiver modules (20 d ) in said second group 
(D). 



19. A receiver according to claim 1, wherein: 

at least one (20 1 ) of said one or more receiver modules both supplies symbol 
20 estimates to said constraint-matrix generating means (42P,43P) for use in deriving said 
constraint matrices and uses said constraint matrices in tuning weighting coefficients of 
its beamformer means (47P 1 ), 

said at least one of said one or more receiver modules further comprising a 
second beamformer means (27P 1 ) having second combining means and second coefficient 
25 tuning means for tuning the weighting coefficients of the second combining means; 

the observation vector deriving means (19,44;44/l,44/2) further comprising 
second reshaping means (44/2) for reshaping said observation matrix (YJ and supplying 
a resulting second observation vector ( y ) to said second combining means (27P 1 ) and 

n 

delay means (45) for delaying the first observation vector Q^_ 7 j relative to the second 
30 observation vector (Y ); 

the second coefficient tuning means also being arranged to tune the weighting 
coefficients of the second combining means in dependence upon the channel vector 
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estimate (f 4 ) used by the first coefficient tuning means in beamformer (47P0; and 
wherein 

the channel identification means (28P') of said at least one of the one or more 
receiver modules derives the channel vector estimate ( y* } from the delayed 
5 first observation vector (y ) and supplies said channel vector estimate to the 

respective coefficient tuning means of the first combining means and the second 
combining means for use in updating their weighting coefficients and supplies to 
the constraint matrix generator means (42P,43P) said channel 
estimates (J-C^) for use in deriving the constraint matrix; 
10 the first combining means (47P j ;5 1 ,52) and second combining means (27P i ;5 1 ,52) 

use their respective weighting coefficients to weight respective ones of the elements of 
the delayed first observation vector and the second observation vector, respectively, and 
combine the weighted elements of the respective first and second observation vectors to 
provide first signal component estimate (§*_ x ) and second signal component 
15 estimate (s^ RCn ), respectively; 

said at least one of said one or more receiver modules further comprises second 
symbol estimating means (29P/20 for deriving from the second signal component 
estimate (s^ RCn ) a symbol estimate (b^ RCn ) and supplying said symbol 
estimate (£ ' ) to the constraint matrix providing means (42P,43P); 
20 the constraint means (42P,43P) comprises constraints-set generating means (42P) 

for generating a plurality of constraints-set matrices {C n _ 1 ) together characterizing the 
subspace of interference attributable to said first group of user station signals and the 
constraint matrix generating means (43P) comprises vector reshaping 
means (4&4 l ,...,4SA Nc ) for reshaping the constraints-set matrices (C^,..., {f n c ) to form 
25 columns, respectively, of the constraint matrix (C j), the constraint matrix generating 
means (43P) supplying the constraint matrix to each of said coefficient tuning means 
(50P d ) of said one or more receiver modules (20 d ). 



20. A receiver according to claim 19, wherein said at least one of said one or more 
30 receiver modules further comprises amplitude estimation means (30?*) for providing 
amplitude estimates (^_ 7 ) of signal components from said first group of user stations 
and supplying the amplitude estimates to the constraints-set generating means (42P) as 
parts of the channel estimates ([K^)- 
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21. A receiver according to claim 19, wherein said at least one of the one or more 
receiver modules is arranged to perform for each frame a plurality of iterations (#1, #2, 
#N S ) to derive each symbol estimate, the arrangement being such that: 

in each of the iterations of a particular frame (n), the constraints-set generator 
5 (42P) uses the channel estimates (JhC n ^), a first symbol estimate (p l ) and a second 
symbol estimate (b l MRCn ) from the first and second beamformer means, respectively of 
said one or more receiver modules, 

in a first iteration, the constraint matrix providing means (42P,43P) generates a 
first-iteration constraint matrix ((7^(1)) in dependence also upon a previous symbol 
10 estimate {b l MRCn _ x ) previously generated by each of said second beamformers of said one 
or more receiver modules and supplies said first-iteration constraint matrix (£ n j(l)) to 
the coefficient tuning means of the first beamformer means (47P 1 ) for use, with the 

spread channel estimate f l to tune the coefficients of the first beamformer means 

— o,«-i' 

(47P 1 ) for weighting each element of the delayed first observation vector to produce 
15 a first-iteration signal component estimate and the decision rule unit (29P/10 processes 
said first-iteration signal component estimate to produce a first-iteration symbol 
estimate (£^(1)), 

in a second iteration, the constraint matrix generating means (43P) uses said first- 
iteration symbol estimate (fcj^l)), instead of the previous symbol estimate {b i MRC x ) to 

20 tune the weighting coefficients and derive a second-iteration constraint matrix (C (2)) 
for use by the first beamformer means and first decision rule means to produce a second- 
iteration symbol estimate (£^(2)), and 

in a final iteration of a total of N s iterations, the constraint matrix providing 
means (42P,43P) uses a penultimate-iteration symbol estimate ^iV -1)) produced by 

25 the first decision rule means (29P/T) in the penultimate iteration to provide a final- 
iteration constraint matrix (C ^(N)) for use by the first combining means and first 
decision rule means to provide a final iteration symbol estimate ($* r (N )) as the target 
symbol estimate of that frame (u) for output as the symbol estimate (b '_ x ) ? and wherein 
the constraints-set generator (42P) buffers the symbol estimate (b 1 ) for use in every 

30 iteration of the next frame (n + 1) instead of symbol estimate b 1 and the constraints-set 

n—2. 9 

generator (42P) uses a new symbol estimate (b^ RCn+l ) from said second beamformer 
(27P 1 ) for all iterations of the new frame, and uses the previous symbol 
estimate (B^ RCn ) from said second beamformer means (27P') in only the first iteration 
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of said new frame, said previous symbol estimate being buffered as required and other 
variables being incremented appropriately. 

22. A receiver according to claim 1, wherein said plurality of receiver modules 
5 (20 1 , . . . ,20 NI ,20 d ) comprises a first set (I) of receiver modules (20 1 , . . . ,20 NI ) for relatively 
strong user signals and that contribute at least respective sets of said channel estimates 
to said constraints-set generator (42) for use in deriving said constraint matrices but do 
not use said constraint matrices to update the weighting coefficients of their respective 
beamformer means, and a second set (D) of receiver modules (20 d ) for relatively weaker 
10 user signals and that use the constraint matrices to update the weighting coefficients of 
their respective beamformer means but do not contribute either channel estimates or 
symbol estimates to the constraints-set generator (42) for use in deriving said constraint 
matrices. 

15 23. A receiver according to claim 22, wherein said plurality of receiver modules 
further comprises at least one further set (Ml) of at least one receiver module (20*) as 
defined in claim 19 for a user signal having a signal strength that is intermediate the 
relatively strong and relatively weaker user signals and that both supplies at least channel 
estimates to said constraint-matrix generating means (42P,43P) for use by receiver 

20 modules in other sets in deriving said constraint matrices and uses said constraint 
matrices derived from constraints supplied by receiver modules in its own set in tuning 
weighting coefficients of its beamformer means (47P 1 ). 

24. A receiver according to claim 22, wherein said plurality of receiver modules 
25 further comprises at least one further set (M2) of at least one receiver module (20 1 ) as 
defined in claim 19 for a user signal having a signal strength that is intermediate the 
relatively strong and relatively weaker user signals and that both supplies at least channel 
estimates to said constraint-matrix generating means (42P,43P) for use in deriving said 
constraint matrices and, in tuning weighting coefficients of its beamformer means (47P'), 
30 uses said constraint matrices derived from constraints supplied by receiver modules in 
either or both of its own set and other sets. 
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25. A receiver according to claim 23, wherein the receiver modules are arranged in 
hierarchical order according to signal power and each lower power receiver module uses 
a constraint matrix formed from constraints-sets supplied by each of the higher power 
receiver modules. 

5 

26. A receiver according to claim 24, wherein the receiver modules are arranged in 
hierarchical order according to signal power and each lower power receiver module uses 
a constraint matrix formed from constraints-sets supplied by each of the higher power 
receiver modules. 

10 , 

27. A receiver according to claim 23, wherein each of the receiver modules in one 
of said sets is different from the receiver modules in the other sets. 

28. A receiver according to claim 23, wherein the receiver modules in one of said 
15 second and further sets of receiver modules differ from each other. 

29. A receiver according to claim 3, wherein said beamformer means of said at least 
one of the one or more receiver modules comprise projector means (100 d ) for multiplying 
a projection (n^) with the observation vector (Y n ) to form an interference-reduced 

20 observation vector (i* 1 ^), and a residual beamformer (27Q d ) responsive to the 
projection (n^) and to the channel vector estimate ( ) to produce said signal 
component estimate and the channel identification means (28Q d ) derives the 

channel vector estimate (yf ) from the interference-reduced observation vector (Y 11 ^). 

25 30. A receiver according to claim 29, wherein said at least one of the one or more 
receiver modules further comprises reshaping means (102Q d ) for reshaping the 
interference-reduced observation vector (y* 1 ^) to form an interference-reduced 

n 

observation matrix (Y^ ,c/ ), and a despreader (19 d ) for despreading the interference- 
reduced observation matrix ( Y^*), with the corresponding user spreading code to form 
30 a post-correlation reduced-interference observation vector (Z n * d ) for use by the channel 

n 

identification means (28Qd) in deriving said channel vector estimate. 
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31. A receiver according to claim 3, for a base station in a CDMA system in which 
at least one of the user stations uses a plurality of different spreading codes to spread 
respective ones of said series of symbols for simultaneous transmission, such that the 
component of the received signal (X(t)) corresponding to that user station comprises a 
5 corresponding plurality of spread signals, and wherein, in said at least one of the one or 
more receiver modules (20 d ), said beamformer means (AlR dA ^. 9 41R d,Nm ) uses different 
sets of weighting coefficients to weight each element of said observation vector ( y ) to 

n 

form a plurality of signal component estimates (s dA 9 ...,s^ m ) corresponding to said 
respective ones of said series of symbols and the symbol estimating means derives from 
10 the plurality of signal component estimates (s dy \...J„' Nm ) a corresponding plurality of 
symbol estimates (b d>l ^..,b^ m ), said observation vector deriving means comprises 
means ( l^ 1 19^ m ) for deriving from the observation matrix a plurality of post- 
correlation observation vectors (Z d ' 1 ... 2^ m ) each corresponding to a respective one 

n 9 5 n 

of the plurality of different spreading codes, the channel identification means (28R d ) 
15 derives from said plurality of post-correlation observation vectors (z^ 1 Z^ m ) a 
corresponding plurality of sets of channel vector estimates ry f 1 ^™} and supplies 
the sets to said beamformer means (47R dA ,... 9 47R dyNjn ) and the coefficient tuning means 
of the beamformer means (47R d >\...,47R dyNm ) uses the sets of channel vector 
estimates (y ,...,y^ m ) to derive said different sets of weighting coefficients, 

0,rt 0,n 

20 respectively. 



32. A receiver according to claim 31, wherein the means (19* 1 ,..., 19**") for 
deriving the plurality of observation vectors comprises despreading means for 
despreading said observation matrix (YJ using one or more of said plurality of different 

25 spreading codes to form a plurality of post-correlation observation vectors 
(Z dA ,...,Z cyV| ") for u se by the channel identification means (28R d ) in deriving said 

n n 

plurality of spread channel vector estimates (y* 1 y^ m ). 

33. A receiver according to claim 32, wherein the observation vector deriving means 
30 further comprises also supplies said plurality of post-correlation observation 

vectors (Z d > l ,...,Z? yNm ) to said beamformer means (47R d > 1 47R d,Nm ) and the 
coefficient tuning means therein uses said sets of weighting coefficients to weight 
elements of respective ones of the plurality of post-correlation observation vectors, and 



139 

each of the user-specific constraint-matrix generator means comprises despreading 
means (55B dA 9 ...,55B d ' Nm ) for despreading the user-specific constraints-set matrices using 
one or more of the plurality of different spreading codes. 

5 34. A receiver according to claim 31, wherein said means for deriving said plurality 
of observation vectors comprises despreading means ( I9 rf > 6 ) for weighting said plurality 
of different spreading codes by the plurality of symbol estimates (B d ' x 9 ... 9 b^ yNm ) 9 
respectively, to form a single spreading code and despreading the observation matrix (YJ 
using said single spreading code to produce a compound post-correlation observation 
10 vector (Z a » 6 ) for use by the channel identification means (28R d ) to derive said plurality 

it 

of sets of channel vector estimates (y y^ m } - 

35. A receiver according to claim 34, wherein the observation vector deriving means 
also supplies said plurality of post-correlation observation vectors (Z d > 1 ... Z^ m ) to said 

15 beamformer means (47/f* 1 ,..., 41R dyNm ), and the coefficient tuning means therein uses 
said sets of weighting coefficients to weight elements of respective ones of the plurality 
of post-correlation observation vectors, and each of the user-specific constraint-matrix 
generator means comprises despreading means (55B dt \...,55B d,Nin ) for despreading the 
user-specific constraints-set matrices using the plurality of different spreading codes. 

20 

36. A receiver according to claim 3, wherein each of said plurality of receiver 
modules (20 1 ,...,20 NI ,20 d ) operates with a frame duration equal to integer multiples 
(F 1 ,...,F NI ,F d ) of symbol periods of the corresponding users and uses a plurality 
(F,,...,F N ,,F d ) of different segments of the same long spreading code equal to said the 

25 number of symbol periods in said frame, in said at least one of the one or more receiver 
modules (20 d ), said beamformer means (475^ 1 ,..., AlS d,Fd ) uses different sets of 
weighting coefficients to weight each element of said observation vector ( y ) to form 

n 

a plurality (F d ) of signal component estimates (§ dA ,...,5^), respectively, and the 
symbol estimating means (29S d > 1 ,...,29S d * Fd ) derives from the plurality of signal 
30 component estimates (s^' 1 ,...,^' d ) a corresponding plurality of symbol 
estimates (b% 1 ,.„ 9 B n * said observation vector deriving means deriving from the 
observation matrix one or more of a plurality of post-correlation observation 
vectors (Z dyl ,...,2?' Fd ) the channel identification means (28S d ) derives from said post- 

n n 
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correlation observation vectors (Z^ 1 ,... 2^') a corresponding plurality of spread 
channel vector estimates (y d,1 5 ... 5 j^ ,Frf ), and supplies said spread channel vector 
estimates (Y^ 1 9 ... 9 Y^ Fd ) y to said beamformer means (47S d -\... 9 AlS dJ?d ), each spread 
channel vector estimate being spread by a respective one of the segments of the long 
5 spreading code, and the coefficient tuning means of the beamformer 
means (47S d ' l 9 ... 9 47S dyFd ) uses said channel vector estimates ry^' 1 y^ ,Frf ) to derive 

V— 0,n 5 "* 9 — 0,/i ' 

said different sets of weighting coefficients. 

37. A receiver according to claim 36, wherein the means (l9S d * l 9 l9S d,F ") for 
10 deriving one or more of the plurality of post-correlation observation vectors comprises 

despreading means for despreading said observation matrix (YJ using one or more of 
said different segments of the same long spreading code such that said plurality of 
observation vectors comprise a plurality of post-correlation observation 
vectors (Z* 1 ,...,^^) one or more thereof for use by the channel identification means 
15 (28S d ). 

38. A receiver according to claim 37, wherein the observation vector deriving means 
also supplies said plurality of post-correlation observation vectors (Z dA , ...,Z?' F ") to said 

n n 

beamformer means (47S* 1 ,..., 47S d,Fd ), and the coefficient tuning means therein uses 
20 said sets of weighting coefficients to weight elements of respective ones of the plurality 
of post-correlation observation vectors, and each of the user-specific constraint-matrix 
generator means comprises despreading means {55B dA , ...,55B d,Fa ) for despreading the 
user-specific constraints-set matrices using the corresponding one of the plurality (F d ) of 
different spreading code segments. 

25 

39. A user station receiver for a CDMA communications system comprising a 
plurality (NB) of base stations (11) and a multiplicity (U) of user stations (10 1 ,...,10 u ), 
at least a plurality (U') of the user stations being in a cell associated with one of said 
base stations and served thereby, said one base station having a plurality of transmitter 

30 modules for spreading user signals for transmission to the plurality (U') of user stations, 
respectively, and a receiver for receiving spread user signals transmitted by the plurality 
(IT) of user stations, the user stations each having a receiver for receiving the 
corresponding spread user signal transmitted by the base station, said plurality (U') of 
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user stations each having a unique spreading code assigned thereto for use by the user 
station and the corresponding one of the base station transmitter modules to spread the 
user signals of that user for transmission, 

the spread user signals transmitted from the base station transmitter modules to 
5 a particular one of the plurality (U') of user stations propagating via a plurality of 
channels (M 1 ,...,^ 0 '), respectively, 

the receiver of a particular one of said plurality (U') of user stations receiving a 
signal (X(t)) comprising components corresponding to spread user signals for said 
particular user station and spread user signals transmitted by other transmitter modules 
10 of said plurality (NB) of base stations for other users, each of said spread user signals 
comprising a series of symbols spread using the spreading code associated with the 
corresponding one of the user stations, 
said user station receiver comprising: 

a plurality (NB) of receiver modules (20 v ') each for deriving from successive 
15 frames of the received signal (X(t)) estimates of sets of said series of symbols 

from a corresponding one of the base stations, 

preprocessing means (18) for deriving from the received signal (X(t)) a series of 
observation matrices (Y n ) each for use by each of the receiver modules (20 v ) in 
a said frame to derive estimates of sets of said symbols, and 
20 means (19,44) for deriving from each observation matrix a plurality of sets of 

observation vectors (Y^ 11 ... Y*' ,ni,Fn " Z"'* u z v '' Aff,Fjw ) and applying each of the 
sets of observation vectors to a respective one of the plurality of receiver modules 
(20 v '); 

each receiver modules comprising; 
25 channel identification means (28T V ) for deriving from the respective one 

of the sets of observation vectors a set of spread channel vector 
estimates (i*' 1 - 1 jr'^ 1 ^" 1 ) based upon parameter estimates of the 
channel between the corresponding one of the base stations and said user 
station; 

30 beamformer means (47r v ' 1,1 47T v '' NI ' Fm ) having coefficient tuning 

means for producing sets of weighting coefficients in dependence upon the 
sets of channel vector estimates, respectively, and combining means for 
using each of the sets of weighting coefficients to weight respective ones 
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of the elements of a respective one of the observation vectors and 
combining the weighted elements to provide a corresponding set of signal 
component estimates (s^ >l ' l 9 s„ yNIyFm ) and 

symbol estimating means ( 29r v/ ' u 29T v ' tNltFm ) for deriving from the 
5 set of signal component estimates a set of estimates ( M , £^ /,Fw ) of 

symbols spread by the corresponding one of the transmitter modules and 
transmitted by the base station; 
said user station receiver further comprising means (42,43) responsive to said 
symbol estimates (B v '' u \... 9 6 v, ^ M iz\i? 9 £) and channel 
10 estimates (IK*') from each of said plurality (NB) of receiver modules, said 

channel estimates comprising at least channel vector estimates (H v ') for channels 

n 

(14 v ) between the user station receiver and said base stations, for providing at 
least one constraint matrix (C n ) representing interference subspace of 
components of the received signal corresponding to said spread signals, and in 
15 each of said receiver modules (20 v ), the coefficient tuning means produces said 

sets of weighting coefficients in dependence upon both the constraint matrix 
(C n ) and the channel vector estimates so as to tune said receiver module (20 v ) 
towards a substantially null response to that portion of the received signal (X(t)) 
corresponding to said interference subspace. 

20 

40. A user station receiver according to claim 39, wherein said observation vector 
deriving means comprises means ( igv'.u ^ i^'^^ for deriving from the observation 
matrix a plurality of post-correlation observation vectors (Z v/ » 1 ' 1 ? ... 5 Z V ,yNI * Fm ) and 

n n 

supplying said plurality of post-correlation observation vectors (Z"'' 1,1 ... z v/yNIJ?NI ) to the 
25 channel identification means (28T V ) for use in producing said sets of channel vector 
estimates (y^u ^'^^ . 

41. A user station receiver according to claim 39, wherein said plurality (NB) of 
receiver modules (20 v ) derive symbols for user signals other than those destined for the 

30 user of said user station receiver, and said user station receiver comprises an additional 
receiver module (20 d ) for deriving symbols from the received signal destined for said 
user of said user station receiver and transmitted by a corresponding serving one (v) of 
said plurality of base stations, wherein each of said plurality of receiver modules 
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(20 v ',20 d ) operates with a frame duration equal to an integer multiple of the symbol 
period and uses a corresponding number of segments (F 1 ,...,F NI ,F d ) of a long spreading 
code, each segment corresponding to a respective one of a plurality (F l5 .. .,F NI ,F d ) of 
different segments of the same long spreading code equal to said plurality of symbol 
5 periods in said frame, and wherein, in said additional receiver module (20 d ), said 
beamformer means (47S* 1 ,..., 47S d ' Fd ) uses different sets of weighting coefficients to 
weight each element of said observation vector ( F ) to form a plurality (F d ) of signal 

n 

component estimates (s d '\...J*' Fd ), respectively, and the symbol estimating means 
(29S d > x ,...,29S d ' Fd ) derives from the plurality of signal component estimates 
10 (s„ l ,...J„ Fd ) a corresponding plurality of symbol estimates (b dA ,...,b* ,Fd ) 9 said 
observation vector deriving means deriving from the observation matrix a plurality of 
observation vectors ( y** 1 Y 4 ^) the channel identification means (28S d ) derives from 

n n 

said observation vectors (Y dfl i ....Y dyFd ) a corresponding plurality of spread channel 

n n 

vector estimates (f ,..,y d,Frf> ) and supplies same said beamformer 
15 means (47S d * 1 ,..., 47S d,Fd ), each spread channel vector estimate being spread by a 
respective one of the segments of the long spreading code, and the coefficient tuning 
means of the beamformer means (47S d * l ,... 9 47S dtFd ) uses said channel vector 
estimates (y^ 1 v .. 5 Y^ Frf ) to derive said different sets of weighting coefficients. 

20 42. A user station receiver according to claim 39, wherein said plurality (NB) of 
receiver modules (20 v ) derive symbols for user signals other than those destined for the 
user of said user station receiver, and said user station receiver comprises an additional 
receiver module (20 d ) for deriving symbols from the received signal destined for said 
user of said user station receiver, wherein each of said plurality of receiver modules 

25 (20 v ',20 d ) operates with a frame duration equal to an integer multiple of the symbol 
period and uses a corresponding number of segments of a long spreading code, each 
segment corresponding to a respective one of a plurality (F l5 ...,F NI ,F d ) of different 
segments of the same long spreading code equal to said plurality of symbol periods in 
said frame, and wherein, in said additional receiver module (20 d ), said beamformer 

30 means (47r v, ^ 1 v .. ? 47T VydyFd ) uses different sets of weighting coefficients to weight each 
element of said observation vector ( y ) to form a plurality (F d ) of signal component 

n 

estimates (s„ d>l ' d ), respectively, and the symbol estimating means 

(29T v > dA ,... 9 29T v ' dyFd ) derives from the plurality of signal component estimates 
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(Sn**' 1 > ->K dFd ) a corresponding plurality of symbol estimates (b^ d ' x 9 ...,b^ d ' Fd ), and 
wherein the coefficient tuning means of said beamformer means 

(47T vA1 ,...,47T vAFd ) derives the weighting coefficients using said constraint matrix 
received from the constraint matrix generating means (43T) and said spread channel 
5 vector estimates (y^ 1 v .. 5 y 1 ^) produced by the channel identification means (28T V ) 

0,n 0,n 

of the receiver module (20 v ) corresponding to the base station (v) serving said user 
station. 

43. A user station receiver according to claim 42, wherein said observation vector 
10 deriving means comprises despreading means ( i9 v '>u ig v/ ^ F Ni^ f or despreading the 

observation matrix to form a plurality of post-correlation observation 
vectors (Z^ 1 * 1 ,„.,Z V yNI ' FNI ) and supplying said plurality of post-correlation observation 

n n 

vectors (Z v/>14 3 ... 5 Z v,w,Fm ) to the channel identification means (28T V ) for use in 

n n 

producing said sets of channel vector estimates ry^'M v*'^ 7 ^ . 

15 

44. A user station receiver according to claim 39, wherein at least one said 
transmitter module of one of the base stations uses a plurality (N m ) of different spreading 
codes to spread respective ones of said series of symbols for simultaneous transmission 
in the same frame, such that the component of the received signal (X(t)) corresponding 

20 to that base station transmitter module comprises a corresponding plurality of spread 
signals, and at least one (20 v ) of the plurality of receiver modules (20U) further 
comprises amplitude estimation means (30U V ) for deriving total amplitude of a set of 
signal component estimates '' u * WVm ) produced by beamformer means 
(47t/ v ' M v .., 47U V /yNIyNm ) thereof, said beamformer means (47£/ v /,u v .., 47f/ v /yNIiNm ) uses 

25 different sets of weighting coefficients to weight each element of said observation 
vector (y ) to form said plurality of signal component 
estimates ,1,1 ? ... 5 i^ ^ 7 ^ m ) corresponding to said respective ones of said series of 
symbols, and the symbol estimating means (29t/ v ' ,u 29U V fyNI ^ m ) derives from the 
plurality of signal component estimates (s^ xx ,WV/n ) a corresponding plurality of 

30 symbol estimates {b v n Ail ,...,b^ IJ4m ), the channel identification means (28U V ) derives 
a corresponding plurality of sets of spread channel vector 
estimates (y v/ * 1,1 5 ... 5 y" ^ 7 ^ m ) each spread by a respective one of said plurality of 
different spreading codes, and supplies the sets to said beamformer 



145 

means (47U y ' x \..., 41 U v '- NI ' N -) and the coefficient tuning means of the beamformer 
means (41U v '' l '\...,41U v '' N '' Nm ) uses the sets of spread channel vector 
estimates {yy'' 1 ' 1 ,...,f ,yNI ' Nm ) to derive said different sets of weighting coefficients, 

respectively. 

5 

45. A user station receiver according to claim 44, wherein said plurality (NB) of 
receiver modules (20 v ) derive symbols for user signals other than those destined for the 
user of said user station receiver, and said user station receiver comprises an additional 
receiver module (20 d ) for deriving symbols from the received signal destined for said 

10 user of said user station receiver, wherein each of said plurality of receiver modules 
(20 v ',20 d ) operates with a frame duration equal to an integer multiple of the symbol 
period and uses a corresponding number of segments of a long spreading code, each 
segment corresponding to a respective one of a plurality (N m ,F d ) of different segments 
of the same long spreading code equal to said plurality of symbol periods in said frame, 

15 and wherein, in said additional receiver module (20 d ), said beamformer 
means (47T v > d > 1 ,..., 47T v ' d ' Fd ) uses different sets of weighting coefficients to weight each 
element of said observation vector (y ) to form a plurality (F d ) of signal component 

n 

estimates (§l Al w ), respectively, and the symbol estimating means 

(29T v > d > 1 ,... y 29T vAFd ) derives from the plurality of signal component estimates 
20 (s* Al ,s^' d,Fd ) a corresponding plurality of symbol estimates ( b vAl , b^ dyFd ), and 
wherein the coefficient tuning means of said beamformer means 
(47T v > d > l ,.„, 47T vAFd ) derives the weighting coefficients using said constraint matrix 
received from the constraint matrix generating means (43T) and said spread channel 
vector estimates (f* Al Y" AFd ) produced by the channel identification means (28T V ) 
25 of the receiver module (20 v ) corresponding to the base station (v) serving said user 
station. 

46. A receiver according to claim 44, wherein said observation vector deriving means 
comprises despreading means (19 V ' U ,..., 19 v for despreading the observation 

30 matrix to form a plurality of post-correlation observation vectors (Z^' 1 ' 1 ... Z v /yNIyNm ) and 
supplying said plurality of post-correlation observation vectors (Z^' 1 ' 1 ... Z v /fNIyNm ) to the 
channel identification means (28U V ) for use in producing said sets of channel vector 
estimates (y^'u ij?'J*Wm\ . 
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47. A user station receiver according to claim 39, wherein at least one said 
transmitter module of one of the base stations uses a plurality (N m ) of different spreading 
codes to spread respective ones of said series of symbols for simultaneous transmission 
in the same frame, such that the component of the received signal (X(t)) corresponding 

5 to that base station transmitter module comprises a corresponding plurality of spread 
signals, and at least one (20U') of the plurality of receiver modules further comprises 
amplitude estimation means (30U V ) for deriving total amplitudes 9 _ 9 ^' NI ) of a 

set of signal component estimates (s^' 1 ' 1 9 ... 9 s^' NI ' Nm ) produced by beamformer 
means (47t/ v '' u v .., 47f/ v ^ /JVffl ) thereof, said beamformer 

10 means (47t/ v ' u ,..., 47 U v /yNIfNm ) uses different sets of weighting coefficients to weight 
each element of said observation vector ( y ) to form said plurality of signal component 
estimates (s*'' 1 ' 1 ' NI ' Nm ) corresponding to said respective ones of said series of 
symbols and the symbol estimating means (29t/ v ' u 29U V /yNIyNm ) derives from the 
plurality of signal component estimates (s„ yl,x ,...J*'' NI ' Nm ) a corresponding plurality of 

15 symbol estimates ( b„' lyX , b„ said observation vector deriving means comprises 

despreading means (19 v/s,1 ,... 5 19 v,s,iVm ) for despreading the observation matrix (Y n ), 
using one or more of said plurality (N m ) of different spreading codes, each of these codes 
being a compound code formed by averaging the codes of all others of said multiplicity 
of interfering users (NI), to form a plurality of observation vectors (Z^* 2 * 1 V ..,Z V ,s ^ m ) , 

n n 

20 the channel identification means (28U V ) derives from said plurality of observation 

vectors (Z v/ ' SJ Z v,s ^ m ) , said signal component estimates ri*' 2 * 1 s^^"") and 

said amplitudes (ifr*' 1 v ..,iJj^ a/ ) a corresponding plurality of sets of channel vector 
estimates (Y v/ > 1 >\„.^ / ' NIyNm ) and supplies the sets to said beamformer 
means (47l/ v ' 1 ' 1 ,..., 47U V yNIyNm ) and the coefficient tuning means of the beamformer 

25 means (41U vfxx ,„.,41U VyNIyNm ) uses the sets of channel vector 
estimates (T^- 1 ' 1 ,..., T / ' NI ' Nm ) to derive said different sets of weighting coefficients, 
respectively. 

48. A receiver according to claim 1, each said receiver module it located in a 
30 user/mobile station and the received signal comprises a plurality of spread user signals 

transmitted by a plurality of transmitter modules at a base station communicating with 
said receiver via said channels. 
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49. A receiver for a CDMA communications system comprising preprocessing means 
(18) for deriving from a received signal (X(t)) a succession of observation matrices (Y n ), 
means (44) for reshaping each observation matrix to form an observation vector (Y n )> 
beamformer means (27N d ) for weighting each element of the observation vector using 
5 weighting coefficients and combining the weighted elements to form a signal component 
estimate ($^) , said beamformer means comprising tuning means for tuning said 

weighting coefficients in dependence upon a channel vector estimate (y 4 ^ , means 

v — ox 

(29N d ) for deriving from the signal component estimate a corresponding symbol 

estimate for output from the receiver module, means (19 d ) for despreading the 

10 observation matrix using the spreading code of the corresponding user to form a post- 
correlation observation vector (z d ) , and channel identification means (28N d ) for 
deriving said channel vector estimate ) in dependence upon said post-correlation 

0>n 

observation vector (z rf ) and said signal component estimate (§ d ) . 

15 50. A CDMA communications system comprising at least one base station and a 
multiplicity (U) of user stations (10\ 10 u ) including a plurality (U') of user stations 
served by said at least one base station, the base station capable of transmitting unique 
space-time encoded signals each dedicated to a corresponding mobile user station, each 
user station having a transmitter and a receiver for communicating with said at least one 
20 base station via a corresponding one of a plurality of channels (14 1 ,...,14 u ), the base 
station transmitter comprising: 

a plurality of transmission antennas; 
means for providing user-destined signals; 

a distribution unit for grouping the user-destined signals into N groups; 
25 a temporal channelization-code unit for spreading each user-destined signal by a 

unique dedicated code belonging to a fixed set of L orthogonal codes and summing the 
spread signals in each group; 

means for adding to the summed spread signals of each group a respective one 
of a plurality of pilot signals each specific to one of the groups and generated by a PN 
30 code generator; 

means for scrambling the summed signal from each of said groups using the same 
long scrambling code specific to the base station; 
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each pilot signal being assigned a fixed fraction of the total power transmitted 
from the base station transmitter; 

means for mapping the signals from the groups (G.(r), -, G N (r)) onto antenna 

branches (A x (f), -,A M (t)) by means of a linear space coding (M) such that signals 
5 assigned to different groups are substantially orthogonal at transmission; 

at least one of the user stations having a receiver for receiving the corresponding 
spread user signal transmitted by the base station, said plurality (U') of user stations each 
having a unique spreading code assigned thereto for use by the user station and the 
corresponding one of the base station transmitter modules to spread the user signals of 
10 that user for transmission, 

the spread user signals transmitted from the base station transmitter modules to 
a particular one of the plurality (U') of user stations propagating via a plurality of 
channels (14\...,14 u ), respectively, 

the receiver of a particular one of said plurality (U') of user stations receiving a 
15 signal (X(t)) comprising components corresponding to spread user signals for said 
particular user station and spread user signals transmitted by other transmitter modules 
of said plurality (NB) of base stations for other users, each of said spread user signals 
comprising a series of symbols spread using the spreading code associated with the 
corresponding one of the user stations, 
20 said user station receiver comprising: 

a plurality (NB) of receiver modules (20 v ) each for deriving from successive 

frames of the received signal (X(t)) estimates of sets of said series of symbols 

from a corresponding one of the base stations, 

preprocessing means (18) for deriving from the received signal (X(t)) a series of 
25 observation matrices (F„) each for use by each of the receiver modules (20 v ) in 

a said frame to derive estimates of sets of said symbols, and 
means (19,44) for deriving from each observation matrix a plurality of sets of 
observation vectors (Y v ' tUl 9 ... > J^' JV/,F «; Z"' ,u , ... j Z v ' Ar/jFw ) and applying each of the 
sets of observation vectors to a respective one of the plurality of receiver modules 
30 (20 v ); 

each receiver module comprising; 

channel identification means (28T V ) for deriving from the respective one 
of the sets of observation vectors a set of spread channel vector 
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estimates (jr' 1,1 f '' ni ' Fn ') based upon parameter estimates of the 
channel between the corresponding one of the base stations and said user 
station; 

beamformer means (47r v ' M 47T V /fNIlFNf ) having coefficient tuning 
5 means for producing sets of weighting coefficients in dependence upon the 

sets of channel vector estimates, respectively, and combining means for 
using each of the sets of weighting coefficients to weight respective ones 
of the elements of a respective one of the observation vectors and 
combining the weighted elements to provide a corresponding set of signal 

10 component estimates (5^ ^ /,Fw ) and 

symbol estimating means ( 29r v ' u 29T vfyNIyFm ) for deriving from the 
set of signal component estimates a set of estimates ( b„' l,1 9 /yNI ^ NI ) of 
symbols spread by the corresponding one of the transmitter modules and 
transmitted by the base station; 

15 said user station receiver further comprising means (42,43) responsive to said 

symbol estimates (i v ' u ,...,£ V /,WF/ "; gKg 1 ,? 3 ) and channel 

f n n n n n 

estimates ( < H v n ) from each of said plurality (NB) of receiver modules, said 
channel estimates comprising at least channel vector estimates (H v ) for channels 

n 

(14 v ) between the user station receiver and said base stations, for providing at 
20 least one constraint matrix (C n ) representing interference subspace of 

components of the received signal corresponding to said spread signals, and in 
each of said receiver modules (20 v ), the coefficient tuning means produces said 
sets of weighting coefficients in dependence upon both the constraint matrix 
(C n ) and the channel vector estimates so as to tune said receiver module (20 v ) 
25 towards a substantially null response to that portion of the received signal (X(t)) 

corresponding to said interference subspace. 



51. A CDMA system according to claim 50, the transmitter further comprising: 

delay means for delaying the signals from the mapping means each by a branch- 
30 specific delay, forming a corresponding signalling pulse modulating the signalling pulse 
with a carrier frequency signal, and supplying the modulated signal to the antenna 
elements for transmission thereby. 
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52. A CDMA system according to claim 50, wherein, in the transmitter, the grouping 
means is arranged to group data signals dedicated for user stations into a predetermined 
number N G of groups and, when the number of users is less than or equal to the 
processing gain, assigning all user signals to the same group; otherwise assigning user 

5 signals pseudo-randomly so as to tend to balance the number of users in each group. 

53. A CDMA system according to claim 50, wherein the transmitter further 
comprises, for each group, a space-time coding means comprising: 

a channel coding unit for spreading user signals by a group specific set of 
10 orthogonal codes; 

a scrambling unit for spreading each group of user signals with the same base 
station specific scrambling code; and 

a spatial coding unit for mapping the total group signals (G (r), G (0) onto 
the antenna branch signals (A At), -,A„ (t)) by a linear transformation (M). 

1 N c 

15 

54. A CDMA system according to claim 53, wherein, in the transmitter, the channel 
coding unit assigns to each user in the group indexed g a chip-code chosen from a fixed 
set of L g < L orthogonal group specific L-chip codes ft), - 9 a iM (t)). , the code-sets 
across groups being chosen subject to minimize the maximum cross-correlation, any code 

20 belonging to any group being orthogonal to any other code within the group while its 
cross-correlation with any out-group channelization code is minimized. 

55. A CDMA communications system comprising at least one base station and a 
multiplicity (U) of user stations (10 1 , 10 u ) including a plurality (U') of user stations 

25 served by said at least one base station, each user station having a transmitter and a 
receiver for communicating with said at least one base station via a corresponding one 
of a plurality of channels (M 1 ,...,^), at least one user station being capable of 
transmitting a user signal comprising a plurality of unique space-time encoded signals 
each carrying different data from that same user, 
30 said at least one user station having a transmitter comprising: 
a plurality of transmission antennas; 
means for providing said user signals; 

a distribution unit for grouping the space-time encoded signals into N groups; 
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a temporal channelization-code unit for spreading each different data stream of 
the user signals by a unique dedicated code belonging to a fixed set of L orthogonal 
codes and summing the spread signals in each group; 

means for adding to the summed spread signals of each group a respective one 
5 of a plurality of pilot signals each specific to one of the groups and generated by a PN 
code generator; 

each pilot signal being assigned a fixed fraction of the total power transmitted 
from the base station transmitter; 

means for scrambling the summed signals from the groups using the same long 
10 scrambling code specific to the user station; 

means for mapping the signals from the groups (G (f), G (r)) onto antenna 
branches (A 1 (f) 9 -,A M (t)) by means of a linear space coding (M) such that signals 
assigned to different groups are substantially orthogonal at transmission; 

the base station having a receiver for receiving a signal (X(t)) comprising 
15 components corresponding to spread signals transmitted by the transmitters of the 
plurality of user stations, each of said spread signals comprising a series of symbols 
spread using a spreading code unique to the corresponding user station, said base station 
receiver comprising: 

a plurality (U') of receiver modules (20 1 ,...,20 NI ,20 d ) each for deriving from 
20 successive frames of the received signal (X(t)) estimates of said series of symbols 

of a corresponding one of the user stations, 

preprocessing means (18) for deriving from the received signal (X(t)) a series of 
observation matrices (Y n ) each for use by each of the receiver modules (20) in a 
said frame to derive an estimate of a symbol of a respective one of said series of 
25 symbols, and 

means (19,44;44/l,44/2) for deriving from each observation matrix a plurality of 
observation vectors (Y ; Y \Z x ..Z NIt Z**) and applying each of the observation 
vectors to a respective one of the plurality of receiver modules (20 1 ,...,20 NI ,20 d ); 
each receiver module comprising; 
30 channel identification means (28) for deriving from one of the observation 

vectors a channel vector estimate (i/ 1 ,...,//^ 7 ;^ ;Y* ) based upon 
parameter estimates of the channel between the base station receiver and 
the corresponding user station transmitter; 
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bearnformer means (27\ . . . ,27 NI ,27 d ; 47 d ) having coefficient tuning means 
(50) for producing a set of weighting coefficients in dependence upon the 
channel vector estimate, and combining means (51,52) for using the 
weighting coefficients to weight respective ones of the elements of a 
5 respective one of the observation vectors and combining the weighted 

elements to provide a signal component estimate ($,),...,£ ^ and 
symbol estimating means (29 1 ,...,29 u , 30 1 ,...,30 u ) for deriving from the 
signal component estimate an estimate (b l ... b u ) of a 
symbol (b^... J b^ / ) transmitted by a corresponding one of the user 
10 stations (10\...,10 u ), 

wherein said receiver further comprises means (42,43) responsive to symbol 
estimates (b^b^ig 1 ^ 2 ^ 3 ^ 1 ^) and to channel estimates 
(Q-C\ ... ( H N n I \ ttl-i) comprising at least said channel vector 
estimates (H l ,...,H NI ) for channels (14 I ,...,14 NI ) of a first group (I) of said 

n n 

15 plurality of user stations (10\...,10 NI ) to provide at least one constraint 

matrix (£ ) representing interference subspace of components of the received 
signal corresponding to said predetermined group, and in each of one or more 
receiver modules (20A d ) of a second group (D) of said plurality of receiver 
modules, the coefficient tuning means (50 A d ) produces said set of weighting 

20 coefficients in dependence upon both the constraint matrix ( c ) and the channel 

vector estimates Q?*) so as to tune said one or more receiver modules (20A d ) 
each towards a substantially null response to that portion of the received signal 
(X(t)) corresponding to said interference subspace. 



25 

56. A CDMA system according to claim 55, the transmitter further comprising: 
delay means for delaying the signals from the mapping means each by a branch- 
specific delay, forming a corresponding signalling pulse modulating the signalling pulse 
with a carrier frequency signal, and supplying the modulated signal to the antenna 

30 elements for transmission thereby. 

57. A CDMA system according to claim 55, wherein, in the transmitter, the 
distribution unit means is arranged to group data signals dedicated for user stations into 
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a predetermined number # of groups and, when the number of users is less than or 
equal to the processing gain, assign all user signals to the same group, otherwise 
assigning user signals pseudo-randomly so as to tend to balance the number of users in 
each group. 

5 

58. A CDMA system according to claim 55, wherein the transmitter further 
comprises, for each group, space-time coding means comprising: 

a channel coding unit for spreading user signals by a group specific set of 
orthogonal codes; 

10 a scrambling unit for spreading each group of user signals with the same base 

station specific scrambling code; and 

a spatial coding unit for mapping the total group signals (G (r),- ■■, G (t)) onto 
the antenna branch signals (AM), -,A N (r)) by a linear transformation (M). 

15 59. A CDMA system according to claim 58, wherein, in the transmitter, the channel 
coding unit assigns to each user in the group indexed g a chip-code chosen from a fixed 
setof l < l orthogonal group specific L-chip codes ( a . x (f) 9 - ,a iL (f)). , the code-sets 
across groups being chosen subject to minimize the maximum cross-correlation, any code 
belonging to any group being orthogonal to any other code within the group while its 

20 cross-correlation with any out-group channelization code is minimized. 

60. A receiver according to claim 1, for use with a transmitter transmitting pilot- 
symbol assisted user signals comprising pilot symbols multiplexed with data symbols, the 
receiver further comprising demultiplexing means (35V d ) for demultiplexing the signal 

25 component estimates from the ISR beamformer to extract pilot signal component 
estimates and data signal component estimates and supplying the data signal component 
estimates to the decision rule unit (29V d ) and the pilot signal component estimates to an 
ambiguity estimation means (31V d ), the ambiguity estimation means (31V d ) smoothing 
or averaging each pilot signal component estimate (5™^) to provide an ambiguity 

30 estimate (^) , conjugation means (32V d ) for deriving from the ambiguity 
estimate (a d ) its conjugate ((a*)*) and multiplier means (15 V d ) for multiplying the 
conjugate with the symbol estimate from the decision rule means (29V d ) to form 

an improved symbol estimate Q> d ) . 



A 
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61. A receiver according to claim 1, further comprising a second ISR beamformer 
(47V/2 d ) connected in parallel with the first ISR beamformer (47V/ l d ) and responsive to 
the same channel coefficients and constraint matrix as the first ISR beamformer to derive 
from the observation vector ( y ) a pilot signal component estimate, ambiguity 

5 estimation means (31V d ) for smoothing or averaging the pilot signal component 
estimate (§*' d ) to provide an ambiguity estimate (d d ) , conjugation means (32 V d ) for 
deriving from the ambiguity estimate its conjugate ((a n V) and multiplier means 

(15V d ) for multiplying the conjugate with the symbol estimate (f>*) from the decision 
rule unit (29V d ) to form an improved symbol estimate (b d ) . 

10 

62. A receiver according to claim 60, wherein the ambiguity estimation means 
comprises a buffer for buffering bit estimates, smoothing means for smoothing or 
averaging the buffered estimates and a further decision rule unit (29V/2 d ) for deriving 
from the smoothed or averaged estimates said corresponding ambiguity estimates (d„) - 

15 

63. A receiver according to claim 61, wherein the ambiguity estimation means 
comprises a buffer for buffering bit estimates, smoothing means for smoothing or 
averaging the buffered estimates and a further decision rule unit (29V/2 d ) for deriving 
from the smoothed or averaged estimates said corresponding ambiguity estimates (d d ) • 

20 

64. A receiver according to claim 39, for use with a transmitter transmitting pilot- 
symbol assisted user signals comprising pilot symbols multiplexed with data symbols, the 
receiver further comprising demultiplexing means (35V d ) for demultiplexing the signal 
component estimates from the ISR beamformer to extract pilot signal component 

25 estimates and data signal component estimates and supplying the data signal component 
estimates to the decision rule unit (29V d ) and the pilot signal component estimates to an 
ambiguity estimation means (31V d ), the ambiguity estimation means (31V d ) smoothing 
or averaging each pilot signal component estimate (sf 4 ) to provide an ambiguity 
estimate (d d ) , conjugation means (32V d ) for deriving from the ambiguity 

30 estimate its conjugate and multiplier means (15V d ) for multiplying the 

conjugate with the symbol estimate (£ <*) from the decision rule means (29V d ) to form 
an improved symbol estimate ($ d ) . 
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65. A receiver according to claim 39, further comprising a second ISR beamformer 
(47V/2 d ) connected in parallel with the first ISR beamformer (47V/ l d ) and responsive to 
the same channel coefficients and constraint matrix as the first ISR beamformer to derive 
from the observation vector ( 7 ) a pilot signal component estimate , ambiguity 

n 

5 estimation means (31V d ) for smoothing or averaging the pilot signal component 
estimate (s* 4 ) to provide an ambiguity estimate (d*) , conjugation means (32V d ) for 
deriving from the ambiguity estimate (d^) its conjugate ((d^)*) and multiplier means 
(15V d ) for multiplying the conjugate with the symbol estimate (S*) from the decision 
rule unit (29V d ) to form an improved symbol estimate (& d ) - 

n 

10 

66. A receiver according to claim 64, wherein the ambiguity estimation means 
comprises a buffer for buffering bit estimates, smoothing means for smoothing or 
averaging the buffered estimates and a further decision rule unit (29V/2 d ) for deriving 
from the smoothed or averaged estimates said corresponding ambiguity estimates (d„) • 

15 

67. A receiver according to claim 65, wherein the ambiguity estimation means 
comprises a buffer for buffering bit estimates, smoothing means for smoothing or 
averaging the buffered estimates and a further decision rule unit (29V/2 d ) for deriving 
from the smoothed or averaged estimates said corresponding ambiguity estimates (d£) • 

20 

68. A receiver according to claim 19, further comprising buffering means (90*), 
deinterleaving means (91% channel decoding means (92 l ), re-encoding means (93 1 ) and 
re-interleaving means (94 1 ) for buffering a first frame of signal component estimates 

(s„_ x (l)) , deinterleaving, channel decoding, re-encoding and re-interleaving the frame 
25' of signal component estimates to provide a frame of improved decided symbol estimates 
(£j) , and supplying same to the constraint set generator (42W), the constraint matrix 
generator (43W) forming therefrom an improved constraint matrix, the ISR beamformer 
(47W 1 ) using said improved constraint matrix to provide an improved signal component 
estimate (£^(2)) for use in a next iteration, the receiver repeating the iterations a 
30 predetermined number of times. 

69. A receiver according to claim 39, further comprising buffering means (90'), 
deinterleaving means (91'), channel decoding means (92 1 ), re-encoding means (93*) and 
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re-interleaving means (94*) for buffering a first frame of signal component estimates 
CVi(l)) > deinterleaving, channel decoding, re-encoding and re-interleaving the frame 

of signal component estimates to provide a frame of improved decided symbol estimates 
(b^) , and supplying same to the constraint set generator (42 W), the constraint matrix 
5 generator (43 W) forming therefrom an improved constraint matrix, the ISR beamformer 

(47W 1 ) using said improved constraint matrix to provide an improved signal component 

estimate (s*_ x (2y) f° r use in a nex t iteration, the receiver repeating the iterations a 

predetermined number of times. 

10 70. A receiver according to claim 1, wherein coefficient tuning means (50A d ) 
produces said set of weighting coefficients in dependence upon both the constraint matrix 
(C ) and the channel vector estimates rfj d ) so as to tune said one or more receiver 

n — n 

modules (20A d ) each towards a response that is bound in magnitude to be close but not 
equal to a null. 

15 

71. A receiver according to claim 39, wherein coefficient tuning means (50A d ) 
produces said set of weighting coefficients in dependence upon both the constraint matrix 
(C ) and the channel vector estimates (H d ) so as to tune said one or more receiver 

n — n 

modules (20A d ) each towards a response that is bound in magnitude to be close but not 
20 equal to a null. 
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ABSTRACT OF THE DISCLOSURE 

A receiver of the present invention addresses the need for improved interference 
5 suppression without the number of transmissions by the power control system being 
increased, and, to this end, provides a receiver for a CDMA communications system 
which employs interference subspace rejection to tune a substantially null response to 
interference components from selected signals of other user stations. Preferably, the 
receiver also tunes a substantially unity response for a propagation channel via which a 
10 corresponding user's signal was received. The receiver may be used in a base station 
or in a user/mobile station. 



